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ABSTRACT 


The 6 MV and 15 MV photon beams produced by linear accelerators 
were studied for electron and scattered photon contamination. Tissue 
maximum ratios with and without accessories in the field indicate the 
presence of contamination. The surface dose was found to be 
attributable to contamination electrons and backscatter photons. The 
surface dose produced by accessories was reduced when electron filters 
made of materials with high atomic numbers were placed underneath the 
accessory. The surface dose for rectangular fields was measured and 
the equivalent squares for the surface dose were derived from these 
measurements. The 15 MV accelerator had a Gaussian lateral distribu- 
Gion, a linear dependence on field width for square fields and an 
inverse square dependence on distance from the bottom of the fixed 
head assembly. This geometrical dependence is consistent with the 
proposal that the field flattening filter is the main source of 
electron contamination at large field sizes when accessories are 
present. The electron contamination at 6 MV could not be localized to 
the same location. A permanent magnet was used to sweep electrons from 
Pies field ae lpavenerd ly. .cl ee tron contamination measured by a probe was 
produced by material close to the probe. A penetration curve for 
electron and scattered photon contamination was produced utilizing the 
linearity of dose with respect to field width. The derived contamina- 
tion curve was similar to the measured build-up curve outside the 
field. The phantom generated tissue maximum ratio, obtained by sub- 


tracting the contamination contributions showed no dependence on 
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field size, source-to-probe distance or presence of accessories. 

The phantom generated tissue maximum ratio was verified independently 
by Monte Carlo simulation. The clinical aspects of the study and the 
requirements for an electromagnet to sweep a large field completely 
free of contaminant electrons were discussed. 

Volume 2 is introduced with a review of dose computation methods 
for radiotherapy. Several deficiencies were revealed, of which, the 
most important was a failure to account for the dose in situations of 
electronic disequilibrium. It was shown that the reduced central axis 
dose in a low-density heterogeneous region was due to electronic dis- 
equilibrium. It was also shown that electronic equilibrium never 
exists near the beam boundary and this gives rise to the beam penumbra 
for high energy linear accelerator beams. 

A dose calculation method based on convolution was introduced 
which accounts for both the transport of charged particles and 
scattered photons. Arrays were generated using the Monte Carlo method 
representing the energy absorbed throughout water-like phantoms from 
charged particles and scattered radiation set in motion by primary x- 
ray interactions at one location. The resulting ''dose spread arrays'' 
were normalized to the collision fraction of the kinetic energy 
released by the primary x-rays. These arrays are convolved with the 
relative primary fluence interacting in a phantom to obtain 3-dimen- 
sional dose distributions. The method gives good agreement for the 
dose in electronic disequilibrium situations such as the build-up 
region, near beam boundaries and near low-density heterogeneities 


irradiated by beams with small field sizes. 
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Preface 


Volume 1 entitled, "Contamination of Megavoltage X-Ray 
Beams By Electrons And Scattered Photons", was originally 
written as a thesis for the partial fulfilment of a Master 
of Science degree. At the defense of the thesis in 
September of 1982, it was felt by the examining committee 
that Volume 1 should be used towards a Doctor of Philosophy 
degree. 

Section 3, in Volume 1, addressed the problem of the 
transport of charged particles in the build-up region of a 
unit-density phantom. The committee felt that this Section 
provided a good foundation for a more complete study of the 
transport of charged particles generated by photon beams in 
an inhomogeneous medium. A decision was made to undertake 
this investigation to complete the thesis. Volume 2 
contains the results of this study as well as a method to 
calculate the primary and scattered dose in heterogeneous 
media. 

It was decided to break the thesis into two Volumes 
because Volume 1 is a self-contained study of photon beam 
contamination.» the, contents of Volume 2;are not) referenced 
in Volume 1. For this reason, the Discussion, Conclusions 
and Appendices concerning contamination are left in Volume 
ig However, in a number of instances, the reader of Volume 
2 is referred to Volume 1, therefore, the pagination in 


Volume 2 is a continuation of Volume l. 
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VOLUME 1 


CONTAMINATION OF MEGAVOLTAGE X-RAY BEAMS 


BY ELECTRONS AND SCATTERED PHOTONS 
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1. INTRODUCTION 


The errors of definitions 
multiply themselves accordingly 
as the reckoning proceeds; 

and lead men into absurdities, 
which at last they cannot avoid, 
without reckoning anew 

from the beginning. 


Thomas Hobbs 


1.1 Rationale and Introduction to the Study 

Megavoltage photon beams are the radiation modality most often 
used in the treatment of cancer by radiotherapy. Photons produced 
by medical linear accelerators have become the most important method 
of delivery of these beams. * 

Megavoltage photon beams become contaminated with relativistic 
electrons and scattered photons. An understanding of the role of 
electron and scattered photon contamination is required to characterize 
dosimetry in the build-up region. Clinicaily, high energy x-ray build-up 
produces a skin-sparing effect and contamination tends to reduce this 
effect. 

The Cross Cancer Institute in Edmonton, Alberta employs three 
Siemens linear accelerators, two of which operate at 6 MV photons and one 
which delivers a 15 MV photon beam. These accelerators have a very 
Similar design so an opportunity existed to compare and contrast the 
contamination at two different nominal energies. 

The dependence of position inside and outside the photon beam, 
distance from the source of primary photons and cross-sectional 
dimensions (field size) of the beam on the nature and amount of contamina- 
tion was investigated. Accessories were placed in the path of the beam to 
determine their effect on contamination. 

The photon beam was passed through the poles of a permanent magnet 
to sweep the contamination electrons out of the beam. The depth of 
penetration of the contamination electrons was determined by comparing 
dose measurements at depth in a phantom with and without the magnet in 
place. 


* Cobalt-60 isotope sources and betatrons can deliver megavoltage photon 
beams as well. 
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The depth of penetration for the electron and scattered photon 
components was determined utilizing the increase in dose with respect 
to the field size of the beam. The tissue maximum ratio (TMR) in the 
build-up region of the phantom generated component was obtained by sub- 
traction of the contaminant contribution. 

A knowledge of the nature of the primary photon component of the 
beam is important to the clinical use of such a photon beam because most 
treatment planning systems isolate the primary from the scattered photon 
component of dose. The validity of the phantom generated build-up curve, 
as obtained by subtraction of the contamination fraction, was verified 
by computer modelling the primary photon contribution. 

The nature and sources of electron and scattered photon contamina- 
tion, based on the experimental and calculated results, is discussed in 
detail. 

The role in clinical radiotherapy of electromagnets to sweep 
electron contamination from the field, and electron "filters" to scatter 


electron contamination, is discussed briefly. 
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1.2 Production and Interaction of X-Rays with Matter 

Megavoltage photons are packets of energy exceeding 1 million 
electron volts (MeV) each.* Although photons have dual wave and 
particle properties, photons (as their name implies) in radiation 
Dhysics, are treated as if they were particles. 

X-rays are photons produced by the interaction of charged particles 
with matter. In commercial x-ray equipment each photon is produced by 
the interaction of a relativistic electron with an atomic electron or 
nucleus of a target material. This process is called bremsstrahlung. 
The differential cross-section with respect to the photon energy, hv, 


for bremsstrahlung production, do >eouOmeaucharcedsnami cle anced 


brem 
d(hv) 
acting with an atomic nucleus is given by: 
Z 
doy em _ LZ- eno ie alae 
~d(hv) Se ee 


where e is the electric quantum of charge 
— Plank's constant h = (h/2n) 
Mm, is the charged particle mass 
c is the speed of light 
Z is the atomic number of the target nucleus 


T is the charged particle kinetic energy 


hy is photon energy 


ive) 
| 


= B(Z,hv/T) is a dimensionless quantity from quantum electro- 
dynanic theory and is of the order of unity (1) 
(Appendix 1 lists numerical values of physical constants). 


* 1 MeV = 1.602 x 102°; 
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The effect of interactions with atomic electrons can approximately 
by taken into account by replacing the factor ie TR ECI Yet ale Dy. 
Z(Z+1). The second order dependence of the cross-section on Z means 
that high atomic number materials are the most efficient targets for 
the production of bremsstrahlung. The inverse square dependence on 
charged particle mass explains why electron beams are always the 
accelerated charged particles used in commercial x-ray devices. 

The number of photons produced by bremsstrahlung per photon energy 
intervals at a given electron energy is proportional to Eqn. ee deand 
sO is proportional to l/hv. Therefore, the value of the photon number 
spectrum for bremsstrahiung produced by monoenergetic electrons is 
greater for lower energies. In the practical situation the electron 
beam interacting with the target is not monoenergetic. The initial 
energy of electrons has a small spread. More importantly, as the 
electrons interact with a target of finite thickness they lose energy 
as they traverse it. This results in fewer high energy photons and 
more lower energy photons than predicted using a "'thin-target"’ bremsstrah- 
lung approximation so the spectrum is strongly peaked at low photon 
energy. Consequently, the nominal megavoltage designation of photon 
sources produced by bremsstrahlung refers to the maximum photon energy of 
the beam. 

With minor exceptions the interaction of photons with matter occurs 
via the photoelectric effect, Compton effect and pair production. 

In the photoelectric process a photon is completely absorbed by 
an atom which becomes ionized by ejecting an electron. The photon 
energy is almost completely shared between the kinitic energy of the 


electron and the binding energy of the electron. Very little kinetic 
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energy is given to the ionized atom but most of the photon momentum is 
transferred to the atom. The probability of interaction is highest when 
the photon energy is near to but greater than an electron binding energy. 
For this reason, the photoelectric absorption probability for an element 
is discontinuous at the atomic shell binding energies. The interaction 


of the K-shell electrons contribute about 80% to the atomic cross- 


section. To a first approximation, the dependence of the atomic cross- 


SeCLi0N, ql» on Z and incident photon energy has the form (1): 


ge 


Gunes 
(hv)° 


ae 


where n is found empirically to vary between 4.0 and 4.6 depending on 
the photon energy. The photoelectric effect is the dominant mode of 
interation at low x-ray energies in high atomic number materials. 

At intermediate photon energies the Compton effect dominates*. In 
the Compton effect an incoming photon interacts with an electron. The 
photon is absorbed and another photon with a lower energy is created. 

The electron recoils with an energy representing almost all of the 
difference between the incoming and outgoing photon. The outgoing photon, 
usually termed the scattered photon, and recoil electron directions are 
in the same plane as the incoming photon direction. The directions of 


the scattered photon, $¢, and recoil electron, 6, shown in Figure l, are 


given by (2): 


* The Compton effect is the most important process from 30 KeV to 
20 MeV in water and from 600 KeV to 4 MeV in lead (9). 
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Figure 1. The Compton effect is an interaction of a photon with energy, 
hv, with an electron at rest. A scattered photon with energy, 
hv*, and a recoil electron with kinetic energy, T, result 
from the interaction. 
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COS *¢ = 1 ~ & (hv-T) Gland) 
tan 9 = (1.2.4) 
(= tan (9/2) es 
where <_ = ae 
m_c 
O 


Equations 1.2.3 and 1.2.4 were derived from considerations of conserva- 
tion of energy and momentum. Equations 1.2.3 and 1.2.4 suggest that the 
scattered photon and recoil electron directions are forward peaked at 
high incident photon energy. Scattered photons, each possessing a 
direction different from the incident primary photon direction, degrade 
the collimation of the primary photon beam. 

The differential cross-section with respect to recoil electron energy 


for the Compton interaction, do/dT, is given by (3): 


s : rr“m,¢°C(hy,T) (1625) 
m AT : : 
C(hv,T) = —>~5} | 2] YT BEE Terme’)? - mye)” | 
(hv-T) (hv) hv (hv) 

(pees) 

ef 

where r_ = 
fe) 2 
mc 


Klein and Nishina in deriving the above equations assumed that the 
electron is originally at rest and the binding energy of the electron 
to the atom is negligible compared to the energies of the primary 
and scattered photon and recoil electron. Both of these assumptions 
are equivalent to assuming a free electron at rest in the laboratory 


frame. These assumptions result in the total Compton atomic cross- 
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section, being proportional to the number of electrons in the 


at ’ 
atom, 1.6. ,° proportional to Z'. 

The interaction of a photon with the Coulomb field of a charged 
particle sometimes causes the photon to disappear with the creation of an 
electron and positron. Consequently, the energy threshold for this 
interaction is equal to the rest mass of the pair of 1.022 MeV. The 
remaining photon energy released is transformed into kinetic energy of 
the charged particles. When the Coulomb field is due to a nucleus, the 
process is called pair production. Triplet production is when an atomic 
electron, rather than a nucleus, is involved in the interaction. In 
this case, the atomic electron recoils with a considerable fraction of 


the available kinetic energy; consequently, the threshold for triplet 


production is higher than that of pair production. 


The differential cross-section for the pair production interaction 
with respect to charged particle (electron or positron) energy, 


dn/dT, is given by (4): 


hp PChy.T) (ine?) 


P(hv,T) is the kinetic energy dependent part of the differential cross- 
Section which can be expressed in a compact form if an energy parameter, 


v, 1s defined (4): 


ake (e283) 


Therefore, v is the ratio of total charged particle energy compared to 


the photon energy. 
P(hv,v) = |v" + (1-v)? + a) ee : : -c(y) 
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where c(y) is a dimensionless quantity usually much less than 1 which 
decreases rapidly with increasing y. y is the dimensionless screening 


parameter in the Fermi-Thomas model of the atom and is given by (4): 


AT Meena mndanlecsh (12870) 


In low Z materials and for megavoltage energies less than 20 MeV, 
the parameter, c(y), may be assumed to be zero. Triplet production may 
be included by replacing Z in Equation 1.2.7 with Z+1. 

The root mean angle between the direction of a secondary charged 
particle and the primary photon is given by (4): 


Sey! Qn os 
(<2) Seis ie (ea) 


O 


where q(T,hv,Z) is of the order unity and is given in Appendix 1. 

The total atomic pair production cross-section, ql 1s proportional 
to Ze so high atomic number materials have a large cross-section. 

When a positron and electron get close enough they annihilate 
liberating their energy. This usually occurs when the positron and 
electron have little kinetic energy. The rest energy is converted to 
two 511 KeV photons which are emitted in opposite directions. 

Other less important interactions of photons with matter occur such 
as photcdisintegration, Cerenkov radiation, Auger electron production 
and fluorescence. Of these, fluorescence is the most important, 
especially in crystal structures. Atoms can be excited due to the pas- 
sage of x-rays. When the atom returns to its ground state fluorescence 
radiation is emitted. When the process is delayed long after the initial 
excitation, the radiation emitted is called phosphorescence. 

The total cross-section for a beam of photons interacting in matter 


depends on the beam geometry. In narrow beam geometry both absorption 
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and scattering of photons are effective in reducing the intensity at the 
detector. Figure 2 illustrates a narrow beam of photons interacting 
in a thin layer of matter, dx, between the incident beam and a small 


detector. The number of photons that arrive at the detector is given 


by: 
N(x) = N(x=0)e Ut” Cal 
Me is the total linear attenuation coefficient and is given by: 
‘9 Ne OL ot 1) (ze) 
ui aed aa 


where No = Avagadaro's number 


A molar weight 


physical density 


D 
i) 


ae ao and q™ are the photoelectric, Compton and pair production 
atomic cross-sections, respectively. 

Broad-beam geometry assumes that the cross-section of the primary 
photon beam is wide. Only complete absorption of photons contributes 
to the attenuation. Scattering of photons away from the detector are 
compensated by scattering of photons toward the detector. The number 


of photons arriving at the detector from a broad beam is given by: 


Ne 
Ha A a 


ie 
ee a) (a7 1 Si) 
where oe is the Compton atomic absorption coefficient. 

g%q Tepresents the average fraction of the primary energy absorbed 


by the electron multiplied by the total Compton attenuation coefficient. 


oa" {hve tv) at (eeata) 


where hv~* is the energy of the scattered photon. The Compton 


atomic scattering coefficient, ee is a complementary quantity to Sone 


t om 

et epee sreptinr: al ue Totes ad erhyasies wre ance) 
eipfaetsiie 2he30ny St fier ae SqURT ails & = one F 
Stade Ene ved srebined ade anpd ease pee pets oe Se 


fevia Ar ot Ads ve eine sed) queen: 40 TA eee 


T toes Te =si 3 748) 6 meiisé shehct 


;a/ 


“exerts 2) cnbep autem , 
trig teat to ee 
Eel le 
TG; eae ets, tes ais a pe eee 
ag Lasts ah ne 
f-33-220%5 = isa) aamiee! CHaETOR 
iF3 nO erated eS nabtaaeZee Fes ok Scat sai» 
Tt Suez re hie Ad athe 
shine? oft, “cragsedsl sds Woe) Seger we gain ha 
d'tere ei 46sd lsgta a fe 4)oWe? tee) ond; iF sarin ¢ 
2 ds 9 a AM hen eae 


tek US gh} or Ven) oo ee ii . 


4 Le * 
$9960 AoLsrt Aa tey aie BEEK va eae 
7 
hide gers vert ses ualiq sce IO rotiaeee Saezers Sat 


‘Teka Pieo toeskonstia eI) Fasee SaF Hi patlqis itm moves : 


et 


Figure 


ro) 


Le 


Narrow 
Beam 
2 Sadao ecard 
f 
Thin Small 
Slabs Detectors 
Broad | 
Beam 
Z cel 
sceaaerniindss 
Pe Sete 
creas eae 
| ee 
a 
ERE 


Narrow beam geometry assumes that all scattered photons are 
not detected. Broad beam geometry assumes that as many 
photons are scattered toward as away from the detector. 
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At low energies a?s @Pproaches Ap and at high energies glq COn- 
tributes most of the attenuation. A listing of total attenuation 
eCOeificients, u,,and total Compton attenuation coefficients, ae , 
for water at selected energies appears in Appendix 2. 

The linear coefficients, Ue and u, are often normalized to the 
physical density, p, and then the term "ass'’ replaces "linear". The 
mass attenuation coefficient, u,/Ps for ice and water, are the same 
even though their linear attenuation coefficients are different. 

Other linear absorption coefficients are often used. The mass 
energy transfer coefficient, uy/ es and the mass energy absorption 
coefficient, Le Os, are defined analogously to the mass absorption 
coefficient, u/p- w/e takes into account the escape of Compton 
scattered photons. uy /'p takes into account the escape of Compton 
scattered photons, fluorescence and annihilation photons. Hen /? takes 
into consideration all the escape photons that uy/p does and includes 


bremsstrahlung photons. The general formulae for all these GOCTELCIENTS 


ner) s8 


Lee) ye) fo} T 
yesOr Tu 1k Gt hacia arty) (1.2.18) 


Where x can be ''k'', "a!" Or en 
Me are conversion factors which weight the average fraction of 
photon energy that is converted to charged particle kinetic 
energy taking into account the amount and form of energy that 
escapes the region of primary interaction. 


: Ab Oo 
For example, we have seen earlier that fo = f, =21 enhire= “tae 


(see Equation 1.2.16) 
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Sometimes it is convenient to compare the linear attenuation 
or absorption properties of two types of matter. The equivalent thick- 


ness of one substance compared to another is defined to be: 


ae 
(t_)° =— (2719) 
1s aeteait 


The primed and unprimed coefficients refer to two different materials. 

Since the linear attenuation coefficient depends on energy, the 
photon spectrum is going to be modified as it passes through matter. 

This especially affects low energy (below megavoltage energies) spectra 
because the photoelectric effect is strongly energy dependent. Low 

. energy photons will be depleted faster than high energy photons so the 
mean energy of the beam will increase. This is known as hardening of the 
beam. 

Hardening is not as important in megavoltage beams because few of 
the interactions (in low Z materials) occur via the photoelectric effect. 
Most of the interactions occur via the Compton Seeceuee Forward 
scattered photons tend to replenish the low energy part of the spectrum. 
In addition, the linear attenuation coefficient, over the energy range of 
megavoltage spectra, does not decrease rapidly with energy (in low Z 
materials). The photon spectrum of a beam of x-rays is often called the 


"quality'"’ of the beam. 


* If very high energies are employed, pair production could become 
important. The pair production cross-section increases with increasing 
energy so the spectrum would be depleted in high energy photons. 
“Softening'' could be said to occur. 
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1.3 Interaction of Electrons with Matter 


For the purpose of medical physics, a free electron can be 
thought of as a point particle with a rest energy of 511 KeV 
(rest mass = 9.11 x 1004 kg) and an electric charge of 1.602 x TO 
Coulombs. ” Free electrons are produced through ionization of atoms. 
The source of electrons in a linear accelerator is from thermionic 
emission from a heated filament or indirectly heated cathode. The 
photoelectric, Compton and pair production interactions also produce 
free electrons with a considerable kinetic energy. 

When an electron traverses through matter the Coulomb interaction 
between it and atomic electrons produce ionization of the atomic 
electrons and excitation of the ionized atom. Thisgresults in the 
kinetic energy of the electron, t, being deposited along its path which 
results in the electron slowing down. The energy lost in any given 
collision may be as low as a few electron volts or as high as one 
haif of the electron energy. The energy lost per unit mass of matter 
is defined to be the radiation dose absorbed by the matter. The 


collision or ionization stopping power is defined to be the average 


amount of energy lost per unit path length traversed, [at/dx]. 0, , and 
is given by: 
4n_z_ { (T+mc)? Ponies Vee 
| _ 27e 90°) O on o— 5 Li 
7o0el oda 1. 23 aes | 
ion =m, T(T+2m.c) moc I 
(Legs) 


* Since electron beams encountered in medical physics are not 
polarized, spin effects of free electron can be neglected. 


arn ad v ; wuss J AL te 


; a 
= 4 * - 
tf fei vet Pie ks | SO 
se : > ee os : 71 
rg : 7 ; 
ad has Fas Get? 6. Sees Seen aos 


fs \Z fniyahe (25° & i 49 ein eyed re 

. i y 7 
6 hor an By eens | iytcele se baa ose a Gi'x wees 
het “Fo noites ho. fea beet eT one xg 208 ot 
JPA untads. eth’ el) YeIoiesanoee se A= | Bn cd 
| ee ceo ORT Sas seed 
sAumoatys tig re hp a 


Pa ami aut # ta 


2 MSIFAN! AguITas eee 


ie ee cae 


mVveis Vue al L Sgr sritecnie wel Branch Serer) ; aig 
gud 2€ dutitees ca erty nips ol geeeet ie Bie Yt i } 
Ta) lia) 10 225i i 1s) =p Yarses sa! “cei haw 
Sir Vveiran, eds AES: Uy ecg eel Fey tA ey 38" 


diviel  -beasavend rane tek crea Saueo ee 
fit] "<bgeroveny Maps etek stead Saab bala 


a pay. hia ie a ’ - F 
i * +. a aiid ce he Pa - 
‘ ‘ ' in g 7. +. ii anf os 1H 
a L. j 0 


el) Paabst) acy, 


Ee Pos eSreiig Sisk ise mn a Gs 


where I is the average ionization potential of the atoms through which 
the electron is traversing. 
The factor (NZ /A) has the dimensions of cn? and represents the number 
of electrons per cm of a material and is called the electron density. 
When T is large (T>>m,c*) the collisionai stopping power increases 
Slowly because the factor in front of the logarithm approaches unity 
and taking the logarithm modifies the approximate T° dependence of its 
argument. At low energies (T<<m,c°) the stopping power increases 
approximately inversely with decreasing electron kinetic energy. 
As was discussed in Section 1.2, bremsstrahlung is produced by 
fast electrons interacting in matter. The stopping power due to 
bremsstrahlung, often called the radiative stopping power, can be 


found from the differential bremsstrahlung cross-section (Equation 1.2.1): 


N_o T 
dT Rp ite) P Bs 
fa ae hs {i PG ae (1.3.2) 
brems " 


The collisional stopping power is much greater than the radiative 
stopping power except at high electron energy or high atomic number 
materials. The total stopping power is just the sum of the collisional 
and radiative stopping powers. Appendix 3 lists the total stopping 
powers of electrons in water as a function of electron energy. 


The average path length, P.L., traversed by the electron is given 


by: 
T. | 
i 
Pale =f [72] ot dT + rel (GPSS .5)) 
iE 
min 


where ive TSeciceinicialecnercy 
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ie is the kinetic energy when the electron can be considered 
to be aimost stopped 
Eb ecidual is the residual path length after the electron has 


reached an energy, {pea 
All electrons do not have the same path length because of the con- 
siderable fluctuation of collisional energy loss. 

The range of an. electron in matter can be defined to be the 
distance between the point of entry of the electron in matter and the 
position where it stops. The mean range is always less than the mean 
path length. The electron having a small mass is easily scattered 
by the Coulomb field of the atomic nucleus. These interactions are 
usually elastic so only the electron direction is changed. The 
multiple scattering of electrons after passing through a path, 2, is 
often assumed to have a Gaussian angular distribution. The mean square 
mass angular scattering power e2/o°k of a material is a measure of the 
angular spread of an electron beam. 2 /p.k is directly proportional 
to had so the scattering increases rapidly with high atomic number 
materials and decreasing electron kinetic energy. Appendix 4 lists 
¢/o+k as a function of kinetic energy for electrons in water and lead. 

If an electron source is located at a point and directed in a 
direction, z, then the electron fluence, ¢, at a distance, r, from the 
point compared to the initial electron fluence, ae after it has 
traversed a distance, z, intO-4 material is (5): 

P 30/28 
6 = _—__,— e374) 
l+tan 6 
where 6 is the angle between the z and r directions 


ibeaer = r?/2° 
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1.4 Medical Electron Linear Accelerators 

The x-ray tube common to all diagnostic x-ray units is the simplest 
type of electron accelerating device. Potentials larger than 1 MVp 
across a space of a few centimeters cannot be maintained because 
dielectric breakdown occurs. In order to produce high kinetic energies, 
a linear accelerator requires an accelerating electric field travelling 
in the same direction and velocity as the electrons. In a linear 
accelerator, an electron can gain several MeV/m, yet the potential 
difference between any two points at any time can be less than several 
hundred kVp. 

The basic accelerating unit is the microwave resonance cavity. 
A cavity is merely a microwave waveguide that has been enclosed at the 
electric field minima points. Figure 3 illustrates the relation 
between the magnetic and electric fields in a microwave cavity for 
one complete oscillation. The cavity oscillation can be set up by 
connecting the cavity to a source of microwave radiation such as a 
klystron or magnetron. 

A klystron comprises two coupled microwave cavities (6) - one 
called the buncher and the other called the catcher (see Figure 4). 
The oscillating electric field from a low amplitude, high stability 
microwave source alternately retards and accelerates a continuous 
Stream of low energy electrons produced by a hot filament. This 
modulates the velocities of the electrons leaving the buncher. The 
buncher and collector ees are joined by a drift tube which allows 
higher velocity electrons to catch up to lower velocity electrons. 
When the electrons enter the catcher cavity they generate strong 


retarding electric fields by inducing charges on the ends of the cavity. 
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Figure 3. One oscillation of an electromagnetic cavity. The electric 
field 1s a result of charges at either end of the cavity. 
The magnetic field is due to charges moving along the cavity 


in response to the electric field?” “Microwaves can induce or 
be produced by such oscillations. 


Iss, 


i oa Oe all | r- 
| ~ e a e é 4 
ee > > —~ 
* . . . 

oo) be iy Ne 6 oO ~ — 
. i i ~j 
Pall —@» (2 =~ a 
' é < uy 
| ee _——- i a 
i e 

— —— — 


| ime te oeematied _——ene 9 
’ — 

cn eg emt ne i 

— } 


apt ee 0 eee ee 


a rr 


\ 


\ 
> ot. 
a eee a 2 
7] 


silenpeV 
. bisit 
agticadd 


ie ent. J ives js dares ued 
wis a. eS 16 =) 


eden oit fF ie 


Low Level Amplified 


Microwaves Microwaves 


Filament 


Buncher Catcher 
Cavity Cavity 


Beam 
Collector 


B) 


Figure 4. A) Retarding and B) accelerating phases of the operation 
of a buncher. A weak microwave cavity modulates an electron 


beam. The kinetic energy of the beam is transformed into an 
amplified microwave signal. 
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This reduces the electron velocity. The modulated beam current arriving 
at the catcher induces high amplitude microwaves utilizing most of the 
kinetic energy of the electrons. The electron beam collector is a 
dump for removing the electrons. The kinetic energy not converted to 
microwaves heats the collector which must be water cooled. Some low 
energy bremsstrahlung is produced in the collector which must, therefore, 
be shielded. The microwave power is produced by a klystron in the 
Mevatron -20 accelerator. 

In the Mevatron -6 accelerator the microwave power is produced by 
a magnetron. A magnetron has a flat cylindrical geometry. A central 
cylindrical cathode is surrounded by an outer anode which is separated 
by a drift space. A D.C. magnet field is applied perpendicularly to 
mene plane of the cylinder. A pulsed-electric field is applied radially 
inward. Electrons spiralling from the cathode to anode induce an 
alternating charge distribution on the anode which produces, and has the 
same frequency as, the microwaves. Most of the electron kinetic energy 
is converted to microwave eres 

The two types of accelerator structures in use are the travelling 
wave and standing wave side-coupled designs. Both accelerators used in 
the study were standing wave side-coupled designs so the travelling 
wave design will be mentioned only briefly. Both types of accelerators 
accelerate electrons that have been produced and grouped in much the 
same fashion as a buncher cavity of a klystron. 

In a travelling wave accelerator, radiofrequency power enters 
the first cavity (see Figure 5). The electric field travels down 
the waveguide which is a series of resonant cavities at the same 


group velocity as the accelerating electrons. Since the electron 
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Schematic diagram of a travelling wave accelerator accelerating 
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velocity increases, the spacing between the cavities are increased in 
length toward the distal end of the waveguide. Upon reaching the 
distal end of the waveguide, the radiofrequency power is dissipated 
in a resistive load. 

In a standing wave accelerator the microwaves are not dissipated 
with a concomitant saving of energy, but are allowed to reflect such 
that a standing wave is set up throughout the waveguide (see Figure 
6). Since half of the cavities have a zero or small electric 
field at all times and so play a small role in acceleration, these 
Cavities may be moved off axis. A waveguide with such a configuration 
is said to be side-coupled. The separation of the roles of the coupling 
and accelerating cavities allows optimization of the size and shape of 
each cavity type. This increases the amount of energy that can be 
gained per unit length which reduces the waveguide length. Figure 
6 illustrates hypothetical longitudinal and cross-sections of a 
standing wave side-coupled accelerator waveguide. The phase velocity 
of the a cna tates are matched to the electron velocity such that 
a group of electrons experiences the same direction and magnitude 
of electric field throughout its passage along the accelerator structure. 

To be effective in treating patients from all directions, the 
accelerator must be able to be rotated about an axis. Since the 
accelerator structure is more than a meter long it has proven convenient 
toplace the accelerating waveguide parallel or nearly parallel to the 
axis of rotation. In order to direct the beam perpendicular to the 
axis, a magnet must be employed to bend the beam. To accomplish this the 
beam could be bent through a right angle, but a 270° achromatic bending 


magnet is usually employed (see Figure 7). 
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Figure 6. Schematic diagram of a side-coupled standing wave accelerator 
accelerating section. The forward wave travels from left to 
right and the reflected from right to left. The waves travel 
between accelerating cavities via coupling cavities. 
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A 270° achromatic beam-bending magnet. 


Figure 7. 
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The relationship between the kinetic energy, magnetic field strength 


and the radius of curvature of a particle in a magnetic field is given 


mc Tém 4 
ene 7 ee = Gers) 
LAS 


where B is the magnitude of the magnetic field. 


by: 


An achromatic bending magnet has a greater magnetic field near the 
outside of the magnet. Electrons with higher than average kinetic energy 
will travel to outer radial parts of the field where they will be bent 
more. Slower electrons experience an opposite effect. The net result 
is that all electrons, regardless of small deviations in their kinetic 
eilerey, arrive at the same point. This is desirable because if a target 
(see Figure 8) is placed at this point the bremsstrahlung produced will 
emerge from a limited area. Photons from a small source when collimated 
will produce sharply defined edges, i.e. a small penumbra (see Figure 9). 
Additionally, an achromatic magnet results in a more stable delivery of 
electrons to the target. 

The treatment head can be divided into the fixed head assembly and 
the movable collimator. The fixed head assembly consists of the target, 
electron absorbers, primary collimator, field flattening filter, monitor 
ion chamber and field light mirror (see Figures 8 and 9). 

The electron absorber is a low Z material placed in the beam after 
the target to stop low energy electrons. The primary collimator limits 
the beams angular width to 203. 

The bremsstrahlung produced by the target is forward peaked in 
intensity so there are more photons directed along the central axis than 


at points off axis. The field flattening filter is machined in a conical 
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Figure 8. Scale drawing of the treatment head of a Mevatron -20 medical 
hineangacceleracor. 
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Figure 9. Diagram of the field defining system. 
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shape to attenuate the photon beam preferentially along rays near the 
Sentralsaxis. The result 1s a photon fluence distribution which is 
relatively uniform over the entire cross-sectional area of the beam. 
Since introducing attenuation into a photon beam tends to harden the 
beam, the lateral distribution of the photon spectrum is not uniform. 
This can produce a nonuniform distribution of dose. 

A multiplate ion chamber is usually placed in the fixed head 
assembly to measure dose rate, integral dose and field symmetry. This 
allows Rena onne of the dose received by patients (timer systems are 
always used as back-ups) and is used as part of the verification of the 
safe operational status of the accelerator. 

A field light, the source of which appears to be at the same 
position as the x-ray source, illuminates the field. This illusion is 
obtained by placing a thin mirror in the field and placing the light 
source off axis (see Figure 9). 

The collimator system can consist of steel, lead, tungsten alloy 
or depleted uranium blocks thick enough to reduce the primary photon 
fluence to a fraction of a percent of the fluence in the field without 
producing photoactivation neutrons. A penumbra at the field edge occurs 
because of the finite width of the electron beam and target thickness 
(see Figure 9). The field boundary is defined to be where the dose falls 
@oe50%5 0f thetdose at the central axis. The amount’ of the collimator 
jaws covering the field are continuously adjustable to produce square or 
rectangular fields with dimensions up to 30 cm x 30 cm. 

A removable thin lucite tray with cross-hairs inscribed on it can 
be placed directly beneath the collimator jaws (see Figure 8). The 


shadow of the cross-hairs in the illuminated field define the position 
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of the central axis. 

The beam axis of rotation for most accelerators is 100 cm from the 
photon source. Accessories such as beam wedges or blocking filters are 
often placed between the collimators and the patient at a distance of 
60 cm from the source in the Mevatron accelerators on a Lucite accessory 
or ''shadow' tray. When there are no accessories or beam modifying 
devices in place the field is said to be open. 

A phantom consists of tissue equivalent material such as polystyrene 
or water, in the case of soft tissue, in which a radiation measuring 
device has been placed. The distance from the source to the phantom 
is specified by the source to (phantom) surface distance, SSD, and the 
source to probe distance, SPD. The phantom material used for all the 
experiments to be reported is polystyrene (electron density = 3. 50x10°> 
e /cm, physical density = 1.08 g/cm”) and the probe was usually a 


thin-window parallel-plate ion chamber. 
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1.5 Ionization Dosimetry 
The absorbed dose, or simply, dose, is the amount of kinetic energy 


deposited in a small mass of material or: 


wits 
Dose = i (Seely) 


Dose is not a differential quantity. The dimension of the energy 
absorbing mass, Am, must be much larger than the average distance 


between ionization-excitation events given by: 


Average Distance Between Ionization Events = sop (eZ) 


where W is the average energy lost in an ionization-excitation event. 


Typically, the average distance between events produced by fast electrons 


is less than a micron so the dose can be defined for very small masses. 
Dose is defined in terms of an amount of kinetic energy absorbed. The 
amount of rest energy deposited is not included so dose cannot be 

defined in terms of AE where E = Tm,c’. 

The Gray (1 Gray = 1J/kg) is the SI unit of dose, However, for 
historical reasons the rad = 1 centiGray is often unofficially used. 

An ion chamber measures the amount of charge accumulated due to 
ionization of air, One of the simplest ion chambers in design and the 
type used in most of the experiments conducted is the parallel-plate 
ionization chamber (see Figure 10). In essence, it is a single element 
parallel-plate capacitor separated by air. A potential is maintained 
between the plates. The potential is sufficient such that all ion 
pairs produced in the air are accelerated to the plates but not large 
enough for the accelerated ions to produce more ionization of the air. 


There is some recombination of ion pairs as they traverse the chamber, 
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Schematic drawing of a parallel-plate ionization chamber. 
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but if the ionization current density is low the amount of recombination 


will be negligible. The amount of charge required is determined using 


an electrometer by measuring the amount of charge placed on the parallel- 


plates required to maintain a constant potential. 

Although an ion chamber reading does not yield an absolute measure 
of dose, it is possible to theoretically relate the ionization in the 
ion chamber probe with the dose received in matter in the absence of 
the probe (7). In practice, the conversion relating the amount of 
ionization in a given probe to dose is established by a standards 


laboratory. The amount of dose is given by (8): 


Dose = CCA K, K, PR Aas 


where R is reading of charge accumulated per unit mass of air 
C is the calibration factor 


CA is the energy dependent conversion factor between charge/mass 
and dose 


K, is a temperature and pressure correction factor 


K» is a ddetor that eorrects tor difference Such asthe: quality 
of the radiation beam 


P is a factor that takes into account the perturbance the probe 
makes in the fluence of the medium 


Dose deposition due to photons is a two-stage process. First, 
Charge particles are set in motion, then the charged particles slow 
down depositing their kinetic energy as dose. The concept of KERMA, 
kinetic energy released in the medium, describes the first step. 
KERMA is given by: 

Vk 


KERMA = hy — 6 
0 ry 


(154) 
where iy /'p is the mass energy transfer coefficient 


%, is the photon fluence in particles per unit area 
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Since the photon fluence is decreasing as a function of depth 
in the phantom, the KERMA decreases as well. The charge particles 
set in motion are mostly forward directed so their fluence as a function 
of depth initially builds up. Once set in motion these charged 
particles slow down and are stopped so their fluence reaches a 
maximum intensity at some depth. At this position, maximim fluence 
will produce maximum ionization in a probe placed there. This point 
is called oan After cue the charged particle fluence decrease 
follows the KERMA decrease. The following differential equation has 


been used to describe the region where the dose builds up:* 


Se bl Pe (1536) 
ai ar ee A ¥ t (X) 


and the solution is 


rs = (Gy,/ OX (ise) 
¥ As) = (ros y/P 
Wy Oo ae View) ODE caltek/iol x 
¥,tQ9 = eee (: y -e re" (10528) 


where ¥ (x) is the total amount of energy fluence due to a full spectrum 
Y of photons and is given by: 


¥ (x) =[¢ (hv,x)dhv (ie) 
ay, sy; 
spectrum 


where COs is the initial photon energy fluence 


» +(x) is the energy fluence of charged particles and is defined 
: analogously to Equation 1.5.7 


* Equations 1.5.6, 1.5.7 and 1.5.8 have been modified from similar 
ones by Cassen, Corrigan and Hayden (10). 
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xis the thickness ot phantom (an units of mass/area) 


u is the energy weighted mean linear energy absorption 
coefficient and is given by: 


» ak aca 


spectrum 
©, (hy ,x)d (hv) 


(ese) 


spéctrum 


u_+ is the mean "energy absorption coefficient" of charged 
particles set in motion at a point (in units of mass/area) 


The above expression is analogous to the decay of a parent-daughter 
isotope if the daughter isotope initiaily has zero concentration. The 
effective energy transfer coefficient and the mean energy absorption 
of charged particles is depth dependent due to hardening of the beam, 
and the production of scattered photons. However, to a good approxima- 
tion they can be considered to be a constant in the build-up region. 
The energy absorption coefficient can be used instead of the energy 
transfer coefficient at low energies in low Z materials. It has 
been assumed that once produced the energy fluence of electrons set 
in motion at a given point in the phantom decays exponentially with 
depth. The first term on the right hand side of Equation 1.5.6 
represents the total KERMA of the beam. 

The depth of the maximum dose, dues occurs at the point of the 


maximum charged particle energy fluence and is given by: 


ee 
ae —— te (055 ahi) 
maxX =, tp ny 

The output factor is defined to be the reading at ch for a 


given field size and SPD divided by the reading at a field size of 
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10 cm x 10 cm at the same field size. The output factor has a small 
dependence on SPD. 

The dose at a given depth is often normalized to the dose which 
occurs at a reference depth. There are two common procedures for doing 
this. They are the percentage depth dose, P(x,X,»W,,SSD, (hv) 4) 
and the tissue maximum ratio, TMR (x,W, (hv) J: 

The percentage depth dose depends on the depth (measured from the 
surface) in the phantom, x, the depth to the reference point, Xo» the 
field dimension measured at the surface, Wo» the source-to-surface 
distance, SSD, and the nominal photon beam energy (Iv), 4 Nets 


ax 
percentage depth dose is defined to be (8): 


_f Dose at x 5 
P(x,xX,,W.,SSD, nee JGeee om GIS 5.12) 


The position of the dose measured at points x and X, are shown in 
Figure 11. The percentage dose is obtained by keeping the SSD 
constant and determining the dose at each depth by moving the measuring 
probe. 

The tissue maximum ratio is a ratio between the dose as measured at 
a depth, x, in a phantom compared to the dose as measured at deer. The 
measurements are made by keeping the source-to-probe distance, SPD, 
constant and the depth from the surface is varied by changing the 
thickness of overlying material. The field dimension is measured 
at the source-to-probe distance. The tissue maximum ratio is defined 
CODE (8) 
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Figure 11. 
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The arrangements for measuring percent depth dose is shown 
on theleft and tissue maximum ratios on the right. The 


position of maximum dose in the TMR measurement, d ax? 1s 
not shown. 
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The dose at point x is shown in Figure 11. The percentage depth dose can 


be obtained from tissue maximum ratios by the expression (8): 


TMR(x,W, (hv). ) / SSDx, \7 
P(x,X,5W,SSD, (hv) J " TRG, 1 OH) = Sspex x 100% 
(1.5.14) 
2 


SSD+x 
takes into account a 1/r? fall-off in the primary fluence 


SSD+x 
between points Xo and x. 


where 


Using Equation 1.5.3, the tissue maximum ratio can be found: 


Cy (OK, ODROD 


Fran Gay) Epa) 


Cia Seeks.) 
0 max 


TMR = 


In the build-up region the quality of radiation does not change appreciably 
so the ratio k,Q)/k, (do) can be taken as unity. The ratio C, (x) /C, (do) 
is approximately unity as well.* Therefore, the tissue maximum ratio 

is just the ratio of readings taken at a depth, x, and at se 


Using Equation 1.5.8, the tissue maximum ratio is predicted to be: 


eae oo et/P)X 
: “TT COL (renee (sero) 


Fquations lwo. Scand 12 5.bo predict: that: the dose at) the proximal 
surface of the phantom (i.e. x=0) is zero. The reason is that no charged 
particles have been set in motion. The low surface dose can have important 
clinical effects. For example, the skin is a fairly radiosensitive tissue 


and often is not affected by deeper seated cancerous tissue. When this 


see =water 
C, (x)/C, (d_4x) iS approximately equal to S) Cia, Fale (7) where 


d over the 
b ccer x) is the stopping power ratio of dete to air average 
Sean) Oe ee at a depth, x, in the phantom. Nahun (53) has shown that 


qwater (.._ 0) is within 2% of 


cS Mv. a 


eee ax) in the energy range between 6 MV 
m 


38 


? ~ 
u - 7 — : ¥. 
<2 tee Wigel’ sgstrs sad SF: Ape orig Ri -mpovs at 
(3) aciccerers Qi Ve cots peng) et thee? 


pe) Fee Z ; if 
a = Stone pane: 1 eae ft) 


as ssi Fe. 
ff 
ef at 4 “iba ohm) i I> Te: ] Aa fe raat ae) 
x8 4 
id 
‘a? 
as 
bien? oti fine « . are LEM SUeeRS pis 
4 
§ = ~~ SP 
a _— Mm —s _ : — a 
7 | ¥ + a : = nd 
a ot | 
SHO SUS 310 i, MOLSELDSS 1 Li 
CG } all af <a - i e [ 


pes Migaut os (+g iced Sant es, 
y i - “a nee , 
~ ¥ ' ‘ 
ef OF Osi OP 2th J) a(r=7rT9 afl? iB 
¢ é 
Mine, n, gh, 
e a? z 
} _ —_—~ =" —_ mila . « - 
ime » ~ a a “— oa rhe 
ae 7 {> at a i F. ¥ 
oo a) =~@ i-" > 


‘bid sales vageais euoreone) beasgoe t8taso 


ise il iarats bz’ ae 
en STAT oy, a 
‘& mat. = 


is the case, the dose to the skin should be as low as possible. 

In practice, the surface dose is non-zero. In part this is due to 
backscattered radiation. The backscatter dose can be divided into two 
components. Some dose is due to forward directed primary photons which 
produce charged particles set in motion near the surface of the phantom 
which scatter backwards. The rest is due to photons scattering backwards 
and producing charged particles set in motion which are then mainly back 
directed. 

The rest of the surface dose is due to contamination. The 
contamination component of the beam is defined as charged particles 
and scattered photons produced by interactions of the primary beam 
with material outside the phantom. Since the dose is directly due to 
the slowing down of charged particles, the surface contamination dose 
is due to charged particles produced outside the phantom. 

The primary photon beam is defined to be photons emerging through 


the collimators which appear to have been produced by bremsstrahlung in 


x 
the target. 


* Other photons which are indistinguishable from bremsstrahlung produced 
photons are also considered primary. These include bremsstrahlung 
photons that have been forward scattered by beam modifying components 
such as the field flattening filter or annihilation photons directed 
parallel to the primary photon beam. 
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1.6 Review of Previous Work on Contamination of Megavoltage Photon 


Beams 

The University of Illinois (Urbana, illinois) betatron was the 
first device used in radiotherapy to produce a megavoltage photon beam. 
This betatron could produce bremsstrahlung photons from electrons 
with kinetic energies of 5 MeV to 20 MeV. The presence of electron 
contamination was evident in the beam. Quoting from a 1942 paper by 
Koch, Kerst and Morrison (11): 

"Experience in taking these [depth-dose] data has 
shown that the surface dose can be greatly affected 
by stray electrons striking the phantom. These 
electrons arise from two sources; some are original 
beam electrons which escape in large numbers from 
the acceleration chamber and others are Compton 
electrons scattered from objects near the x-ray 
beam. To obtain consistent results the primary 
beam was stopped by absorbers placed close to the 
target and the secondaries [contamination electrons] 
were avoided by keeping scattering objects from 

the vicinity of the beam. The magnetic field of 
the betatron undoubtedly removes a great number 

of secondary electrons." 

The surface dose, when compared to the maximum dose produced by the 
Illinois betatron, decreased as a function of increasing nominal 
beam energy. 

A 1941 paper by Johns, Darby, Haslam, Katz and Harrington (13) 
described the depth-dose and isodose distributions from a 22 MeV betatron. 
There is no mention of the production of contamination. They claim depth- 
dose curves are independent of field size or source to surface distance 
(retrospectively, there was some small dependence on these factors). 
The measuring equipment employed did not allow measurements to be made 
less than 0.5 cmMfrom the surface and the largest field was a 10 cm 


diameter circle measured at 70 cm SPD. These conditions are not 


conducive to observing contamination. 
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The absence of a build-up dose dependence on field size found on 
the Saskatchewan betatron was contradicted by Laughlin, Beattie, Lindsay 
and Harvey who found a dependence on field size when the Illinois 
betatron was operating at 25 MV (12). 

The introduction of the oe” isotope therapy devices clearly 
indicated that megavoltage primary photon beams were contaminated. 
Johns, Epp, Cormack and Fedoruk (14) determined that the beam was 
being contaminated mainly by electrons. They found that the surface 
dose, as measured with a thin window parallel-plate chamber, increased 
as a function of field width and decreased with greater source to 
surface distance. The contamination resulted in a shift in d ax towards 
the surface and produced an elevated value of the maximm dose. They 
recommended a distance of 20 cm between the diaphragm (secondary 
collimator) and the patient's surface. An aluminum absorber placed in 
the beam following the primary collimator reduced the amount of 


contamination. 


The results of the measurements of Johns et al have been confirmed 
fom Gove beams by many other investigators. The dose in the build-up 
region increases as a function of increasing field size (15-21). The 
dose at depths below dan increases with decreasing distance to the 


a beams produces a dramatic 


Sourcec(iS): Di2e20" 22) in Contamination “ofece 
shift in the position of the maximum dose (15, 20, 21). 

Following the theoretical calculations of Wilson and Perry (23) in 
1951 which showed that medium Z materials could be used to reduce 
electron contamination, many authors experimentally confirmed their 


findings. The best electron filter materials were claimed to be 
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authors used low or high Z materials as electron filters, for example, 
lead, Z=82, (17); leaded glass (25) and perspex (21). 

In contrast to the findings that materials placed in the beam re- 
duced the amount of contamination, it was found by authors that accessories 


ae field increased the amount of contamination (26). 


placed ina Co 
In linear accelerators a great deal of contamination is produced 

by the introduction of accessories into the beam. Rao, Pillai and 

Gregg (28) measured contamination at many beam energies produced 

by accessories. They found that the amount of contamination 


increased with accessory tray thickness up to a certain 


thickness (about 8 mm for a lucite accessory 29 cm above the phantom 
when the field size was 25 cm x 25 cm in a 6 MV beam) and was constant 
with increasing thickness for greater thicknesses. The amount of 
accessory tray contamination increases with larger fields and smaller 
source to surface distances. Other authors have confirmed these 
results for linear accelerators (19, 25). 

Rao, Pillai and Gregg (28) found that a lead sheet placed under 
the accessory tray reduced the dose at all depths below or i Gray 


et contaminant filtering, suggested 


(25), based on the success of Co 
using a leaded glass accessory tray in linear accelerators and Wu (29) 
recommended placing a tin sheet beneath the accessory tray in 10 MV 
photon beams. 

Scrimger and Kolitsi (30) studied the effect outside the beam 
when scattering layers were placed inside the field. An 8 MV linear 
accelerator was used. The relative dose outside the field fell off 


rapidly with depth suggesting that the scattered radiation is mostly 


electrons. The amount of scattered radiation outside the field increased 
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with field area and decreased with increasing distance to the scatterer. 

A field flattening filter is designed so that the lateral 
distribution at ora for megavoltage photon beams is as flat as possible. 
However, it has been found that the surface dose is generally higher 
near the central axis compared to the dose near and beyond the beam 
poses e(19, 21, 25,50, 250). 

As was observed for (0°? beams, megavoltage photon beams from 
linear accelerators have an increased dose due to contamination when 
the field size is large or the source to probe distance is small 
Wee eo 5255 O45 to 5). 

Contaminated accelerator beams produce a shift in dag through 
the entire range of energies currently employed in radiotherapy; for 
E2emple dt cv (Zoey poem (18, 25) 2 LOM (58)5° 25 MV (18,935 /,0 359)5 
34 MV (34) and 45 MV (40). 

Penetration characteristics of the contamination component have 
been measured in an number of different ways. Dawson (34) produced 
what he called ''difference curves". The difference in dose between 
neighboring field sizes at a fixed depth and source to surface 
represents the effect of an increase in field size. At shallow 
depths the difference curves decrease rapidly. Attenuation of the 
dose due to contamination electrons is cited as the cause. These 
measurements were done in a phantom greater in extent than the field 
so some of the difference between the dose at neighboring field sizes 
is due to photons scattered in the phantom arriving at the measuring 
point. This is evident in Dawson's data because the difference curves 
increase with depth when measurements are made at a depth beyond d ax’ 


Marbach and Almond (37) placed a lead pyramid-shaped block on the 
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accessory tray to shield a detector from the primary beam. Contamina- 
tion was then claimed to be the only component of the beam arriving at 
the phantom. They measured depth-dose curves for the contamination 
component. Unfortunately, the shallowest depth-dose measurement obtained 
was 2 cm from the surface. A key to the validity of such an approach is 
to confirm that the lead block does not shield or produce a significant 
amount of contamination. Surface dose measurements taken before and 
after the placement of the lead blocks would have determined if there 
was any change in contamination. 

Marbach and Almond subtracted the contamination curves from the 
total curves to yield the primary depth-dose curves. The primary 
depth-dose curves did not have a shift in aa and when normalized to 
their own maxima were independent of field size. 

Marbach and Almond went on to hypothesize that electrons were 
not the cause of the cee shift. They placed a 1/4 inch aluminum 
plate in the field below the collimators in the anticipation that this 
would remove all of the electron contamination. Since there was little 
change in the depth dose curve, they conciuded there were no 
contamination electrons in the beam. Their interpretation was not 
necessarily valid. As many electrons may have been produced as were 
stopped and scattered. Indeed, even though their contamination depth- 
dose did not include shallow depths, a rapid increase in dose with 
decreasing depth is evidence for an electron contamination component. 

Biggs and Ling (41) repeated the method of primary beam blocking 
developed by Marbach and Almond (37). They produced contamination 
attenuation curves that included the dose at the surface and at other 


depths below oe These curves exhibit a rapid fall off at shallow 
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depths and tend to be an asymptotic constant at deeper depths. They 
did not subtract the contaminant component to obtain the primary 
build-up curves. 

Clifton Ling and colleagues (41, 43, 44, 45) have employed electro- 
magnets with maximm central magnetic fields between 1.8 kG and 
3.9 kG to sweep electrons from the field. Padikal and Deye (42) 
used a permanent magnet with variable gap pole pieces. The smallest 
central field strength still capable of sweeping all of the 10 MV 
contamination electrons produced between the target and the magnet 
was 0.5 kG. 

Biggs and Ling (41) showed that the build-up curve foray 25° MV 
photon beam was field size independent when the magnet was sweeping 
contamination electrons from the field. 

Padikal and Deye (42) defined a parameter, «, to be the difference 
between the reading of the ionization chamber without the magnet and 
the ionization reading with the magnet normalized to the reading at 
d ax’ Since the only difference between the readings was an absence of 
electrons produced between the target and magnet, « represents the 
component of the TMR curve attributable to those contaminant electrons. 
The electron contaminant TMR curve had a peak between 1 and 2 mm and 
rapidly decreased to zero at about 2.5 cm. Despite evidence to the 
contrary, they claimed that the shift in ee with field size is due to 
a relative enhancement in the soft x-ray component as the collimator jaws 
are opened. 

Ling (43) has shown that a magnet in a 4 MV photon beam produced 
noeshiftsin diye There wseauslicht cield size dependence of, the 
surface dose and the build-up curves at shallow depths. This is 


interpreted as due to electrons originating from the air volume between 
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the magnet and the phantom. Ling, Schell and Rustgi (45) repeated the 
measurements at 10 MV and found similar results. They measured the 
build-up curves with the magnet in place at two different source to 
surface distances to see if the contamination that was previously 
attributed to air was SSD dependent as theory would predict (37). 
Their results were inconclusive so they have stated that they will 
repeat the measurements with a helium-filled bag to replace the air 
volume. 

Ling, Rustgi and Gromadzki (44) have measured the production of 
secondary radiation by 10 MV photons from scatterers placed in the beam.’ 
The measurements were done outside the primary beam field when a magnet 
was off and on. The depth of ionization curves for scattered photons 
and electrons were obtained separately. They investigated the amount 
of electron production as a function of the atomic number of the 
Scaccerer.. Fora target thickness. of 1.7 gm/em? and a scattering 
angle of 10; high Z materials such as lead produced less ionization. 

Recently, Nilsson and Brahme (46) have done calculations predicting 
the nature and amount of contamination due to scattered photons using 
the Monte Carlo method. The scattered photon spectrum for both 6 MV 
and 21 MV does not change appreciably between the central axis and the 
field edge for a 20 cm circular field. The maximum energy of the 
scattered photons is about 85% of the maximum energy of the primary 
spectrum. The peak in the scattered and primary spectrum occur at the 
same energy although the relative number of photons in the peak are 
greater in the scattered photon spectrum. The absorbed dose at oa 
due to scattered photon contamination for both 6 MV and 21 MV photons 
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field. The main source of scattered photons was the collimator for a 
6 MV photon beam and the field flattening filter for a 21 MV photon 
beam. 

The main disagreement in the literature concerns the type and 
origin of contamination. Some authors say electrons and others say 
scattered photons are the main cause. Every component between the 
source of primary radiation and the phantom has been suggested as the 
main source of contamination. Table 1 summarizes the conclusions 


of a number of authors on the nature and source of the contamination. 
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2. EXPERIMENTAL TECHNIQUES AND RESULTS 


The documentation 

of experimentation 

often has too much elaboration 
and not enough elocution. 
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je Open Field Build-Up Curves 


Tissue maximum ratio curves were measured for 15 MV and 6 MV photon 
beams from Siemens Mevatron -20 and Mevatron -6 linear accelerators. A 
schematic diagram of the beam defining head of the Mevatron -20 is 
shown in Figure 8. The Mevatron -6 has a similar design. All 
determinations were made at constant source-to-probe distances (SPD) 
using 0.16 cm (1/16'') thick and 0.64 cm (1/4'') thick square polystyrene 
siabs.. the phantom slabs used in the determination of the open field 
build-up curves each had cross-sectional dimensions of 25 cm x 25 cm. 

The detector used in this section and others was a Capintec 192A 
electrometer with a PS-033 thin window parallel-plate ionization 
chamber which has an effective volume of 0.5 mt. The entrance 


window was aluminized polyester film 0.5 mg/cm” thre kk. 


The field size is referred to as the lateral field dimensions 
defined at 100 cm. Only square fields were used. 

The central axis build-up curves at SPD = 100 cm for a 15 MV beam 
for various field sizes are illustrated in Figure 12. These build-up 
curves demonstrate a field size dependence. There is little backscatter 
expected at 15 MV so the dose at the surface should be very small and 
only weakly dependent on field size. Therefore, most of the field 
Size dependence and the elevated surface dose can be attributed to con- 
tamination. 

The open field build-up curve at SPD = 100 cm for a 6 MV beam are 
shown in Figure 13. The build-up curves were measured at the central 
axis for various field sizes. Greater field sizes result in larger 
tissue maximum ratios (TMR) for the same depth. Since a very similar 


effect was observed at 15 MV, the open field build-up curves at the two 
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Figure 12. 


Central axis open field build-up curve at 100 cm SPD for a 
15 MV beam. 
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Figure 13. Central axis open field build-up curve at 100 cm SPD for a 


6 MV beam. 
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38) 


energies have been compared by normalizing the depth with respect to 
ae of each curve. The result is illustrated in Figure 14. At 

field sizes of 10 cm x 10 cm and 30 cm x 30 cm, the TMR curves for the 
6 MV photon beam have a similar form to those for the 15 MV beam at 
depths less than daa However, at any normalized depth the TMR 
values for 6 MV slightly exceed those for 15 MV. The 6 MV surface 


dose exceeds that of the 15 MV beam for the same field size. 
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Zz.2 Accessory Build-Up Curves 


Various thicknesses of Lucite (density = 1.18 g/cm ) alone or 
Lucite and lead (density = 11.4 g/cm”) slabs were placed in the beam to 
determine their effect on the build-up curves. The distal surface of 
the slabs was placed at the accessory tray holder position 56 cm from 
the source. In all cases, a Lucite surface was facing the probes 

Figures 15, 16 and 17 show the build-up at an SPD of 1408 100 and 
joc, respectively. The field s#ze-in’all cases is 30: cm-x 30° cm: 
When accessories are placed in the field at the tray holder position, 

a larger TMR for all depths less than d ax is observed. There is 
also a shift in the point of cae to shallower depths. This is more 
pronounced at SPD = 75 cm (Figure 17). There is virtually no 
difference in the build-up curves between a 3.2 cm (1-1/4'') Lucite 
accessory and a 0.30 cm lead slab on top of a 0.64 cm (1/4"') Lucite 
tray. 

At smaller source-to-probe distances there is a greater difference 
between the tissue maximum ratio with and without accessories present. 


Therefore, as the distance to the source of contamination decreases, 


the dose increases. The tissue maximum ratio for SPD = 75 cm (Figure i7) 


is greater than the tissue maximm ratios for SPD = 100 cm and SPD = 
120 cm (Figures 15 and 16) at all depths. The difference between 
the tissue maximm ratios as a function of source-to-probe distance is 
most pronounced at the surface. 

The accessory build-up curves were measured on the Mevatron -6 for 
SPD = 75 cm and SPD = 100 cm. They show very similar features when 
compared to the 15 MV curves (see Figures 18 and 19). Accessories 


placed in the field increase the TMR at all depths below den There 
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Figure 15. 
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Accessory build-up curves at 140 cm SPD for a 15 MV beam. 
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Figure 16. 
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Accessory build-up curves at 100 cm SPD for a 15 MV bean. 


ay, 


i 


: > 
ed 
Bayt 


+ is a EOS see ee 


: ee & 
Redeye wa ya nh Awe 
Pe ary 7 i 3 7 ; ; or " : . 


Figure 17. 
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Accessory build-up curves at 75 cm SPD for a 15 MV beam. 
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Figure 18. Accessory build-up curves at 75 cm SPD for a 6 MV bean. 
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Figure 19. Accessory build-up curves at 100 cm SPD for a 6 MV bean. 
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MomOrdiiterence between a 3.2 cm and a .64 cm*thick lucite accessory 

at an SPD of 100 cm. However, a small difference is observed in TMR 
between the same two thicknesses of Lucite when introduced at SPD = 75 cm. 
The depth at which she occurs is shifted from 1.6 cm to 183,em with 

the introduction of the Lucite accessories at both SPD = 75 cm and 

100 cm. A comparison between 6 MV and 15 MV was made by normalizing 

the depth to dant The comparison is illustrated in Figures 20 and 

21. As observed for the open field build-up curves, the accessory 
build-up curves at 6 MV exceed those at 15 MV at the same relative depth. 


This applies for both source-to-probe distances of 75 cm and 100 cm. 
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Figure 20. Comparison between 6 MV and 15 MV accessory build-up curves 
at 7 5wcm SPD. 
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Figure 21. 
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Comparison between 6 MV and 15 MV accessory build-up curves 
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2.3 Filtering the Electron Contamination 


When at the surface, the chamber used has a sufficiently thin 
window to ensure that the build-up depth is approximately zero (5 um 
polystyrene equivalent). The window thickness has an equivalent thick- 
ness of less than a centimeter of air. Only very low energy Snehane 
would have had a reasonable probability of interaction with such a 
window. In order to be detected these photons would have had to have 
been produced within a few centimeters of the chamber, otherwise they 
would have been rapidly attenuated in air. Primary photons with a low 
enough energy to have had a high probability of interaction with the 
window would have been completely attenuated when the beam emerged from 
the beam defining head. Therefore, the surface dose is not due to low 
energy forward directed photons. Some of the dose at the surface is 
attributable to backscatter, but most of the dose is due to contaminant 
electrons. 

Filtering of electron contamination produced by accessories is 


ot The 


rarely practiced for megavoitage energies above those of Co 
efficacy of electron filtering for materials of various atomic numbers 
at 6 MV and 15 MV was studied. 

Placing a filter beneath an accessory produces a number of effects; 
the intensity of the primary beam is reduced, the number of interactions 
to produce electron contamination is increased and there is increased 
attenuation and scattering of electrons. The reduction of the primary 
beam can be taken into account by normalizing the readings taken at the 
surface to readings taken at oleae under the same conditions of filtration; 


in other words, by taking a tissue maximum ratio under conditions of 


filtration. 
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Figure 22. Reduction in the 15 MV surface dose when electron filters 
with materials of various atomic numbers are placed beneath 
the accessory tray. In all cases, an equilibrium thickness 


Of filter was used. 
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It is unlikely that filtration changes the amount of electron contami- 
nation present in the beam because the amount of electron contamination 
did not continue to increase with thickness of accessory after 0.64 cm 
ef Lucite at 15 MV (see Section 2.2 and Figure 16). This was verified 
by placing the filter material above the accessory tray (thickness = 
0.64 cm of Lucite). The reading at the surface with the filter material 
above the accessory compared to day was no different within experimental 
error from that with the accessory tray in the field alone. 

There was some reduction in the surface TMR when a filter with an 
atomic number higher than that of Lucite was placed beneath the accessory 
tray. The ratio of the TMR of the accessory tray alone to the TMR with 
a filter underneath the tray is shown in Figure 22 for 15 MV at field 
Sizes Of LO0icm x 10jcm and 15 can x 15 cm with an SPD of 100 cm. High 
atomic number materials produce the most filtering although the effect 
is not very significant. Lead will reduce the surface dose by about 
20%. An increased amount of electron scatter and absorption can be the 
only cause of the filtering effect. 

Similar results were obtained at 6 MV for a field size of 15 cm x 
15 cm at an SPD of 100 cm (see Figure 23). Again, lead proved to be 
the mee effective filtering material at this energy. It reduced the 


surface dose by about 20% as well. 
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Figure 23. 
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Reduction in the 6 MV surface dose when electron filters with’ 


materials of various atomic numbers are placed beneath the 
accessory tray tor a field*size of 15 cm-x [5 cm. 
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2-4 Distance Distribution of the Surface Dose 


aay | 


LeAne 


SPD Dependence of the Surface Dose 


The surface maximum ratio is defined to be the tissue 
maximum ratio measured with the chamber at the surface 
(i.e. d=0). Figure 24 illustrates that the surface maximum 
ratio in the open 15 MV beam increases rapidly as the source- 
to-probe distance decreases for both 20 cm x 20 cm and 
50 cm xX 30 cm field sizes. 

The curves for field size collimated to 20 cm x 20 cm 
show the effect of a thin accessory. The Siemens Mevatron 
-20 is equipped with a removable 1 mm thick Lucite cross- 
hair tray which can be inserted at 40 cm from the Shee 
When the source-to-probe distance is less than 75 cm, the 
surface dose with the cross-hairs in place is greater than 
when it is removed. However, the cross-hairs, when in 
place, reduce the surface dose in the clinical region of 
SPD greater than 75 cm. This is an exception to previous 
observations that accessories in the field increased the 
surface dose. 


Determination of the Apparent Source of Contamination 
Evéectrons 


As the source-to-probe distance increases, the volume 
60 


of air between the target and probe increases. For Co, 
Nilsson and Brahme (46) have predicted an increase in 
absorbed dose as a function of source-to-probe distance 


i! photon interactions with air. In these 


due to Co° 
experiments the reverse appears to be the case anda 


decrease in the dose as a function of increasing source- 
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Distribution of the surface dose as a function of source-to- 


probe distance at 15 MV. 


Figure 24. 


to-probe distance has been observed. This implies that at 
L5-MV, tor avtieldesize o£-30 cm x 30 cm, the interaction 
of the primary photon beam with air is not the major source 
of electron contamination. 

The surface maximum ratio is normalized to the dose 
measured at ee within a polystyrene phantom. The 
maximum dose in a phantom is due almost entirely to primary 
photons for which the variation with distance falls off 
as 1/ (SPD). However, the electron contamination does 
not necessarily arise at the source of primary photons. 
Instead of normalizing the surface dose at a given SPD to 
the dose at d ax? each surface dose at the same SPD is 
normalized to the dose measured at a convenient point 
inside a phantom using a standard field size and photon 
source-to-probe distance. The isocentric normalized dose, 
IND, is defined as: 

IND = Dose delivered to a point at an arbitrary d, SPD, field size 
Dose delivered at the isocenter 
(22555) 
The isocenter point was at che ateorl.=s LOO; cm with aster 
Sizer0r 1 0remn x 10s. 
The apparent source of contamination electrons can 
be located assuming that it is a point source and there is 
no attenuation of electrons by air. Therefore, the isocenter 
normalized dose at the surface is directly proportional to 


the inverse square of the distance to the source of the 
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IND (d=0)« ————,, 
(SPD-d.) 


CAD 


where qd, is the distance of the source of contamination 
electrons below the primary photon source. 

Figure 25 is a graph of V/1/IND(d=0) versus SPD. The 
graph yields a straight line except for very large source- 
to-probe distances. Therefore, the source must be small 
and the inverse square approximation valid for source-to- 
probe distances less than 120 cm. The x-intercept of the 
graph yields a distance, do of 12 cm below the primary 
source of photons. This corresponds approximately to the 
position of the bottom of the fixed head assembly which 
consists of the field flattening filter and the beam 
monitor ion chamber. The bottom surface of the field 
flattening filter is 3.1 cm in width, which at distances 
greater than 50 cm, subtends an angle no larger than ings 
which would approximate a point source. 

The dependence of the surface isocentric normalized 
dose, IND(d=0), on SPD, was investigated for a 6 MV beam. 
The inverse root of the isocentric normalized dose, 
1/VYIND, versus SPD is plotted in Figure 26. The intercept 
appears to be between 30 and 40 cm beneath the target, 
however, the line is not straight (in contrast to that 
obtained at 15 MV) indicating that a l/r? dependence from a 
localized source of contamination is not the entire explana- 


tion for the presence of the surface dose. 
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Graph illustrating the inverse square dependence of the 
surface dose on the distance to the source of contamination 


at 15 MV. 


Figure 25. 
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Figure 26. 
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Graph illustrating the failure at 6 MV of the inverse square 
dependence to a localized source of contamination. 
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2.5 Lateral Distribution of the Surface Dose 

The magnitude of surface dose as a function of lateral displacement 
from the central axis was measured for several source-to-probe 
distances. The central axis angle, C.A.A., was defined to be the angle 
between the central axis and a line joining the probe to the apparent 
source of electrons. This established parameters for the lateral 
direction in order to compare the distribution at various source-to- 
probe distances (see inset on Figure 27). In each case, the surface 
dose at any central axis angle and source-to-probe distance is 
normalized to the dose at d ax on thecentral axis (ClA.A,=0). The 
curves in Figure 27 follow a Gaussian distribution. 

At all points across the field, a smaller SPD results in a larger 
electron contamination dose. This agrees with the surface dose 
dependence measured along the central axis. Both the lateral dis- 
tribution and the dependence on the distance from the source of the 
electron contamination agree with the work of Almond (5). 

The lateral surface dose distribution for 6 MV and 15 MV beams at 
SPD=100 cm is shown in Figure 28. The distribution for each beam has 
been normalized to the maximum surface dose. Since the distribution 
at 6 MV was mapped for only one source-to-probe distance the abscissa 
is left as the distance from the central axis. The ordinate is 
defined as the surface dose at some distance away from the central axis 
normalized to the dose at ce measured at the central axis. There are 
two differences between the 6 and 15 MV lateral distributions. First, 
the maximum dose does not occur at the central axis but from 4 to 6 cm 
Ditethe? central, axis... fhese, !’horns''}in the.surface, dose, are a very 


small but real effect. Secondly, there is a more rapid decrease in 
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The 


inset is a schematic diagram of the experimental arrangement. 


Gaussian distribution of the surface dose at 15 MV. 


migure 27. 
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Comparison between 6 MV and 15 MV of the lateral distribution 
Olathe surtace dose at. 100 cn: SPD...” The lateral distribution 
isnot Gaussian at 6 MV.. The field size is 50° cmix. 30) cm: 
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the surface dose of the 6 MV beam near the field edge compared to the 


15 MV beam. The lateral distribution at 6 MV is not Gaussian. 
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2.6 Surface Dose Measured for Rectangular Fields 

The surface isocentric normalized dose was measured for rectangular 
and square fields. Both collimator openings were varied by 5 cm 
increments at 100 cm SPD so that the dose at all such fields between 
5 cm x 5 cm and 30 cm x 30 cm were determined. Table 2 illustrates the 
FESULTS. 

The surface isocentric normalized dose as measured at SPD = 100 cm 
depends on which collimator defines the long or short axis. If Field A 
(upper collimator) is the long axis (lower left hand part of Table 2) 
the surface dose tends to be greater than if Field B (lower collimator) 
is the long axis for the same field dimension. For example, the 
surface isocentric normalized dose at A = 10 cm and B = 25 cm is .106 


whereasuthe doseat’ Bs = 10.cmand A=125 cm as 2109. 


In general, the larger the field area, the greater the surface 
dose. If a square and rectangular field have the same area the 
square field will have the greater surface isocentric normalized dose. 


For example, if A = 10 cm and B = 10 cm, the surface dose is .071 but if 


> 
I 


encima B= 207cm «the surtace doseuis...065 (12.A.=920)cm-and 


ee) 
ll 


s7cem,, the dose is 064), 

The equivalent square field dose for a rectangular field was 
determined by plotting the dose as a function of square field width 
and finding the square field that has the same dose as the rectangular 
field. Table 2 also lists the results of the equivalent square 
determination. For large fields (greater than 10 cm x 15 cm) the 
equivalent square field area is approximately equal to the geometric 
field area. For example, if A = 20 cm and B = 10 cm, the equivalent 


2 
square area is (13.8 + 4)cm x (13.8 # A)cm = (190 + 10)cm and the 
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peometricsarea 15 (20.0 21) cnx (10-0e4+-.1) on = (200 4 3) cm’. 
However, for small fields the equivalent square field area is less than 
the geometric field area. The surface equivalent square field is very 


different from the equivalent square field at ak 
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2.7 Penetration Curves Outside the Field 


2.7.1 Transmission Through the Collimators 


2. Med 


The transmission of primary photons through the 
collimators was studied by placing a Capintec 0.6 cc 
cylindrical ion chamber (Model PR-06C) at ea A measure- 
ment was taken at the central axis with a field size of 
20 cm x 20 an. Ail subsequent readings were normalized 
to this measurement. 

Transmission measurements were made with one or both 
collimator sets closed. Figures 29a) to d) illustrate 
the position of the upper and lower collimators and the 
relative dose received. The percentage transmission at the 
field boundary through one collimator is about 1% and 
about 0.1% through both collimators. 

Determination of the Penetration Curves Outside the Field 

By definition, except for transmission through the 
collimators, the primary photon dose outside the field 
boundary is zero. Only contamination electrons and 
scattered photons can contribute to this dose. This was 
studied experimentally. 

Nilsson and Brahme have shown that the scattered photon 
spectrum at the field edge is almost identical to the 
photon spectrum at the central axis for a 20 cm x 20 cm 
field for nominal beam energies of 6 MV and 21 MV. It is 
unlikely that the electron contamination spectrum changes 
appreciably between the center and edge of the field. 


In order to get as close as possible to the field 
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Figure 29. Transmission measurements on the Mevatron -20. a) Lower 
collimator closed, upper open; b) Upper collimator closed, 
lower open with collimator assembly rotated 90°; c) As in b) 
but collimator assembly not rotated; d) Both collimator sets 
closed. 
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boundary, the thin-window parallel-plate probe was placed near the edge 
of the phantom (see Figure 30). The phantom was tilted at an angle of 
8.5°*% so that the position of the phantom surface would not move closer 
or further from the field boundary with the addition of build-up layers. 
This also reduced the number of particles incident at oblique angles 
upon the phantom at oblique angles. 

Figure 31 illustrates the penetration curve 3 cm outside a 
30 cm x 30 cm field. The penetration curve has a peak at a shallow 
depth (1-2 mm) then falls off rapidly at depths greater than 1 cm. The 
curve has a relatively slow fall-off at depths beyond about 2 cm. The 
general shape of the curve agrees with the work of Scrimger and Kolitsi 
(30) for 8 MV x-rays. 

The penetration curve 3 cm outside the field was measured at 6 MV 
for the same field size (30 cm x 30 cm) and SPD (100 cm) as the 
comparable experiment carried out at 15 MV. The phantom was again 
tilted at 8.5°. Figure 32 illustrates a rapid decrease in the penetration 
curve at shallow depths with a levelling out at greater depths. The 
maximum dose is at the surface at 6 MV rather than at a shallow depth 
beneath the surface. This suggests the dose at shallow depths is again 
due mainly to contamination electrons while the dose at deeper depths 
is due to scattered photons. 

The relative dose in the plateau region (dad ox) for both che 167 MV. 
and 15 MV curves outside the field is approximately 0.3 of the maximum 
dose received. 
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Figure 30. Diagram of the experimental arrangement to determine the 
penetration curve outside the field. 
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Figure 31. Penetration curve at 15 MV obtained 3 cm outside of a 
30 cm x 30 cm field (see Figure 30 for the experimental 
arrangement). 
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Figure 32. 


ml 
© 
= 
O 
= 
Oh 
O 
ee) 
= 


oe top Oe SAN RNS 


(ay) Oude © WNUWIXD\A) SNSsl | 


Penetration curve at 6 MV obtained 3 cm outside of a 
SUSE] SOc cmet 1 eld. 


Li) 
j 


; 


he 


a2 


oe, 


87 


2.8 Increase in Dose with Increasing Field Width 


Weel 


Zesuz 


Measurement of the Stem Effect 

Exposure to ionizing radiation may cause a spurious 
Signal attributable to the ion chamber cable such as a 
leakage current which is mistaken for ionization in tie 
chamber volume. These phenomena are called the stem 
effects. A stem effect, if present, should increase with 
increasing field size. The manufacturer's claim that the 
stem effect for their thin window ion chamber is less than 
1% 

The stem effect was measured by exposing the chamber 
at chee in a phantom to a rectangular field as wide as the 
chamber volume and as long as the phantom. One measurement 
was taken with the cable completely in the field and another 
with the cable completely outside the field. The stem 
effect was found to be less than 0.5%. 

Increase in Dose with Increasing Field Width 

A number of authors (15-21,25,28,29,31-35) have 
commented on the dependence of the surface dose on the 
square root of field width both with and without accessories 
present in the beam. Many of these authors calculated TMR 
surface doses, in other words, normalized the surface dose 
at each field size to the dose at da for the same field 
size. However, since the dose at CR, increases with field 
size (as expressed by the output factor), such a plot would 
not unambiguously express the dependence of surface dose on 


field size. 
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A more general form of the isocentric normalized 
dose can be used to more clearly describe the dose dependence 
on field size. The ionization readings at some field size 
and a given source-to-probe distance are normalized to the 
maximum ionization readings (at de obtained at a field 
size of 10 cm x 10 cm at the same source-to-probe distance. 
This quality will be called the normalized dose, ND. It is 


given by: 


Dose) delivered.to,. a,point.at.an,. arbitrary d, SPD,. 

2 field size 
~ Dose delivered to d ax at the same SPD at a field size 
i Orel 0. cma LL Usem 


(Eyes s) 


The normalized dose is equal to the isocentric normalized 
dosemat 100¢cm SPD: 

The tissue maximum ratio is directly proportional to the 
normalized dose. The constant of proportionality is the 


output factor. The relationship is given by: 


ND (d,SPD,W) = Output Factor (SPD,W) x TMR (d,SPD,W) 
ears) 
where W represents the field size. 
The normalization reading (the denominator of Equation 2.8.1) 
is not field size dependent so ND clearly illustrates 
dependence of the reading (the numerator of Equation 2.8.1) 


on field, siize.* 


* Data acquisition is simpler using the concept of normalized 
dose. Only one set of measurements at dmax for a field 
size of 10 cm x 10 cm need be taken. A set of readings at 
Graxsatsalletield sizes is not required. If an output factor 
curve has been determined, the TMR value may be obtained, if 
desired, using Equation 2.8.2. 
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The normalized dose, ND, for a variety of conditions at 
d=0 is plotted as a function of field width in Figure 33. 
Field width is defined as the square root of field area 
for square fields. Under these conditions, at 15 MV 
there is a linear increase in the surface dose with 
respect to field width. Note that either decreasing the 
SPD or placing the accessories in the field increases the 
slope of the curve. Since the surface dose for any given 
field size is proportional to the slope of the line, the 
Slope is also a measure of the electron contamination, 
but with the advantage that it is not dependent on the 
field width. 

A small phantom was placed compietely within the 
field so that an increase in detected signai due to in- 
creasing field size could oniy be attributed to beam con- 
tamination. To obtain the dose as a function of depth for 
as many field widths as possible, the phantom lateral 
dimensions were kept as small as possible (i.e. 10 cm x 
10 cm) which in turn dictated a minimum field size used in 
the determination of the contamination depth-dose curve of 
15 cm x 15 cm. The maximum square field size attainable 
for the machine studied was 30 cm x 30 cm. These field 
sizes defined the range of the study. 

The dose as a function of field width was measured 
for a number of depths with the limited polystyrene 
phantom of the type illustrated in Figure 34. Since 
there was build-up material on top of the ion chamber, it 


could not be assumed that ali of the contamination 
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Linearity of the surface dose with respect to field width 


under a variety of conditions at 15 MV. 


Figure 33. 
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Figure 34. 


10 cm SCALE 


Diagram of the experimental arrangement to measure the 
increase in contamination with increasing field width. The 
phantom slab cross-section was 100 cm2 and the slabs were 
always within the field boundary. 
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contributing to the dose was due to electrons. 

Figure, 55 illustrates that at 100 cm SPD for various 
depths in the phantom less than oe there is a linear 
dependence of normalized dose on field width. Qualitatively, 
the slope at shallow depths is greater than the slope at 
deeper depths. As has been discussed for the dose at the 
surface, the slope of normalized dose versus field width 
is a measure of the magnitude of the contamination 
penetrating to the depth specified. Therefore, there is a 
greater magnitude of contamination at shallow depths than 
at deeper depths. 

Figure 36 is a graph of normalized dose versus field 
width at 75 cm SPD for various depths into the limited 
phantom. The linear fit at 75 cm SPD is nearly as good as 
that obtained at 100 cm SPD. 

Normalized dose measurements were made on the 15 MV 
accelerator at#field sizes of 15 cm*x"15 cm, 20 cm x-20 cm 
and 25 anexe25 cmeat -l00°cm SPD when Lucite accessories of 
0.64 cm and 3.18 cm were placed in the field at the accessory 
tray position. There is approximately a linear increase in 
dose as a function of field width with a 3.18 cm accessory 
in place (see Figure 37). 

Figure 38 is a graph of normalized dose versus field 
width for various depths at 100 cm SPD for the 6 MV 
accelerator. Like the curves at 15 MV, the linear fit is 


good. 
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Figure 35. 
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The increase in normalized dose with respect to field width 
for various depths beneath the surface of the phantom in an 
open field at 100 cm SPD for the 15 MV beam. The increase 

can only be attributed to contamination because the phantom 
was completely within the field boundary (see Figure 34). 
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Figure 36. The increase in normalized dose with respect to field width 
for various depths beneath the phantom surface in an open 
Pieidcat./SecmeobD ati. 15. MV, 
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The increase in normalized dose at 15 MV with respect to 
field width for various depths beneath the phantom surface 
with the probe at 100 cm from the source with a 3.18 cm 
accessory in place. 
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Figure 38. The increase in normalized dose with respect to field width 
for various depths beneath the phantom surface in an open 
field at 100 cm SPD at 6 MV. 
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2.9 Determination of the Amount of Photon Backscatter from the Phantom 

The surface dose is the sum of backscatter from the phantom and 
contamination produced outside the phantom. By the definition of 
phantom generated radiation (Section 1.5), backscatter is part of this 
component of the beam. Backscatter photons are any photons scattered 
through 90° or more so that they have a sense in the direction component 
normal to the phantom that is opposite to the sense of the primary 
photons. A determination of the amount of backscatter dose is essential 
in order to isolate the dose due to contamination at the surface. 

The amount of backscatter due to photons generated in the phantom 
was measured. The thin window parallel-plate chamber was first placed 
on a thick siab of polystyrene (>20 cm thick) which had a cross-section 
larger than any field size used. The surface dose was measured for 
field sizes between 3.5 cm x 3.5 cm and 20%cm % 20:cm and source-to- 
probe distances between 80 cm and 120 cm. The normalized dose, N.D. 
(Equation 2.8.1) was obtained. 

The full phantom was replaced with various thicknesses of back- 
scatter. The thinnest layer was approximately 1.5 mg/cm? and the 
thickest layer was 480 mg/cm? of unit density material. There was 
virtually no difference between the readings at 1.5 mg/cm? and 480 mg/cm? 
indicating that there is little contribution from charged particles set 
in motion near the surface of the phantom which scatter backwards. The 
normalized dose difference was obtained between the normalized doses 
with full phantom and with 1.5 mg/cm’. The normalized dose difference 
represents the normalized dose due to backscatter photons. The results 
are shown in Figures 39 and 40 as a function of field width for 6 and 


15 MV at various source-to-probe distances. 
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Figure 39, 
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The normalized dose due to backscatter photons as a function 
of field width for source-to-probe distances of 80 cm and 
120 cm at 6 M. 
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Figure 40. The normalized dose due to backscatter photons as a function 
of field width for source-to-probe distances of 80 cm, 90 cm 
and 120 cm at 15 MV. 
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Between field widths of 5 cm and 20 cm, the amount of backscatter 
increases with field width for all source-to-probe distances. The 
amount of photon backscatter at 6 MV is about 50% greater than at i5 MV 
for the same field width and source-to-probe distance. At all field 
sizes and source-to-probe distances for both accelerators, the amount 
of dose due to photon backscatter was much less than the amount of 


dose due to contamination electrons at the surface of the phantom. 
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2.10 Sweeping Electrons from the Field Using a Permanent Magnet 


2.10.1 Determination of the Field Distribution of the Magnet 


ZaLOVeZ 


A survey of the literature (41, 42, 43, 44, 45) has 
shown that the central field magnetic strength must be 
somewhere between 1 and 3 kG in order to sweep a 10 cm x 
10 cm field nee of contamination electrons produced by 
photons in the energy range of 4 MV to 25 MV when the 
source-to-probe distance, SSD, is about 85 cm. 

The magnetic field strength distribution of a dis- 
carded magnetron magnet was determined using a Rawson- 

Lush Rotating Coil Gaussmeter Type 820 Probe and a Type 
501 Indicator manufactured by Rawson Electrical Instruments 
Co. of Cambridge, Massachusetts. 

The uncertainty in the strengths at maximum field 
strength was + 30 G. The magnet was a permanent type with 
the pole gap fixed at 7.5 cm. The maximum central field 
was found to be 1340 G. Figure 41 illustrates the measured 
magnetic field distribution. The central distribution 
was mapped every 2 cm until the magnetic field fell to 


60 G. 


Verifying Electrons are Swept by the Magnet 


If the amount of bending is small, to a first 
approximation, the electrons will travel along the central 
axis in and near the magnet. The angle of deflection 
(in radians) for an electron travelling a distance, ay 
in a uniform magnetic field is given by (see Figure a2) 
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Figure 41. 


The measured magnetic field distribution along the direction 


of the photon beam. 
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Constant Magnetic Field 
Direction Out Of rhe Plane 
if d Is Small, d-d 


Figure 42. The angle of deflection of an electron in a uniform magnetic 


field (above) and the measured deflection of 18 MeV electrons 
(below). 
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The above equation uses the definition of a radius of 
Curvature. (Equation 1.4.1). The total amount of 
deflection can be found by summing the amount of angular 


deflection occurring along the central axis: 
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The field, B; , WaS measured every 2 cm so d; is a constant 
(d, =const=2 cm). The Mevatron -20 an produce mono- 
energetic 18 MeV electrons. The electron contamination 
of the 15 MV beam will not have kinetic energies this 
high so that if the 18 MV beam can be deflected then the 
electron contamination can also be swept. The amount of 
deflection predicted by Equation 2.10.2 for 18 MeV 
electrons in the measured magnetic field distribution is 
approximately ye 

The magnet was mounted 65 cm from the source, so the 
largest field size (measured at SPD=106 cm) that can be 
swept is 10 cm x 10 cm. This defines a central axis 
angle of ot Clearly, this simple theory predicts 
enough deflection will occur. 

The analysis was verified semi-qualitatively by 
experiment. A sheet of film (Kodak XV-2) was placed 
parallel to a beam of electrons. An equal amount of 


exposure was delivered with an open electron applicator 
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with and without the magnet beneath the applicator. 
Figure 42 illustrates the results. The angle of 
deflection between the image of the beams is about 20°° 
The deflection of contamination electrons was 
directly verified by placing a half cylinder of filn 
beneath the magnet. Figure 43 illustrates the setup. 
The contamination electrons produced a broad hazy patch 
on the film (see Figure 44). The width of the image 
increased as a function of angle swept. This is to be 
expected because the amount of bending and the amount 
of scatter eee with decreasing electron energy. 
The overexposed central patch is mainly due to photons. 
Transmission through the collimator can be discerned. 
There is no direct evidence of positron contamination 


which would be deflected in the opposite direction. 
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Figure 43. The set-up to verify that contamination electrons were 


being swept from the field of the 15 MV beam. 
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Figure 44. Image of the swept electron contamination. 
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The Surface Dose With and Without the Magnet in Place 


The magnet was mounted on a holder so that it could be 
placed just beneath the accessory holder position. The 
holder had a square hole cut in it enabling a 10 cm X 
10 cm beam to pass through without interacting with it: 

A photograph of the magnet attached to its holder is shown 
in Figure 45. Figure 46 shows the magnet mounted from the 
accessory holder position. 

With the magnet in place the electron contamination 
produced between the target and magnet is swept from the 
beam. The surface dose under this condition is due to 
electron contamination produced between the magnet and 
phantom and photon backscatter produced from the phantom. 
Since air is the only material in the beam between the 
magnet and phantom, the electron contamination with the 
Magnet in place is produced by interactions of the photon 
beam with air between the magnet and phantom. 

The dose was measured at the surface of a full phantom 
(>20 cm thick and with a cross-sectional area greater than 
the width of the beam). The measurements were taken 
under various conditions of field size and SPD and presence 
of accessories. The small magnet pole separation meant 
that the target field could be no greater than 10 cm x 10 cn. 

Figure 47 illustrates the dependence of the normalized 
dose at the surface, when the magnet is in place, on field 
width for various accessory configurations at an SPD of 


100 cm. The maximum doses used to normalize these curves 
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Figure 45. The magnet attached to its holder. 
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Figure 46. The magnet mounted from the accessory holder position. 
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Figure 47. 


Ox 
S 
> 16 " 
= palates 
= 
53 2 So 
Tes Sac 
YD + — 
Cale (s 
U Te) = 
£ g O 
a = 
al Mea) 
Bey = 
oO Va) 
a a 
c \) 
Orsay 
a O = 
ore 
4q 
—— ae 


S&S > 


eee 1 riley Ole Ne) 
JSOd GAZITVWYON 


. 


The surface normalized dose with the magnet in place for 
various conditions at 100 cm SPD for the 15 MV bean. 
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were Measured at an SPD of 100 cm at a depth of d ax for a 
field size of 10 cm x 10 cm under the same accessory con- 
figuration. The normalized dose curves are the same 
within 0.005. This is to be expected because electron con- 
tamination produced by the accessories is swept from the 
field by the magnet. The normalization procedure takes 
into account the reduction in the primary dose available for 
the production of electron contamination in air and photon 
backscatter. 

The dose due to photon backscatter can be subtracted 
from the surface dose with the magnet in place to leave 
the dose due to electron serniet or produced in air* 
as the only component. Figures 48 and 49 illustrate the 
normalized dose due to contamination produced in air as a 
function of field width for various source-to-probe 
distances when an open field is used (with cross-hair tray 
in place) for a 15 MV and 6 MV photon beam, respectively. 
At each source-to-probe distance the normalized dose 
increased as a function of field width. This increase was 
larger for smaller field widths. 

The measurements suggest that curves taken at various 
source-to-probe distances converge to a zero normalized 
dose at a zero field width. The behavior at field sizes 
less than the chamber diameter cannot easily be determined. 

The difference between readings with and without 
the magnet in piace can only be due to what is being 


swept by the magnet when it is in place; namely electron 


* This refers to air between the magnet and phantom. 
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Figure 48. 


6 / 


5 
FIELD WIDTH (cm) 


4 


3 


(oe) ON —~ 
eee OO Oe ae 


3900 CazivWWYON: 


Normalized dose due to contamination produced in air as a 
function of field width for various source-to-probe distances 
for an open field at 15 MV. 
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Normalized dose due to contamination produced in air as a 
function of field width for various source-to-probe distances 


for an open field at 6 MV. 
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contamination produced between the target and magnet. * 
Figures 50 and 51 are graphs of the normalized dose 
difference between magnet in and out for 15 MV and 6 MV 
photon beams, respectively. The normalized dose is very 
small at small field sizes, but increases rapidly as a 
function of field size. The normalized dose due to 
contamination produced between the magnet and target in- 
creases approximately in direct proportion to the beam 
cross-sectional area. 

There is considerable difference between the electron 
contamination contributions due to air and due to inter- 
actions between the target and phantom at small field sizes. 
(less than 10 cm x 10 cm). The contributions due to air 
exceed those of contributions produced between target and 
magnet for both 6 MV and 15 MV photon beams at all source- 
to-probe distances. The increase in duse as a function of 
field width at very small field sizes (less than 5 cm x 
2 cm) is greater for electron contamination produced in air 
than for electron contamination produced between the target 
and magnet. The most interesting difference between these 
Curves is that the normalized dose for air-produced 
electron contamination is greater for large SPD whereas 
the target-to-magnet electron contamination is greater for 


smaller source-to-probe distances. The source-to-probe 


* Tt did not matter if normalization measurements were 
done with the magnet in or out. The difference at dj,, 
between readings taken with the magnet in place compared 


to readings with the magnet not in place was less than 
0.5% This indicates that few contamination electrons 


produced between the target and magnet penetrate to dj,,. 
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Figure 50. Normalized dose difference between magnet in place and out 
of place for a 15 MV beam as a function of field width for 
various source-to-probe distances for an open Pied 
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Normalized dose difference between magnet in place and out 
of place for a 6 MV beam as a function of field width for 


various source-to-probe distances for an open Piel 


Figure 51. 
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dependence of the surface dose is better illustrated in 
Figures 52 to 55. Figures 52 and 53 are graphs of normalized 
dose versus SPD at 15 MV for air contamination and contamina- 
tion between the target and magnet, respectively, for a 
family of field widths. Figures 54 and 55 are for contamina- 
tion produced in air and between the target and magnet, 
respectively, under the same measurement conditions as 15 MV 
but done at 6 MV. 

If an assumption is made that the location where the 
contamination is produced between the ieee and magnet is 
a point source, then the same procedure used in Section 
2.4.3 to locate the apparent source of contamination 
electrons can be used again. The difference in the absolute 
readings taken with the magnet not in place and magnet in 
place is a measure of the contamination dose produced 
between the target and magnet. If the readings were taken 
at nearly the same time, then the difference in absolute 
readings need not be normalized. Figure 56 is a graph of 
the square root of the difference in absolute readings 
versus SPD at 15 MV for a number of field sizes. The 
apparent source of electron contamination is at 40 cm SPD. 
The electron contamination appears to be coming from the 
cross-hair tray. A similar result was obtained at 6 MV 
(see Figure 57). 

The cross-hair tray was removed and the readings 
repeated at 15 MV. Figure 58 illustrates that the apparent 


source of contamination is shifted to between 30 cm and 


fue bene sau he Beer es : 
bostivavon lo ceynrk ewe Ge threes 3 vol a oF Sc.¢ 
grvaptin: leer chpethekisa she 34 a ar 8 Pte ecm 
5) whet domgee meant Nk Paget te enced mat 
-einegnes vet see ES Rew We - swereelt eter bis? “T tinea 
1 bes Jogunt de sespeerg sd iaeaty: rie Fol riser ots 
‘ae 2! picibeas Sages" iad seute old “ohne cbesoogear: 3 
1 hese engh int sh 


now Lae TAL <2 iJ tan bitag aj 7 icy is 


amipn bos Seyves wl aero hanchund od: wid pn a 
nuitom @l feay oti aed att detia: ses tog 8 . 

EH Seine to sede arena ‘ent su et a | , 
ottoman add at iniadiied ort .tiragh nr i al 


saeiay cr ane a aan het 


ee > ar ‘J rt Pe RN : 
sidan tain is bie 


S is71.¢ at fix soni 88H 8g aft i al tl 
To Aneta ®, Bf ob en et ri 
! 4 


oy hei: < 
ae 


2 
y iL ‘ : : 
amin beet protons. tia, i] 


= 


J 


= ; : arp 4G 
5 cor vir rey i mel He 


oo Sal | ie 


> ? Me 8 
en ya e: 
of ne 7 + - : 7 


> 


Figure 52. 
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Normalized dose due to contamination produced in air as a 
function of source-to-probe distance for various field — 
widths for an open field at 15 MV. 
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Figure 53. 


1 d magnet in a 
Normalized dose produced between the target an 
15 MV beam as a Paneeien of source-to-probe distances for 
various field widths for an open field. 
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Figure 54. 
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Normalized dose due to contamination produced in air as a 
function of source-to-probe distance for various field widths 


for an open field at 6 MV. 
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Figure 55. 
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Normalized dose produced between the target and magnet in a 
6 MV beam as a function of source-to-probe distances for 
various field widths for an open field. 
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Figure 56. 


35cm xX 35cm 


Sem Scm 


Jom X /7cm 


lOcm X 10cm 


de 


50 ]OO 150 
SPD (cm) 


Dependence of the dose produced between the target and magnet 
on the inverse square of the distance from the apparent source 


at 15 MV when the cross-hair tray was in place. 
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Figure 57. Dependence of the dose produced between the target and 
magnet on the inverse square of the distance from the 
’ apparent source at 6 MV when the cross hair tray was in place. 
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Figure 58. 
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Dependence of the dose produced between the target and 
magnet on the inverse square of the distance from the 
apparent source at 15 MV when the cross-hair tray was not 
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35 cm. The range of x-intercepts corresponds approximately 


to the position of the collimators. 


Pin re oe ue me hs 


‘ie EML Ase plain 
swomatie aah ter 


| eOrt a 
Ve! i Lat ais 
<< ie ie : hy 


Mog 
2.10.4 Determination of the Penetration Curves of Electron 
Contamination Produced Between the Target _and Magnet 

Electrons produced by interactions between the target 
and magnet Should penetrate some distance into a phantom. 
The penetration curve can be obtained by measuring the 
difference between readings taken without and with the 
Magnet in place when the same thickness of overlying 
material is placed over the ion chamber. 

Figure 59 illustrates the central axis penetration 
curve for electron contamination produced by the 15 MV beam 
between the target and magnet as measured by a probe at 
100 cm SPD ina 10 cm x 10 cm field. “ This penetration 
curve has its maximum at about 0.15 cm in the phantom and 
a relative plateau in dose of about 1 cm. From 1 cm to 
3.5 cm the dose drops rapidly as a function of depth. At 
2.9 cm the dose deposited is only 20% of the dose deposited 
at 0.15 Gm 

An almost identical penetration curve was produced 
when the probe was 3 cm outside a field with a size of 
16 cm x 10 cm at 100 cm SPD for a 15 MV beam (see Figure 
60}. The phantom setup was similar to the inset in Figure 
30 except the phantom was tilted at Re 

Figure 61 illustrates the penetration curve for 
electron contamination produced between the target and 
magnet for the 6 MV beam. The shape of the curve is 
similar to the result at 15 MV except that the maximum 
dose occurs at the surface at 6 MV and this beam 
reaches its 20% level at 1.3 cm instead of 2.9 cm. Figure 


62 compares the penetration Curves for electron 


sit Ne 


p ered sk Mar 


‘ 
1'% 
co 


haa tse % 
ch 


F CoP wii 
i'd ceanr AONE Be 4283 aed suet TH a 
5) ro ace aT 

| jootendea of ASS SOT RIF tae ad? \ 
i tote? egelbhes teextet @ empattih led 

| : Re 
oat wit met oneiq sf angen 


rT 


~ a] 


wt une arte ww bagels el Losec : 
A 5 p 


wettatlk Ve -emgrt | 


ee a ae iT 
1 Ve beoulrta ime » PATIsete 102% OVIUD 
' a 
Log! sagret ety aceeted 

” —_ ity ie 

a 7 


Go ibs ~ aac = OGL 


ot A} uote In m pe esi far ov? ae’ : 3 
i Ai H ites FH S206) 1 cr ths: ~" avitalen 5 me 
; 
vh titers 20 cere teiegs «2p nh eect tt 2 Ac 
‘ ie ie ot a “ao 
tet 2 ety fy oft Yo eqy Wim ero bepe bea! ee) pdt a 
L' - ‘ie © ; 7 r i] ed i. . 
i eal fiw . 
i> a 7 r 


= af as 
erayt toi een m | ‘saline pi te os : 
| Rita: ce aa Rend = 


ate as Rip ot ay 
—_ si Cae 
meu 


ah ; cae aL 
oa “te 


ie 9 ? 


wr 
ae Gs 


ssh, 


~ 


sO 
‘CN 
& 
are 
< 
©. 
ees ® 
C) 
re ars 
© ‘ «> - @ & 
© cO ~O = ON 


Es0(] aA} 


Figure 59. Central axis penetration curve for electron contamination 
produced between the target and magnet by the 15 MV beam at 
100 cm SPD for a field size of 10 cm x 10 cm. 
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Figure 60. Penetration curve for electron contamination between the 
target and magnet measured 3 cm outside a 10 cm x 10 cm 
field at 100 cm SPD for a 15 MV beam. 
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Figure 61. Central axis penetration curve for electron contamination 
produced between the target and magnet by the 6 MV beam at 
100 cm SPD for a field size of 10 cm x 10 cm. 
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Figure 62. 
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Comparison between the 15 MV and 6 MV contamination 
penetration curves produced between the target and magnet 
in a 10 cm x 10 cm beam at 100 cm SPD. 


TS 


contamination produced between the target and magnet for 
6 MV and 15 MV by normalizing the depth of penetration with 
respect to the depth required to reach the 20% level. 


The curves are almost identical. 
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2.li A Method for Separating Contamination from Phantom Generated 
Penetration Curves 


The graphs of normalized dose versus field width 
allow one to separate the dose due to contamination from 
the dose due to photons interacting in the phantom. The 
total dose, T(d,SPD,W), is a sum of the phantom generated 
dose, P(d,SPD), and the contamination dose, Cidade SPDewW). 
The phantom generated dose is not a function of the field 
width, W, when measured with a limited phantom (see Figure 
34). The phantom generated dose is due to primary photons 
interacting with the phantom and scattered photons 
“generated by the primary interactions in the phantom. 

The phantom generated dose is due to the primary photon 


component of the beam, Mathematically. 


T(d,SPD,W) = P(d,SPD) + C(d,SPD,W) xa wlene 


Total Reading 
Standard Reading 


P(d,SPD) = see Due To Primary Photons 


Standard Reading 


where T(d,SPD,W) = 


c(d,spp,w) = Reading Due To Contamination 


Standard Reading 


The standard reading used is the same as in Equation 2.8.1. 


oT (d, SEDW) > Pen SeD) i; ate Sr) (2iates 
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Therefore, 


dT(d,SPD,W) = dC(d,SPD,W) meng 
Seca ere ee C2 )) 


where oot, is the slope of the line of normalized dose 
versus field width at a depth, d, source- to-probe 
distance, SPD, and field width, W. The normalized dose 
due to contamination, C(d,SPD,W), is found by integrating 
Eyuation 2.8.3. This makes an assumption that the 
nature of the contamination does not change as a function 


of field width, but the magnitude of contamination does. 


Mathematically, 
uy, 
ciaeeruW)) = 9). See (2.11.4) 
i 
If aT (d, SPD,W) is not a constant then, 


dw 


C(d,SPD,W,,W,) = T(d,SPD,W,)- T(d,SPD,W,) (2.11.5) 


dine ard DN) is a constant then, 


C(deSPD) Slope of Line (Aot.6) 


where the last equation does not have any dependence on 


field width. An assumption is implicit that dose increase 


of small field size is the same as dose increases for 


larger Sizes. 
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The slopes in Figures 35 to 38 were linear so 
Equation 2.11.6 can be used. Therefore, the amount 
of contamination is independent of field width. 
Table 3 lists the amount of contamination as a 
function of depth normalized to the maximum amount 
of contamination for the 15 MV beam. The penetration 
Characteristics of the contamination are relatively 
independent of the thickness of accessory and source- 
to-probe distance. Therefore, to a good approximation, 
the relative amount of contamination only depends on 
the depth of measurement in the phantom. A plot of the 
contamination penetration curve at SPD = 100 cm is 
shown in Figure 63. 

The open field penetration behaviour at 100 cm 
SPD was obtained for the 6 MV accelerator also by 
using the increase in normalized dose versus field 
width. The normalized contamination penetration 
curve is shown in Figure 64. 

The contamination penetration curves are 
equivalent to the tissue maximum ratio of the con- 


taminant component if they were the only component in 


the beam. 


Determination of the Phantom Generated Build-Up Curves 


The phantom generated build-up contribution is 


obtained by subtracting the contamination component 


from the total build-up curve. The total build-up 


curve is measuréd with a full phantom (that is, the 
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Table 3 


Maximum Contamination Dose At 15 MY 


SPD = 100cm 


No 
Accessory 


SPD = 190@cm 
8.64cm Lucite 


Accessory 


Contamination Dose Normalized To The 


SPD = 100cH 
3.18cm Lucate 


Accessory 


SPD 
No 


Accessory 
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Figure 63. 
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Figure 64. Contamination penetration curve determined by the increase 
in dose in a limited phantom as a function of field width 


at 100 cm SPD for a 6 MV beam. 
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cross-sectional area is greater chen the field 
width). The magnitude of contamination under 
any condition is determined by equating the 
difference between the total build-up curve 
and the photon backscatter at the surface with 
the magnitude of the contamination dose at the 


surface. Mathematically, 


C(d=0,SPD,W) = T(d=0,SPD,W) - PBS(d=0,SPD,W) 


(2.1437) 
where PBS(d=0,SPD,W) is the amount of photon 
backscatter at the surface of the phantom. The 
amount of contamination at any depth is determined 
from a knowledge of the amount of contamination 
there is at the surface and the penetration 


characteristics. Mathematically, 


C(d,SPD,W) = C(d=0,SPD,W) Aa) (2.11.8) 


where CIMR(d) is the value of the normalized 
contamination penetration curve at a depth, d, in 


the phantom. 

Having established the magnitude of contamina- 
tion at any depth in the phantom, the contamination 
component is subtracted from the total build-up 


component to yield the phantom generated build-up 


curve. 
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A computer program called Seprimary.For* (for 
a listing and cociientapeen of Seprimary.For see 
Appendix 5) calculates the phantom generated build- 
up curves. The input data lists the normalized 
contamination penetration data as a function of 
depth, the photon backscatter normalized dose at 
the surface as a function ake field size as well as 
SPD and the total build-up normalized dose as a 
function of field size, SPD, depth and amount of 
accessories in the field. 

The "pseudo-unit'"' of the output is normalized 
dose. This can be converted to the tissue maximum 
ratios of the phantom generated build-up curves 
by dividing the curves by the maximum normalized 
dose value of each curve. 

An estimate of the limit errors in the 
calculated quantities is also computed. The errors 
in the measured quantities are entered with the 
input data. 

Tables 4 to 15 list the phantom generated 
normalized dose and tissue maximum ratio curves 
as a function of depth along with their associated 
limit errors for a variety of conditions of field 
size, SPD and presence of accessories for the 15 MV 
and 6 MV accelerators. The phantom generated 


i t 1 as.a function of 
normalized dose at d,, increases as 
* The suffix, '.For', indicates an executable Fortran 
program. The suffix, '.Dat",, indicatescandara 
file storing constant values. 
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Table 4 


.15 MY 


SPD = i18@cm 


No Accessories 


Field Size = 3.0cm X 3.8cm 
Photon Backscatter N.D. 


Depth 
(cm) 
8.8008 
@.16 
@.32 
8.48 
@.64 
@.95 
ag 
1.94 
2.54 
2.86 
3.18 
3s50 
3.61 


15 MY 


ND 


@.001 
@.168 
@.298 
8.388 
@.456 
8.577 
@.663 
8.759 
@.802 
@.8i1 
8.813 
@.813 
O.et2 


SPD = 


No Accessories 


Field Size = 4.@cm X 4.8cm 


Error 


8.808 
6.888 
8.088 
8.888 
8.808 
- 807 
. 807 
. 886 
. 805 
. 885 
.- 885 
. 885 
. 885 


Qaeaao0n0o0ng 


188cm 


Photon Backscatter N.D. 


Depth 
(com) 
2.290 
@.16 
@.32 
@.48 
9.64 
@.95 
1227 
1.91 
2.54 
2.86 
3.18 
3.52 
3.81 


ND 


@.001 
9.173 
8.303 
@.482 
@.486 
@.614 
@.786 
@.818 
@.862 
8.872 
9.881 
@.882 
8.880 


Error 


Q.800 
8.989 
@.889 
9.088 
8.@98 
8.908 
@.007 
@.006 
8.906 
8.985 
9.085 
@.8@95 
8.9985 


in ND 


in ND 


@.8@01 


TMR 


Pprererreaagsdqanqoqonda 


@.001 


TMR 


HMB#rerrergaagaaqnn9qaondaened 


Error 


Error 


8.82 
0.91 
8.01 
@.01 
8.01 
@.01 
@.01 
8.01 
@.01 
@.01 
@.01 
8.01 
@.01 


- 8B 
-@1 
- 81 


gOagoaoaoaoada 
@ 
ren 


in TMR 


in TMR 
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Table S 


15S MV SPD = 188cm 
No Accessories 
Field Size = 5.8cm X §5.@cm 


Photon Backscatter N.D. = 8.002 
Depth ND Error in ND TMR 
(cm) 

@.88 @.882 @.001 8.80 
@.16 @..173 @.819 Wie Bem 
@:. 32 @.384 @.018 @.33 
8.48 8.407 @.91¢ @.44 
@.64 8.584 @.669 g.55 
OB. 95 @.633 @.889 8.69 
Ten @.728 8.808 o.79 
2.91 @.848 6.887 @.91 
2.54 @.897 @.006 8.97 
2.86 8.9190 @.886 6.99 
@.i18 02.917 @.006 1.982 
3.58 @.922 8.806 1.82 
ait 8.928 @.885 1.88 


rs MY SPD = 18@6cm 
No Accessories 
Field Size = 8.@cm X 8.@cm 


Photon Backscatter N.D. = @.004 
Depth ND Error in ND TMR 
(cm) 

8.80 @.884 @.061 6.88 
@.16 0.194 @.812 8.20 
@.32 @.332 @.812 @.34 
@.48 @.439 @.011 6.46 
@.64 @.532 @.@11 O.5> 
O.95 6.669 @.8018 @.69 
a yar d @.765 @.0198 8.79 
ee Th @.881 8.088 @.91 
2.54 0.948 @.807 6.98 
2.86 @.953 @.807 @.99 
3.18 @.957 @.007 @.99 
3.508 6.964 @.807 1.08 


3.81 @.963 @.086 1.08 


Error in TMR 


@.80 
@.@1 
@.@1 
6.01 
@.01 
@.@01 
@.@1 
@.@i 
@.@i 
6.81 
@.61 
@.81 
@.01 


Error in TMR 


@.8@ 
@.61 
8.81 
@.8e2 
@.82 
8.82 
8.82 
@.01 
6.61 
@.01 
@.01 
@.61 
8.81 
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Table 6 


15 MY SPD = 18@cm 
No Accessories 


Field Size =1@.@cm X1G.@cm 


Photon Backscatter N.D. = @.@86 

Depth ND Error in ND TMR Error in TMR 
(cm) 

8.00 @.006 @.881 @.@01 8.88 
@.16 @.191 8.013 @.19 @.@1 
G.. 32 8.333 @.813 8.34 @.@2 
8.48 @.443 @.812 @.45 8.82 
8.64 8.548 8.012 o..55 @.82 
a. 95 8.677 @.011 8.69 8.82 
1.27 8.774 @.018 Oo. 73 @.82 
Looe 8.897 8.089 @. 92 @.82 
2.54 @.955 @.988 8.97 @.01 
Zang 8.978 8.887 O..99 8.01 
3.18 o. 3G¢ 8.007 8.39 2.01 
3.58 8.983 @.807 1.98 @.01 
3,82 8.983 8.207 1.80 @.01 


1S my SPB = 1@8cm 
No Accessories 
Field Size =12.@cm xX12.@8cm 


Photon Backscatter N.D. = 8.887 

Depth ND Error in ND TMR Error <in. THR 
(cm) 

8.820 Q.007 8.801 @.@01 8.82 
8.16 0.205 @.015 By 22 @.82 
9.32 8.349 @.814 o.35 9.02 
8.48 @.462 @.014 8.47 9.82 
8.64 @.556 @.013 @.56 8.82 
o.95 @.699 @.012 8.78 @.82 
1.27 0.796 8.012 8.89 8.02 
2.98 9.914 @.019 8.92 8.982 
2.54 8.973 8.089 8.98 8.62 
2.86 8.987 @.888 @.99 @.82. 
3.18 0.9393 9.808 1.09 @.02 
3.58 0.992 @.088 1.00 8.82 
3.81 8.993 8.088 1.08 8.82 
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Table 7 


1S BY 


SPD = 1@8cm 


No Accessories 


Field Size 


Depth 
(cm) 
@.80 
@.16 
9.32 
@.48 
@.64 
8.95 
Wena d 
tee Je 
2.54 
2.86 
3.18 
3.59 
pao 


15 MY 


ND 


@.809 
8.207 
@.354 
@.463 
8.568 
8.704 
6.806 
@.921 
8.975 
6.991 
8.997 
1.882 
1.884 


SPD: = 


No Accessories 


Field Size 


Error 


=15.0cm X15.8cm 
Photon Backscatter N.D. 


@.001 
8.817 
@.016 
@.016 


Qeaoaoaa0n go0o0 a 
ey) 
as 
® 


Q 
g 
w 


180cm 


=28.0cm X20.8cm 


Photon Backscatter N.D. 


Depth 
(cm) 
8.88 
@.16 
8.32 
@.48 
0.64 
oO. 95 
py -ae 
pa de 
2.54 
2.86 
2.10 
3.58 
3.81 


ND 


@.812 
8.205 
9.353 
0.462 
@.565 
@.7808 
8.882 
@.922 
@.982 
0.997 
1.661 
1.8084 
1.807 


Error 


8.882 
@.828 
@.826 
8.019 
9.818 
0.817 
8.816 
@.614 
6.812 
@.011 
@.@11 
@.011 
@.811 


in ND 


in ND 


8.889 


TMR 


@.@1 
9.21 
8.35 
@.46 
@.56 
8.78 
@.8@ 
6.92 
8.97 
G:.93 
as99 
1.08 
1.68 


8.8012 


TMR 


@.@1 
8.28 
8.35 
@.46 
8.56 
a.7@ 
8.88 
0.92 
8.98 
oO. 99 
O.99 
1.80 
1.68 


Error in TMR 


@.68 
8.82 
@.82 
@.62 
@.82 
@.62 
@.82 
8.82 
@.82 
6.62 
@.82 
@.@2 
@.62 


Error in TMR 


8.88 
8.02 
8.82 
8.82 
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Table 8 


1S mv SPD = 108cm _ 
No Accessories 
Field Size =25.@cm X25.@8cm 


Photon Backscatter N.D. = 8.015 

Depth ND Error in ND TMR Error sn thnk 
(cm) 

@.88 @.015 8.802 @.@1 @.08 
@.16 8.268 0.824 @.2i @.83 
@.32 @.357 @.023 Q@.36 8.93 
8.48 @.471 @.022 8.47 8.863 
@.64 @.565 @.022 @.56 @.93 
8.395 8.710 8.828 @.71 8.83 
Poer 8.802 9.019 @.88 @.93 
1.94 @.919 @.016 @.92 @.93 
2.54 8.977 @.9014 8.97 8.83 
2.86 8.994 @.014 o.99 8.03 
3.18 @.999 8.813 @.99 @.983 
3.58 1.001 8.013 1.08 @.83 


3.81 1.884 @.8013 1.08 @.83 
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Table 9 


15 MV SPD = 75cm 
No Accessories 
Field Size =1@.8cm XK1@.8cm 


Photon Backscatter N.D. = 9.007 

Depth ND Error in ND TMR Error 210° iMk 
(cm) 

8.80 0.007 @.0@1 @.@1 8.88 
@.16 @.282 @.014 @.21 @.82 
@.32 8.339 @.013 O. so 8.802 
@.48 @.449 @.813 8.46 @.82 
@.64 Q@.539 @.8i2 8.55 @.62 
0.95 8.682 @.01ii1 8.69 8.82 
1.27 0.777 @.@1i1 8.79 @.82 
2 CA 9 @.896 8.889 @.91 6.62 
2.54 @..395S 8.088 8.98 @.82 
2.86 @.971 8.088 @.99 8.82 
3.18 @.978 8.808 1.08 @.@2 
3.58 @.981 6.807 1.00 8.82 
Sant @.982 @.8007 1.00 @.@1 
15 MY SPD = 7Scm 

@.64cm Lucite Accessory 

Field Size =1@.@cm X1@.@8cm 

Photon Backscatter N.D. = 8.807 

Depth ND Error in ND TMR Error in TMR 
(cm) 

@.80 @.007 @.a01 @.81 8.88 
@.16 @.2802 @.019 OQ. er @.82 
a Pat @.342 6.018 @.35 @.82 
@.48 @.449 @.018 8.47 8.62 
@.64 8.538 6.017 @.56 8.82 
@.95 8.679 @.016 8.78 8.02 
SAr-ate @.765 6.815 G.79 8.82 
lade Ie a.88sa 8.813 O.g1 @.82 
2.54 @.942 @.011 a.3% 8.82 
2.86 @.9S5 @.611 @.99 @.82 
3.18 @.961 @.@11 6.59 @.82 
3.58 @.963 @.818 1.08 @.82 
3.81 @.966 @.012 1.08 8.82 
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Table 18 


15 My SPD = 75cm 
No Accessories 
Field Size =30.8cm X38.@8cm 


Photon Backscatter N.D. = 9.8018 
Depth ND Error in ND TMR 
(cm) 

8.88 @.018 8.882 @.02 
@.16 @.21@ @.829 @.21 
Grse Or3asa @.628 @.35 
@.48 8.462 8.027 @.46 
@.64 @.551 @.026 @.55 
@.95 @.692 @.024 @.69 
Lee? 0.782 8.823 8.78 
1492 @.908¢8 8.920 @.98 
eeo"g @.967 @.017 @.97 
2.86 0.984 @.616 8.99 
hs Gs9S3i @.016 Grss 
3708 @e395 8.6816 1.02 
3-81 @.998 @.8015 1.80 


15 MV SPD = 75cm 
@.64cm Lucite Accessory 
Field Size =30.8cm X3@.0cm 


Photon Backscatter N.D. = 8.818 
Depth ND Error in ND TMR 
(cm) 

@.00 @.018 @.802 @.02 
@.16 @.18@8 @.845 0.19 
8.32 8.382 @.843 G6 .uce 
8.48 @.398 @.842 @.41 
@.64 @.463 @.041 @.49 
@.95 @.599 6.837 8.63 
eaeG @.688 @.@35 a RA dad 
pA @.816 @.831 @.86 
2.94 @.901 8.827 6.95 
2.86 0.927 @.826 @.98 
ee AS @.936 @.826 6.99 
3.5U @.946 6.825 1.00 


3.81 8.948 @.824 ce 


Error in TMR 


8.80 
8.83 
8.83 
8.83 
8.03 
8.83 
6.83 
@.03 
8.83 
@.@3 
8.63 
2.03 
8.83 


Error in TMR 


8.88 
8.85 
8.85 
8.05 
8.06 
6.86 
@.06 
@.85 
8.85 
8.85 
8.85 
@.85 
8.85 
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Table ii 


15 MY SPD = 188cm 
3.16cm Lucite Accessory 
Field Size =30.8cm X3@.@cm 


Photon Backscatter N.D. = 8.8013 
Depth ND Error in ND TMR 
(cm) 

@.80 @.813 8.082 @.61 
8.16 Dd? @.828 O.17 
8.32 @.319 0.027 Q@.31 
8.48 @.433 @.026 @.42 
@.64 8.537 8.825 Osce 
@.95 @.686 @.823 @.67 
Lace @.782 8.822 8.76 
1.91 @.914 @.819 @.89 
2.54 @.987 @.017 8.96 
2.86 1.804 @.016 6.98 
2. 18 1.014 @.016 8.99 
3.50 2arS @.015 1.08 
c= § 1.024 @.015 1.08 


15 MY SPD = 14@8cm 


No Accessories 
Field Size =38.@cm X3@.0cm 


Photon Backscatter N.D. = 80.017 
Depth ND Error in ND TMR 
(cm) 

@.80 @.017 8.902 @.82 
@.16 @.194 6.835 @.28 
9.32 @.332 6.834 8.34 
8.48 @.445 8.833 8.45 
0.64 @.540 8.032 @.55 
6.35 @.678 8.829 @.69 
AN g @.749 8.828 @.76 
12.91 @.876 8.824 8.89 
2.54 8.949 @.82i @.96 
2.86 @.965 8.828 8.98 
2.16 @.969 8.828 6.98 
3.358 @.976 @.819 o739 


3.62 @.986 @.019 1.08 


Error in TMR 


@.88 
@.83 
8.83 
@.03 
8.983 
8.03 
@.83 
@.83 
@.83 
@.83 
6.83 
8.83 
@.83 


Error in TMR 


8.80 
8.04 
6.04 
8.04 
6.84 
8.84 
8.04 
8.04 
@.04 
8.04 
8.04 
@.84 
8.84 
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Table i2 


6 MV SPD = 100cm 
No Accessories 


Field Size = 5.@cm X 5.@cm 


Photon Backscatter N.D. = 
Depth ND Error 
(cm) 

@.80 Q@.802 @.601 
@.16 @:391 @.011 
0.32 @.581 @.818 
@.48 8.698 @.3989 
@.64 8.785 @.2@@8 
@.95 0.878 @.088 
T.i22 @.917 @.007 
Last @.930 @.007 


6 MY SPD = 108cm 
No Accessories 


8.062 


in ND TMR 


FP aggagaaa0ndd 
aN | 
Ui 


Field Size =18@.8cm x¥10.8cm 


Photon Backscatter 


N.D. 


Depth ND Error 
(cm) 

@.00 @.009 8.001 
@.16 @.428 @.814 
6.32 6.617 8.012 
@.48 o.739 @.011 
8.64 @.831 @.8128 
8.95 @.926 @.009 
avert @.968 8.869 
°° @.978 8.809 


8.089 


in ND TMR 


Error 


Error 


8.80 
8.82 
Q@.82 
9.02 
@.8e 
8.82 
0.82 
@.82 


@.88 
@.82 
@.82 
@.82 
@.82 
@.82 
@.82 
@.82 


in TMR 


in TMR 
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Table 13 


6 MV SPD = 16@6cm 


No Accessories 


Field Size 


Depth 
(cm) 
8.80 
@.16 
@.32 
@.48 
@.64 
@.95 
a b= 5a 
2%.91 


6 MY 


ND 


@.918 
8.437 
0.648 
@.766 
9.859 
@.954 
9.994 
1.885 


SFPD §= -160cm 


No Accessories 


Field Size 


Error 


=20.8cm X28.0cm 
Photon Backscatter N.D. 


@.802 
@.619 
@.817 
@.016 
@.814 
@.013 
@.8012 
@.8i2 


=3@.8cm X36.8cm 


Photon Backscatter N.D. 


Depth 
(cm) 
8.88 
Q@.16 
8.32 
8.48 
9.64 
@.95 
Lea? 
ie S1 


ND 


@.826 
@.445 
@.647 
8.772 
8.865 
@.957 
@.996 
1.885 


Error 


8.983 
9.824 
@.822 
8.828 
0.018 
@.816 
@.015 
@.015 


in ND 


in ND 


8.818 


TMR 


@.82 
@.43 
6.64 
8.76 
8.85 
8.95 
8.99 
1.06 


8.626 


TMR 


Fogo gqagagaga 
a | 
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Error in TMR 


- 88 
- 82 
82 
02 
82 
Be 
82 
82 
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Error in TMR 


8.82 
@.83 
8.83 
@.83 
8.83 
8.83 
@.03 
8.83 
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Table i4 


6 MV SPD = 7Scm 


No accessories 


Field Size 


Depth 
(cm) 
8.82 
@.16 
8.32 
8.48 
@.64 
@.95 
a ead 
Ure Jak 


6 MV. 4SPD»= Sem 


ND 


8.029 
8.468 
@.667 
8.789 
2.878 
@.966 
1.903 
ac OL2 


Error 


=30.8cm X38.8cm 
Photon Backscatter N.D. 


@.083 
8.826 
8.823 
@.821 
@.019 
@.017 
@.016 
@.616 


@.64cm Lucite Accessory 


Field Size 


=3@.8cm X30.6cm 


Photon Backscatter N.D. 


Depth 
(cm) 
@.80 
0.16 
@.32 
8.48 
8.64 
8.95 
eh 
i Ot 


ND 


@.829 
@.481 
8.658 
2.763 
@.852 
0.937 
8.973 
@.981 


Error 


@.883 
8.838 
@.834 
0.032 
@.029 
@.026 
@.825 
8.824 


in ND 


in ND 


@.029 


TMR 
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. 83 
.49 
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. 88 
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Error 
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Table i5 


6 MV SPD = 75cm 
3.18cm Lucite Accessory 
Field Size =3@.@cm X3@.8cm 


Photon Backscatter N.D. = 8.829 

Depth ND Error in ND TMR ERO mal ieee 
(cm) 

@.08 @.829 8.803 @.a3 @.aa 
@.16 @.497 8.848 8.49 @.85 
O.o2 8.678 @.8036 0.67 6.85 
g.48 @.784 8.833 8.78 8.85 
@.64 8.877 8.838 @.87 @.85 
@.95 @.963 @.827 8.96 8.85 
tee 2.93399 8.826 @.99 8.85 
4.92 1.886 8.825 1.00 @.85 


6 MY SPD = 108cm 
3.18cm Lucite Accessory 
Field Size =10.8cm X10.0cm 


Photon Backscatter N.D. = 8.889 

Depth ND Error in ND TMR Error in TMR 
(cm) 

@.80 8.003 8.801 @.@1 8.80 
@.16 @.420 @.014 @.43 8.02 
@.32 8.621 @.813 8.64 8.82 
@.48 8.742 @.012 @.76 @.82 
@.64 @.832 @.@11 @.85 8.82 
@.95 9.925 @.818 8.95 8.82 
Pisce @.965 8.809 8.99 8.82 
Laos @.976 @.a09 1.08 @.@2 
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field size. The ratic of the phantom generated 
normalized dose at oe to the phantom generated 
normalized dose at a field size of 10 cm x 10 cm 
are the phantom generated output factors. Figures 
65 and 66 illustrate the phantom generated output 
factors compared to the total output factors for 
6 MV and i5 MV, respectively, as a function of 
field width. The phantom generated output factor 
increases as a function of field width due only to 
increased scatter to ae as a function of field 
width. The increase in output as a function of 
field width due to contamination has been removed. 
There is only a few percent difference between the 
total and phantom generated output factors because 
there is little contamination present at ee 
However, the phantom generated output curve is 
greater than the total output curve at small field 
sizes and reaches a plateau at a field size yreater 
than 20 cm x 20 cm. 

Tables 4 to 8, 12 and 15 indicate that the phantom 
generated TMR curve is independent of field size. 
Tables 9 to 11, 14 and 15 illustrate that the phantom 


generated TMR curve is independent of SPD and the 


presence of accessories. 
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Figure 65. The phantom generated output factor curve compared to the 
total output factor curve at 6 MV. 
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Figure 66. 
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The phantom generated output factor curve compared to the 
total output factor curve at 15 MV. 
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Zale Blocking the Primary Beam 


A conical ingot of cerrobend or Lippowitz's metal (an alloy of 
lead, tin and bismuth which will attenuate half the dose in a thickness 
of 1.05 cm at 6 MV and 1.7 cm at 15 MV) 10 cm in length was produced 
such that the sides of the ingot were parallel to the beam direction 
when it was placed on the accessory tray. The ingot blocked the primary 
beam from a conical region of space underneath it. At 100 cm SPD the 
maximum diameter of the blocked region was 4.3 cm. 

If the surface dose with the block in and out is constant, then 
the electron contamination and photon backscatter are not perturbed by 
the addition of the beam block. Table 16 lists the surface normalized 
dose at 100 cm SPD both with and without the beam block in place and 
the difference between these conditions under a variety of circumstances. 
There is a difference between the surface dose with the block in and 
out so the block does modify the surface dose. 

A field size at 100 cm approximately 4 cm x 4 cm* is being blocked 
by the ingot. The surface normalized dose for a beam of 4 cm x 4 cam 
with an accessory tray in place is listed in Table 16. The surface 
dose for a 4 cm x 4 cm field size is almost the same as the difference 
between the dose with the beam block in and the beam block out. The 
remainder is the contribution to the surface dose due to photon back- 
scatter between 4 cm x 4 cm to the field size listed. Therefore, the 
beam block stops an amount of electron contamination corresponding to the 
projected area of the block at SPD = 100 cm as well as blocking the 
primary beam. 


The normalized dose as a function of depth in a phantom was 


* The area of a 4.3 diameter circle is approximately the same as the 
area Of a square with sides measuring 4 cm. 
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Table 16 
Surface Normalized Dose At 1@8cm SPD 
With And Without A Cerrobend Primary 
Beam Block Present On The Accessory Tray 
Block Block 
Energy Field Size Out In Difference 
6 MY 18cm X 18cm @.169 8.071 @.a98 t 9.2007 
6 MY 28cm X 28cm 8.380 @.2084 @.096 = @.a07 
1Syny 28cm X 28cm OB. 271 8.228 @.951 7 98.904 
Surface Normalized Dose At 10@cm SPD 
With The Accessory Tray In Place For 
A dom X 4com Field $1ze 
Energy Surface Normalized Dose 
6 MV @.a@92 ¢ 38.0@5 


15 MY 9.948 * @.993 
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obtained with the block in and the block out. The difference between 
these readings is due to both a lack of the primary beam and a partial 
loss of contamination. 

The loss of electron contamination is reflected by the ioss of 
surface dose. If it is assumed that the scattered photon component of 
the contamination is attenuated by the beam block to the same extent 
as the electron component, then the contamination TMR curve obtained 
in Section 2.11.1 can be used to obtain the phantom generated N.D. 
and TMR for the portion of the beam blocked. This was obtained using 
the program Seprimary. For described in Section 2.10262 The:result as 
shewn in Table 17. The phantom ‘generated TMR values agree with the 


phantom generated TMR values obtained in Section 2.11.2. 
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Table 17 


15 MY 


Field Size 


Difference Between The Primary Beam 


Blocked And Not Blocked With The 


Contamination Dose Removed 


SPD = 


1@8cm 
@.@64 Lucite Accessory Tray Left In When Block Removed 
=28.8cm X20.8cm 


Photon Backscatter ND’. 


Depth 
(cm) 
8.08 
@.16 
yee Pa 
8.48 
0.64 
im be 
L.27 
2 im 
2.54 
2.86 
3.18 
a OU 
2.61 


ND 


8.003 
8.288 
Ou356 
@.447 
0.543 
GO. 675 
8.765 
8.874 
0.923 
0.935 
@.939 
0.941 
Ga339 


Error 


8.083 
8.017 
@.6016 
6.016 


‘86.815 


@.814 
8.013 
@.611 
6.809 
8.889 
8.889 
8.088 
6.888 


in ND 


8.883 


TMR 


1.82 


Error 
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2.13 Penetration Curve Outside a Magnetically Swept Field 


The techniques described in Section 2.7 (Penetration Curves Outside 
the Field), Section 2.10 (Sweeping Electrons from the Field using a 
Permanent Magnet) and Section 2.11 (A Method for Separating Contamination 
from Phantom Generated Penetration Curves) can be used to determine 
the penetration curve of the scattered photon contamination component. 

An ion chamber in a partial phantom as iliustrated in Figure 30 
was placed 3 cm outside an open 10 cm x 10 cm field. Readings as a 
function of the thickness of overlying phantom were obtained with the 
Magnet in place. The probe was placed on the side of the field opposite 
to the direction of deflection of the electrons. The only contribution 
to the dose under these conditions is due to scattered photons and 
electrons produced in air between the magnet and the partial phantom. 
The readings as a function of overlying phantom normalized to the 
maximum reading obtained is shown in Figure 67. This curve is 
characteristic of photons. If it is assumed that the surface dose is 
entirely due to electrons produced in air (in other words, neglecting 
the backscatter surface dose produced by the scattered photons), then 
this contribution can be subtracted using the technique of Section 2.11. 
The result is also shown in Figure 67. The scattered photon build- 
up curve is characteristic of a megavoltage photon beam between 6 MV 
and 15 MV. The depth of the maximum dose is at approximately 2.5 cm 


compared to 1.6 em for 6 MV and 3.5 cm for 15 MY. 
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Figure 67. 


Oo <— 
® 
| zs) 
_ 
1S) 
3 
Ss 6 
=e¢é 
ec 
= ~ 
s 2 
edzas 
25 
cS 
= £€ 
= a6 
CGC = 
= C¢ - 
6 0 E 
Ores pes 
5 5 
gS 
33 ir 
he O 
oo 
aa 
ts ot 
6 5 
e 8 _— 


Qisioy 4) Ih 30. 0 - S, oe a 


a 4JSOd 3AUV134 


The scattered photon contamination build-up curve obtained 
3 cm outside a magnetically swept 10 cm x 10 cm field at 
100 cm SPD for the 15. MV beam. 
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3. DISCUSSION 


It is truth very certain that, 
when it is not in our power 
to determine what is true 

we ought to know, 

what is probable. 


Rene Descartes 
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3.1 Independent Modelling of the Phantom Generated Build-Up Curves 


Sig li 


SediarZ 


introduction 

The phantom generated build-up curves were determined 
using a computer model. This was done to verify the phantom 
generated build-up curves produced by the subtraction of 


the contamination contribution from the total build-up curves. 


The computer program called Buildup3.For was written in Fortran 


-77 to run on the VMS operating system of a VAX 11/760 
computer. Buildup3 is listed in Appendix 6. Buildup3 
models the generation of charged particies set in motion in 
water from megavoltage photon spectra and has the capability 
to harden the spectra. The charged particles are followed 
through the phantom using the Monte Carlo method and the 
dose absorbed by the phantom is determined. 

Data Preparation Programs 

Two additional programs prepared the look-up Tabies 
for the program, Buildup3. The program, Quadint3.For, 
interpolates the total scattering and stopping power data. 
The program, Normdist.For, generates the value of the error 
function, erf(x). 

The amount of energy lost per unit path iength, stopping 
power, and the root mean square angular scattering power for 
representative energies are stored in a data file called 
Master.Dat. Quadint3.For interpolates using the quadratic 


method (47), the scattering and stopping power data from 


Master.Dat. The output of Quadint3.For is a listing of the 


stopping and scattering power data for 0.1 MeV increments 
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in energy from 0.2 MeV to 18.0 MeV. This listing is stored 
in a file called Quadint3. Dat. Buildup3 reads the value of 
Quadint3.Dat before execution. 

Normdist.For generates the error function, erf (x): 

s 2 
erf(x) = al ent aae Gr 

-X 
Normdist.For is listed in Appendix 7. The value of 
erf(x) represents the probability of a value of t being 
within the limits of integration, [-x,x], if t follows a 
Gaussian distribution with a mean equal to 0 and a standard 
deviation equal to l. 

The file named Gaus.Dat lists the inverse error func- 
tion from 0.00 to 0.99 in increments of 0.01. It was 
generated by interpolation of the inverse error func- 
tion from the above values of the error function. - Gaus. - 
Dat is used to determine a normally distributed random scat- 
tering angle from a uniformly distributed random number 
between 0 and 1. Figure 68 graphically illustrates the 
procedure used in a number of places in the program and 


will be discussed further in the following section. 


Generation of Charged Particles Set In Motion by Photons 


Buildup3 requires the photon number spectrum to be 
entered as a function of photon energy. The number spectrum 
was not measured but was interpolated from spectra pub- — 


lished in a paper by C.W. Sandifer and M. Taherzadeh of 
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Figure 68. 
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The probability of an angle, 9, occurring when the root 
mean scattering angle, o, is equal to R. Conversely, ita 
random number, R, between 0 and 1 is chosen, a ratio, 8/o, 
can be found. Knowing o, then 9 is determined. This figure 
illustrates that a random number of 0.71 results in a value 


Of0/o-.0f; 1.00. 
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E.G. and G. Inc. (48). They measured thick-target 
bremsstrahlung spectra produced by electrons impinging 

on a 0.2 radiation-length gold-tungsten (Z of gold = 79, 
Z of tungsten = 74) target for electron energies in the 
range of 5.3 to 20.9 MeV. The electrons were accelerated 
by a 90° bent beam linac. An aluminum block was placed 
behind the target to absorb low energy electrons emerging 
from the target. 

The Siemens Mevatron -20 and Mevatron -6 are similar 
except for a number of features. The target is pure gold 
in both the Mevatron -20 and Mevatron -6. The Siemens 
accelerators have a 270° bend instead of the research 
accererator' s 90° bend. 

Sandifer and Taherzadeh measured spectra produced by 
electron energies very close to the nominal photon energies 
of the Siemens accelerators; 5.3 MV compared to the Siemens 
6 MV and 14.8 MV compared to the Siemens 15 MV. The desired 
Siemens nominal energy was determined by normalizing the 
nearest measured differential energy spectra greater than 
and less than the desired energy to unity and linearly in- 
terpolating between these spectra to obtain the spectrum 
at the desired energy. Table 18 summarizes the determina- 
tion of the 6 MV differential spectrum. The number spectrun, 
which is the number of photons in discrete bins as a 
function of energy, is determined by multiplying the value 


of the differential spectrum of an energy bin by the energy 
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Table 18 


Maximum 
Bin 
Energy 
(MeV) 


8.64 


Determination of the 6 MU Differential 


Spectrum from Measured Spectra 


Experimental 


5.3 my 18.8 MV 
dgN/dhvu agN/dhy 
Ci/MeYyV) C1i/MeYV) 
-1 -1 
5.34 * 18 2.53 X 18 
1.88 1.80 
-1 -1 
7.88 X 18 6.76 X 12 
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-1 -i 
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6.0 MV 
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CivMeV) 
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width of the bin: 


N, = (S-). (shy), (3.1.2) 
ak 


where (Ahv) ; is the width of the bin 


(-) is the value of the differential energy 
i spectrum at the bin energy 


N. is the number of photons in the energy bin 


The number spectrum is smoothed graphically (see 


Figures 69 and 70) and is normalized so that the sum of 


the bins. issunity: 


a (391.3) 


\ 
where N. is the smoothed value of number spectrum at the 
bin energy 


F. is the smoothed normalized number spectrum at the 
bin energy 


The mean eriergy of the photon number spectrum is 
determined by a linear weighting with respect to energy: 


hv = = =F. (hv). al 
Dy Maen pans (See) 


because LF. = 1 
at 
where hv is the mean value of the photon number spectrum 


The mean value of the 6 MV photon beam is 1.5 MeV and the 
mean value of the 15 MV photon beam is 3.9 MeV. 
The central energy of the photon bin, the value of the 


number spectrum, F.; the total attenuation coefficient, 
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The smoothed and unsmoothed unnormalized primary photon 


Figure 69. 
number spectrum at 6 MV. 
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40 


fine Yariab 2 16 
PHOTON ENERGY (MV) 


Figure 70. The smoothed and unsmoothed unnormalized primary photon 
number spectrum at 15 MV. 


Te 


u,.,/p for water, and the total Compton absormtion coef- 


tot 


ficient, o,_,/p for water, are listed in files named 


tot 
Photon6 for 6 MV and PhotonlS for 15 MV. 


The subroutine, Impulse, generates the charged parti- 
cle set in motion spectrum due to photon interactions in a 
1 mm slab of water. The photons are assumed to be 
initially directed normally into the slab. 

The maximum photoelectric attenuation coefficient for 
the spectrum occurs at the lowest energy bin (0.13 Mev) and 
is the only 2% of the value of the Compton attenuation 
coefficient. The ratio between the photoelectric and 
Compton attenuation coefficients for the center of the 
second energy bin (0.21 MeV) is only 0.5% and the ratio 
for the third energy bin (0.31 MeV) is less than 0.1% (9). 
Therefore, the photoelectric effect need not be and was 
not taken into account. 

A flow chart of the subroutine, Impulse, is shown in 
Figure 71. The loop at the left of the figure is repeated 
for both the Compton and, if required, the pair production 
interactions. 

The number of photons from each photon energy bin in 
a 1 m slab, Nis (the variable used in the program for 
N,” Hsi nterace(1),..that, 1S, Ni” + nteract(i)*), is deter- 
mined knowing the number fluence incident on the slab, F., 
and the total mass attenuation coefficient, (u,/e); 


((u,/0)5 = mu(i)*): 


* The arrow indicates a transformation from a mathematical 
variable to a variable or array element defined in a 
program if the variables are different. 
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PAIR 
PRODUCTION 
TAKEN 
INTO ACCOUNT? 


INITIALIZE TO ZERO 
THE ENERGY ANGLE 
DISTRIBUTION ARRAY 


ASSIGN NUMBER OF ELECTRONS 
TO ENERGY ANGLE 


ALL 
PHOTON 
BINS? 


DISTRIBUTION ARRAY Yes 


CALCULATE ELECTRON 
RECOIL ANGLE FOR 
EACH AVAILABLE 
ENERGY BIN 


No 


CALCULATE NUMBER 
OF PHOTONS INTER- 
ACTING IN SLAB FOR 
EACH PHOTON BIN 


REDUCE NUMBER OF 
PHOTONS IN BINS 


CALCULATE X-SECTION 
FOR EACH AVAILABLE 
ENERGY _BIN 


DIVIDE AVAILABLE 
ENERGY INTO BINS 


CALCULATE NUMBER OF 
COMPTON ELECTRONS 
AND PAIRS FORMED 


CALCULATE AVAILABLE 
ELECTRON ENERGY 


Figure 71. A flow chart for the subroutine, Impulse. The loop at the 
left is repeated for both Compton and pair production 


interactions. 
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If the Compton scattered photons do not interact again 
(corresponding to narrow beam geometry), then the only 
photon bin which needs to be modified is the bin from 
which the primary photon came: 

/o) (0.1 g/cm) 


- 1G), 
OT atten Oy betore pre Geel: 


slab slab (3.1.6) 


The number of the interacting photons at a given 
energy bin due to the Compton interaction , NCOMP.. , 1S 
given by: 


GE Ne (30137) 
NCOMP i = (OS OE i e ° 


The number due to pair production, NPP. , is given by: 


_ ez 1.8 
NPP, = 2(N,” - NCOMP;) (3.1.8) 


The above expression assumes that positrons can be treated 


like electrons. Additionally, the total attenuation coef- 
ficient takes into account pair production in the field of 
atomic electrons, "triplet production'', however, the factor 
7 does not take this completely into account. MAt =the 


mean photon energy of the 15 MV beam (14 MeV), the 
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triplet cross-section is about 3% of the pair production 
cross-section and the pair production cross-section is, 
in turn, only about 6% of the total attenuation cross- 
section (9). This results in slightly fewer charged 
particles with a slightly larger mean energy. The 
failure to take into account triplet production rigorously 
wili not significantly affect the distribution of charged 
particles set in motion. 

The kinetic energy available to a charged particle 
is always less than the photon energy, hv(hv > xnrg(1)). 


The maximum energy of a Compton electron, see 


> 


Shere + nrgmax), can be derived from the conservation of 
energy and momentum: 


ae 
iD 


taaen Ze +1 hv er) 


The electron may have any energy up to and including 
the energy, peed The difference between the photon 
energy and the kinetic energy of the recoi electron; T, 
is the energy of the scattered photon. 


The available kinetic energy of the pair production 


interaction, Th (Ty + nrgmax), must be shared between the 


electron and positron: 


Ty = hy, ~. 2.022 Mev. (56: leek OH) 


where 1.022 MeV represents the combined rest mass of the 


electron and positron. 


The available kinetic energy is divided evenly into 
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20 bins. The kinetic energy dependent part of the 
differential cross-section with respect to kinetic energy 
is calculated at each of the 20 kinetic energy bins. The 
kinetic energy differential cross-section for the commecn 
effect do/dT is given by the Kiein-Nishina equation (see 
Equations 1.:2..5, "andy 132= 6)". 

The absolute differential cross-section is not 
calculated because it is not needed. The number of photons 
interacting by the Compton effect is known (see Eqn. 3.1.7). 
The relative number interacting at each available kinetic 
energy bin is determined by normalizing each cross- 
section factor, C(hy,T,), to the sun, rC(hv,T.), where 
the sum is carried out over all 20 octane kinetic 
energy bins. 

The differentiai cross-section for the pair produc- 
tion interaction , dr/dT, is given by Equations 1.2.7, 
1.2.8 and 1.2.9. When considering the. Compton interaction 
calculation, the kinetic energy dependent part of the 
differential cross-section, P(hv,T), needs to be calculated. 
The sum, uP (hy,T.), is used to normalize the cross-section 


‘s - 
so that the appropriate number of electrons 1s set in 


motion at each energy, T;. 

jareeeheeeen zenith direction of the Compton recoil 
electron, with respect to the primary photon direction, 
8, is given by Equation 1.2.4. If Equations 1.2.3 and 


1.2.4 are combined with Equation 3.1.9, the Compton 
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electron, with respect to the primary photon direction, 
8, 1S given by Equation 1.2.4. If Equations 1.2.3 and 
1.2.4 are combined with Equation 3.1.9, the Compton 


electron recoil angle, 6, can be determined: 


ES 
Spi GalorD 


6 = arctan = 
) VT 


An element of the energy-angle distribution array 
is incremented by an amount corresponding to the number 
of electrons set in motion at that kinetic energy and 
zenith angle: 


C(hyv,T, ) 


Pp ZC(hv,T, J 
1 


Sf Bas (3.1.12) 


ke ~ kg * NCOMP 


where Ske (Srp + SPECTR(k,2)) is the energy-angle distribu- 


tion array 

k is the energy bin number 

2 is the angle bin number 

The available kinetic energy bins, represented by i 

for the calculation of the number of interactions at each 
available energy, is not the same as the energy-angle 
distribution kinetic energy bin, represented by k, in the 
above expression. There are 90 energy bins for the energy- 
angle distribution array spaced evenly between 0.0 and 18.0 
MeV. Each energy bin of the energy-angle distribution is 
0.2 MeV wide. Therefore, if the primary photon energy is 
low, many available kinetic energy bins may be mapped into 
one energy bin of the energy-angle distribution array. 


There are 36 angle bins, represented by 2, in 
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Equation 3.1.12 to cover any zenith angle between 0 and 


10. 


The 36 bins are distributed uniformly so that each 
bin is 5° wide. It is assumed that the initial zenith 
angle for the bin equals the center of the bin. Therefore, 
the possible initial zenith angles are Cue hase igs ae 
iyeisoa The maximum zenith angle for an electron 
produced by a Compton photon interaction is 90°, however , 
there is a small possibility of a charged particle having 
an angle greater than 90° when produced by the pair 
production effect. 

The root mean square zenith angle, Gyr, between the 
trajectory of a charged particle and that of the primary 
photon is given by Equation 1.2.11. It is assumed that 
the distribution is Gaussian (5). 

The zenith scattering angle is chosen by the direct 
or inversion method of random sampling. The cumulative 
probability, P(t), associated with a random variable, T, 
is: 

P(t) = Probability (T<t) CS. 3) 
If P(t) is normalized, then 0<P(t)<l. By generating a 
random number, R, such that 0<R<1, and identifying R with 
P(t), a random sample, T, may be obtained from chendis- 
tribution, P(t), by inversion (49). The cumulative 
probability for a particle having a zenith angle, ¢ (¢ > g), 
when the root mean square angle is (9°)3, is given by: 


2 
P(t) = erf(t) = =| Bas: (Sale) 
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where t = 6/Ko°)? 
A numerical method was chosen to invert the above 
expression. P(t) was generated for values of t between 
0 and 3.0. The value of t as a function of P(t) is listed 
in the file, Gaus.Dat. A random number, R, is chosen and 
is identified with P(t). A value, t, corresponding to 
P(t) is obtained from a look-up table (see Figure 68). 
After the charged particles have been set in motion 
by photons from all the photon bins, the energy-angle 
distribution is passed back to the main program of 
Buildup3. The angular distribution is shown in Figure 72 
for a 15 MV beam interacting in a slab 0.1 cm thick. The 
mean charged particle energy (in MeV) is indicated as 
the numbers in the bins. The characteristic bilobal 
forward peaked distribution is evident. The forward 
directed charged particles have a much higher energy 
than the charged particles set in motion at larger angles. 
There are very few electrons set in motion at angles 
greater than 60° from the initial photon direction. 
3.1.4 The Monte Carlo Program 
Monte Carlo methods are applied to the solution of 
physical problems when the microscopic behaviour is 
better understood than the macroscopic behaviour or the 
geometry of the problem does not lend itself to a closed 
analytic solution. The process to be studied must have 
some intrinsic random nature. In the case of charged 


particle interactions with matter, the slowing down and 


imewted 4 to sie 3 wt spl or pani om 
hoters ar (7)2 to nok sant capa uote ot” Wie SD ni hl 
ban abeois at ot cocoa well ie eye tit it “ 
of gttitnogearsty 4 ta Ae faI% tate beittaneBi co 
‘64 sweet ope) otvies quid « emt bentstdo at (3)9 
mites at 292 need Svad. eufattihey'begiad edt nl 
ai heat se aoe 
lo merguty tem sf Od Hope bozeaq et moisidizsel 
avceit chi awode ef noi mediraeib tehugls aft oot 
od! Said? om 1.0 dete a me galeseoannt rood ‘vw 2 8 10h 
-¢ betaarhat ai” (¥eM tt) ‘plone sioitaig bogiad> esa 
"radians ated sate oA ta Pao ssh) 
inewo? oft habia ar iitele=is bing inet 
cyseco Tee aga avit eetobaeeg ila 


2eiunn veeisl Fe seks’ at i x 


eoloms $e ooavim me! tse 
sini toasts nnn — a nadld seer 


| ot: er qe: 

att) 46 susorenied cicoomoraae wf 
beavis 8 of Meath BheL son ah 

avait 92th ‘bothie a8 od 209 


NAS, 


100 


0.80 


Figure 72. The angular distribution of charged particles interacting 
in a 0.1 cm slab at 15 MV. The mean charged particle energy 
(in MeV) is indicated by the numbers in the histogram. 


Each angle bin is 5° wide. 
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the scattering direction is microscopically indeterminate. 
The number of Couiomb interactions resulting in 
ionization and excitation of water molecules per centi- 
meter of path length is as high as ie It would be 
impractical to model every interaction. The total path 
length of an electron (the distance required to come to 
rest in the phantom) is approximated by a number of sub- 
path lengths or steps. Each step, therefore, is part of a 
random waik. At the end of each step the combined effects 
of many collisions are'taken into account. The amount 
of energy lost and, therefore, the dose deposited in the 
phantom, the amount of scattering and the position of the 
electron is evaluated for each step. The electron is con- 
sidered to have stopped when its kinetic energy falls below 
the cut-off energy which was chosen to be 0.2 MeV.* The 
remaining energy is deposited in the slab where this occurs. 
This Monte Carlo algorithm employs two major 
approximations: the continuous slowing down approximation, 
C.S.D.A., and the Gaussian scattering approximation. in 
the C.S.D.A., fluctuations of energy loss are disregarded 
and the particle is taken to be a deterministic function 
of the path length travelled. In the Gaussian approxima- 
tion, the net angular multiple scattering defieamons are 
the result of the combined effect of many small individual 
defiections. This leads to a Gaussian distribution of 


the zenith scattering angle. 


* The mean range at this energy is much less than the slab 
thickness of 0.1 cm. 
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The flow chart of the Monte Carlo algorithm is 
illustrated in Figure 73. The number of histories 
followed at a given energy and angle depends on the 
value of the energy-angle distribution as calculated by 
the subroutine, Impulse. The value of the energy-angle 
distribution, a floating-point array, is converted to 
an integer. This integer is decremented by one for 
each history followed. 

The initial position and the cosine and sine of the 
zenith and azimuthal angles of the charged particles are 
determined. It is assumed that all the charged particles 
generated in a slab are set in motion half way through 
the slab. The initial sine and cosine of the zenith scat- 
tering angle are determined by the angle bin number that 
specifies the energy-angle distribution array. The 
initial azimuthal angle is randomly distributed through 
2n radians. 

A history can end only if one of two conditions are 
met; the energy of the particle falls below 0.2 MeV or 
the particle leaves the phantom. The first condition is 
top tested in the Monte Carlo step loop (see Figure 73). 
The test to see if the particle has left the phantom is 
bottom tested. Since the phantom is semi-infinite in 
extent, a particle can only leave the phantom from the 
top surface. 

The energy lost in each step is a constant fraction 


of the energy at the beginning of the step. Such an 
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THE FOLLOWING IS REPEATED FOR ALL HISTORIES AT EACH 
ELECTRON ENERGY AND ANGLE 


ASSIGN et 
REMAINING 
| DOSE TO SLAB 


ENERGY LOST AND PATH LENGTH 


ZENITH AND AZIMUTH SCATTERING ANGLES 


IN THE ELECTRON COORDINATE SYSTEM st) 
ANGLES WITH RESPECT TO] ~ 
THE PHANTOM SYSTEM 


POSITION AT THE 
ENDS OP THe orer 


HAS THE 
EEPCTRONSTGEFT © oe 
THE PHANTOM? 


ASSIGN 
PARTIAL 
DOSE 10 
SLAB 


ASSIGN DOSE 
TO SLABS 


Figure 73. A flow chart for the Monte Carlo section of the main program 
of Buildup3. For. 
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energy loss criteria is called logarithmic spacing 
because the length of the Step, AS, decreases logarith- 
mically as a function of path length traversed through 
the phantom. Logarithmic spacing has the advantage that 
the average zenith angle changes little from step to 
Step. The smaller the fraction of energy lost per step 
the greater the number of steps required to model the 
history. This will result in a more accurate depiction of 
the history, but it will require more time to execute. A 
fractional energy loss factor of 0.1 was chosen as a com- 
promise between the contradictory demands of accuracy and 


speed: 
Bolg = 0.1 Be (5.1.25) 


where (AT). is the particle kinetic energy lost in the at 
step’ ((AT), > EPS 


T is the particle kinetic energy at the beginning 
ot the, step ae els) 


The number of steps required to complete a history is: 


Number of Steps _ 15.34 ay 
History : 0.105 


(3.1.16) 


where a (Ty + E0) is the initial energy of the particle 
(Note: the number of steps/history is an integer. The 
result of the previous expression must be rounded up to 
the nearest integer.) Therefore, a particle with an 


initial energy of 1.0 MeV will require 16 steps and a 
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particie with Ne = 10 MeV will require 38 steps to com- 
plete its history. 


The path length traversed in a step is given by: 


(AT) ; 


El 
ae 
; ee 


where Sele is the path length 


(AS), = ele) 


dT is the stopping power for the 


theustep 
The zenith scattering angle is the cumulative 
amount of scatter that occurs in a step so the angle 
is measured with respect to the direction of the parti- 
cle at the beginning of the step. 


The scattering angle is obtained by: 
1, 
6= (ot (3.1.18) 


where t is determined from random sampling (see Figure 68, 
Sections’ 3.1.2 and 31.3 for more details). 

The azimuthal angle of scattering in the step, 
with respect to the direction at the beginning of the 
step, is chosen at random to lie in any 2n direction. 

The subroutine, Geom, transforms the zenith and 
azimuthal angle in the coordinate system travelling 
with the particle to the coordinate system "attached" 
to the phantom. The algorithm will be discussed in 


more detail later (see Section 3.1.5). 
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184 


heer. ne 
maT. . 


“ire OF acts Bt withupet st ini, 
yeh univ ehoqede a st Re gana 


x ag ie 


{ dong ~~ : ; 
J o »oo “erry we 
re 
: ay 7 (thee 8 


et’ ie ind ; aye 
ot) nit resis ereiengadis =: . 
figuotht Yi Tie yrs ef SyG.ora t ci 
| ‘Wtines eat 

avissicmu ait £f oignn geaasonse: 


gigns ed? o2 Gade oe entyaae st vosdE52 se 
-ittsq si Yo spbsoorth omg = aeetpes Asiw fewenan et 


apitverth #8 08: wt, it 
bits 5 iste — 


ett mime of lui F ait i ; pio ve i a. ii iy nny A if ’ ine by ' 
“se | tnreb ste 


t sy tf t i" 
7 } ye) A 
pa! ht) ay Oy 
\ ss ae iy ; 
ney) 14 Pi 
Pie yee 


ht iP: ry ee uae 


185 


The distance travelled in the z, or depth 
direction, (the direction perpendicular to the phantom 
face with the sense directed into the phantom) in the 


step is given by: 


z = (AS)cos(¢) (del sel 9) 


where @ (cos@ + cal) is the zenith angie in the phantom 
coordinate system 


The depth traversed in the step is added to the cumuia- 


tive depth, z (the phantom surface is 2. = 0.0). 


Ss)? S 


A test is made to see if the particle has left the 
phantom. The particle has left the phantom ia tine 
cumulative depth is negative. If the particle has not 
left the phantom, the dose it deposited in slabs 
through which it traversed are assigned. If it did 
leave the phantom, the dose in the slabs before it 
left are deposited. The subroutine, Assign, assigns 
the absorbed dose to the phantom slabs (see Section 3.1.6 
on Assign). If the particle did not leave the phantom, 
the control of the calculation is passed to the test 
to see if the energy is above the cut-off energy. If it 
did leave, control is passed to the statement control- 
ling the number of histories at each initial energy 
and angle. If the particle does come to the end of its 
range its remaining energy is assigned to the siab it 
last entered. 

After all the histories are followed, for all charged 


particles’ energies and angles generated in all the slabs, 
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the dose is deposited as a function of depth and then 
normalized to the maximum dose received. The normalized 
dose as a function of depth is printed as the output of 
the calculation. 

The Monte Carlo program was run, without the charged 
particle set in motion part of the program (in other words, 
withouth Subroutine Impulse) with a set of initially 
monoenergetic, monodirectional electrons. The energy of 
the electrons was varied between 2.0 MeV and 10.0 MeV. 

The penetration profile obtained was compared with results 
of a well-established analytic technique due to Spencer 
using the moments method (50). Spencer's method also 
employed continuous slowing down. This comparison is made 
in Figure 74. 

Scattering Angle Transformation 

The subroutine, Geom, transforms scattering angles 
from the coordinate system travelling with the particle 
to the coordinate system fixed to the phantom (see Figure 
75). The direction with respect to the particle reference 
system in polar coordinates is specified by a zenith angle, 
o, (¢ > g), and an azimuthal angle, ~, (y>r). 

Directions in the phantom reference system are also 
in polar coordinates. The zenith direction in the 
phantom system is « and the azimuthal angle is 8. It is 
useful to specify the direction sine and cosines. The 
direction of the particle in the phantom system at the 
beginning of the step is specified by cos « (cos < + oldcal) 


or sin « (sin « + Oldsal) and cos g (cos g + oldcbeta) or 
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Figure 74. Comparison between the Monte Carlo program of Buildup3.For 
and Spencer's analytic formalism to predict the dose 
deposited by various monoenergetic electrons in water. 
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Sin 8 (sin 8 > oldsbeta). The direction of the particle 
in the phantom system at the end of the step (after 
scattering has taken place) is given by cos «“ 
(Cos i= & Cal) or Sin = (sin-= .+‘saljeand Cosme. 
(cos 8“ socbeta) or San’ 8° 1(Sin. Gr Saepeta)n Tihesre— 
lationship between the directions in the particle and 
the phantom system are shown in Figure 75. 

To determine the position of the charged particle 
inside the phantom, a transformation must be made to the 
coordinate system fixed with respect to the phantom. 


The transformation is given by (51)*: 
cos(«*) = sin(¢)cos()sin(«<) + cos(¢)cos(«) (Se be2) 


sin(=*) = V1 - cos2(«*) (3.1.21) 


sin(o)sin(v) awe 


Ssim@>.-C p= Se 


cos(o) - cos(«<)cos(«“) ies) 


cos(8 ge. = sin(«<)sin(<“) 


8° can be found by taking the inverse sin or Cos of the 


above expressions but it 1s more convenient to define 


sin(g*) and cos(8~) as: 


sin(8*-8)cos8 + cos(8~-8)sinB (Selze) 


sin(6~) 


cos(8*-8)cos8 - sin(8~-8)siné (Zo) 


cos(8~) 


Having obtained sin«", cos=-, sing” and cose ,, these 


values are passed back to the main program. 


* The error in the formulation of Equation 3.1.22 in Bruce 


and Johns (51) has been corrected. 
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Phantom Coordinate System xX¥Z 


Particle Coordinate System x,y,z 


Relationship between scattering angles in the particle and 
phantom coordinate system. 
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Assigning Dose 

When called, the subroutine, Assign, deposits the 
required dose in the slab or slabs traversed between 
the beginning and end of the step. Dose is assigned 
for each of the following three situations: The particle 
energy is greater than or equal to the step cut-off energy 
(0.2 MeV), the particle energy is less than the step 
cut-off energy at the end of the step or the particle 
has left the phantom during the step. 

The simplest possible case and the most often 
occurring case is when the particle does not leave the 
slab it was in at the beginning of the step. The total 
dose received by the slab, up to and including the step, 


Dose*(k), (k + kay), should be given by: 


Dose’ (k) = Dose(k) + apt (376) 


where Dose(k) is the total dose received by the slab up 
to and not including the step. Since the phantom is 
semi-infinite, the mass of the slab, Pi dah? PS inriniee: 


This difficulty is eliminated by redefining Dose*(k) as: 


Dose~(k) = Dose(k) + (AT)s (5 lean) 


Now Dose’(k) is not equal to the absolute dose. However, 
since the slabs have the same thickness, Dose~(k) is 
directly proportional to the absolute dose. The final 


output of the program is the dose at each slab normalized 
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to the maximum dose received by any slab in the phantom, 
so absolute dose is not required. 

If the particle traverses more than one slab 
during a step, the energy lost is assigned to the slabs 
in proportion to the fraction of the path length spent 
in each slab.* 

If the particle traversed completely through one or 
more slabs, the dose deposited in those slabs will be: 


~ NCAT 
Dose’(k) = Dose(k) + emai kaci (Zeulee 3) 
j 


where 0.1 cm is the slab thickness. The denominator in the 
last term represents the number (floating point ) of 
slabs traversed. 

If the particle traverses only part of two slabs 
the fraction of the total dose deposited, Dose Fraction, 
(Dose Fraction > Dosefrac), in one of the slabs due to 


a step is calculated. The remaining energy is deposited 


in the other slab. 


Similar procedures are followed if the particle 


* The situation for which this will occur is at the top 
of the phantom and the particle has a, large initial 
kinetic energy. For example, the distance travelled 


by a 15 MeV electron in iS Tipsewscep (AS) is given 


by: 
7 1.5 MeV 


(43), 29396 MeV 


Therefore, the dose imparted to 7 slabs will have to 
be calculated. The mean square angle of scattering 


will be: 


(02)! = /3295 x 10-2 rad@/cm) (0.65 cm) = .16 rad 


which is about sine 
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leaves the phantom, except that dose is not assigned 
for that portion of the path outside the phantom. 
Results 

The program was run to model 6 MV and 15 MV 
accelerators under narrow beam (or small field) and 
broad (or large field) beam conditions. The program 
was designed to calculate the energy-angle spectrum of 
charged particles set in motion in thin phantom slabs. 
The Monte Carlo part of the program follows the charged 
particles through the rest of the phantom. This pro- 
cedure is common to both the small field and large 
field approximations. Figures 76 and 77 illustrate 
the dose delivered to a phantom due to charged particles 
set in motion in a 0.1 cm slab of water at the surface 
of the phantom by a 6 MV and 15 MV photon beam, 
respectively. To a good approximation, the dose due to 
charged particles set in motion in a thin slab attenuates 
exponentially with distance from the slab. The charged 
particle energy absorption coefficient, os can be 
found by finding the slope of the line of a graph of 
Qn(1.0/D(x)) versus depth where D(x) is the relative dose 
at a depth, x (the maximum dose = 1.0). The charged 


1 


particle energy absorption coefficient is 3.1 cm” for 


6 MV and 1.15 cm for 15 MV photon beams. 
In the broad beam, or large field approximation, 


the scattered photons generated in the phantom are not 
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Figure 76. Distribution of dose ina Ss 
. he top 0.1 cm slab by a 6 MV 


electrons set in motion in t 
photon beam. 
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Fioure 77. Distribution of dose in a semi-infinite water phantom due to 
ah electrons set in motion in the top 0.1 cm slab by a 15 MV 


photon beam. 
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lost from the beam. Scattered photons have a spectrum 
that has a similar shape, but an effective nominal 
energy lower than the primary beam (36). It was also 
seen in Section 2.13 that the depth of d.39 of the 
scattered radiation is less than that of the primary 
radiation. The primary photon spectrum tends to harden 
as it penetrates into the phantom. Lower energy 
scattered photons tend to replenish the low energy 
photons lost from the spectrum due to hardening. Large 
field geometry is approximated by assuming that the 
primary photon spectrum is unaltered as the beam 
penetrates the phantom. 

With the primary photon spectrum remaining the same 
as a function of depth in the phantom, the penetration 
characteristics of the charged particles set in motion 
remain the same. The number of charged particles set 
in motion at any depth decreases exponentially with 
depth. The appropriate attenuation coefficient is the 
mean energy absorption coefficient which is an energy 
fluence weighted average normalized to the total energy 
fluence of the beam. Mathematically, the mean energy 
absorption coefficient is given by: 


ZF (hv) | Men) j 


(a as 29) 
v 


Equation 3.1.29 has the same form as Equation 1.5.10. 


The mean attenuation coefficient for the 6 MV spectrum 
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is 0.029 cm ~ and 0.021 cm} for the 15 MV spectrum. 


The dose deposited by the charged particles set in 
motion attenuates exponentially in the phantom and the 
primary photon beam attenuates exponentially. The 
assumption, therefore, that large field geometry can be 
modelled by keeping the primary spectrum constant with 
depth is identical to the build-up model presented in 
Section 1. 5. 

The method of computing the dose in the build-up 
region for broad beam geometry is to attenuate the 
dose due to charged particles set in motion in a thin 
slab at a depth, x, by a factor » Hen, This is done 
for 100 thin slabs (each 0.1 cm thick). The dose re- 
ceived at any point in the phantom is a sum of the dose 
reaching the point from the charged particles set in 
motion from all the slabs above the point. 

This procedure is done by the program, Addose. For. 
The input for this program is the thin slab charged 
particle set in motion dose penetration profiles shown 
in Figures 76 and 77. The depth of maximum dose, 

d ax? as predicted by the large field model computed by 
Buildup3. For and Addose.For, is the same as the ae 
point predicted by Equation 15.112 ey hesoucpute Lor. 
Addose. For is normalized by dividing the dose by the dose 
The result is the broad beam TMR curve in the 
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build-up region. Appendix 8 contains the listing and 
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documentation for Addose. For. 

Small field geometry assumes that all the scattered 
photons produced by the beam in the phantom do not 
interact again. ‘This is modelled by hardening the 
beam. Every photon that interacts in the phantom is 
removed from the number spectrum. ‘This is accomplished 
by the subroutine of Buildup3.For called Impulse (see 
Equation 3.1.16). Small field geometry requires more 
computer time to model than large field geometry because 
the charged particle set in motion energy-angle dis- 
tribution has to be recalculated for each slab. More 
importantly, each charged particle set in motion in the 
slab has to be followed through the rest of the phantom. 
The relative dose was calculated to a depth of 10 cm so 
the small beam geometry requires approximately 100 times 
(10 cm/0.1 cm) as much time to model as the large beam 
geometry. 

The build-up curves of the small and large field, 
together with the experimental results from Section 
2.11 are shown in Figures 78 and 79 for 6 MV and 15 MV 
photon beams, respectively. As was found experimentally, 
there is little difference between small and large field 
geometries in the build-up region. There is good agree- 
ment between the computer modelled and experimental 


build-up curves. 
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Comparison between experiment and computer modelling of 


the phantom generated TMR curve at 6 MV. 


Figure 78. 
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d TMR curve at 15 MV. 


Comparison between experiment and computer modelling of the 


phantom generate 


Figure 79. 
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3.2 Nature and Source of Contamination 

The definition of contamination stated in the introduction implies 
that radiation reaching the phantom, other than primary radiation*, is 
to be considered as contamination. Primary radiation is the photon 
component of the beam that is directed away from the source. 

The utility of this definition of contamination is that it allows 
one to determine the dose deposited by primary photons in a phantom. 
This is called the phantom generated dose. The phantom generated 
dose is deposited by charged particles set in motion by primary 
photons and scattered photons produced in the phantom. Therefore, 
contamination dose is exogenous whereas the phantom generated dose is 
endogenous to the phantom. Scatter, measured by an increase in dose 
with increasing field size using a full phantom in the build-up region, 
cannot be confused with contamination generated outside the phantom. 
Changes in the phantom generated dose can only be due to changes in 
the phantom size, density or configuration. A modification of the 
phantom in these respects will not affect the nature of the contamina- 
tion component (although it will affect the interaction of the con- 
tamination component in the phantom). 

The TMR build-up curves in Section 2.1 demonstrate a field size 
dependence, especially at or near the phantom surface, greater than 
can be expected from lateral scattering of photons and electrons 
generated in the phantom. 

Adding small additional thicknesses of accessories increased the 


TMR in the build-up region. As the thickness of accessory increases, 


* The dose due to charged particles set in motion by neutrons is not 
considered. 
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the amount of contamination produced by the accessory increases. 
However, contamination produced at shallow depths within the accessory 
is shielded by deeper layers of the accessory. At an equivalent thick- 
ness less than 3.2 cm of Lucite, there is a balance between production 
and absorption of contamination. 

The surface dose displays the characteristics of electron con- 
tamination. High atomic number ''filters'' can reduce the surface 
normalized dose due to an accessory when placed beneath the accessory. 
When placed above the accessory the high atomic number material did 
not alter the surface normalized dose. The number of interactions 
produced by the high Z material above and below rhe accessory are 
approximately the same. The mass stopping power is relatively in- 
dependent of atomic number and, indeed, the mass stopping power of 
Lucite, for example, is slightly higher than lead (5). Therefore, 
the only possible difference between the configurations of high Z 
material above and below Lucite resulting in a reduced surface nor- 
malized dose is the difference between the amount of scattering of 
contamination electrons outside the beam. This is not surprising con- 
sidering that the scattering power of lead is nearly an order of magnitude 
higher than water* (see Appendix 4). Therefore, nearly as many contamina- 
tion electrons emerge from the bottom surface. of the jicgh,.2 (i Gaiter’ as 


emerge from the bottom surface of the accessory, but fewer of the elec- 


trons are forward directed. Therefore, a high atomic number material 


placed in the beam should be called an electron scatterer rather than an 


electron filter. 


The distance from the source distribution and lateral distribution 


* The scattering power of water and Lucite are within 10% at energies 
less than 15 MeV (5). 
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of the surface dose for large open fields at 15 MV are typical of 
electron distributions. The distance distribution of the surface dose 
falls off as 1/r* from a point 12 cm below the target which corresponds 
to the position of the field flattening filter. The field flattening 
filter has, at its thickest point, 5.4 cm of stainless steel in the 
beam. It is the only component present in an open field which has a 
sufficient thickness to produce a maximum fluence of contamination 
electrons*, so it is likely to contribute to the dose due to contaminant 
electrons. 

The lateral surface distribution .is Gaussian. This indicates that 
the surface dose is due to electrons and indicates that the electrons 
could have been produced from a smail source. Therefore, the surface 
dose at 15 MV for large fields behaves in accordance with Equation 
L334. 

The situation at 6 MV is quite different. The assumption that 
the surface dose decays as i/r? at large field sizes is not valid.** 
Additionally, the lateral distribution is not Gaussian. The surface 
dose inside the field is iarger than can be expected from a normal 
distribution and the surface dose outside the field is smaller. Both 
of these are contraindications that localized source is responsible 
for most of the electron contamination for large fields at 6 MV. 

The source ot electron contamination at small fields is better 


understood due to magnetic separation of the electrons produced by 


1 1 1 ill not 
* A Lucite accessory with a thickness greater than 3.2 cm wi 
contribute to greater contamination. The equivalent thickness of 


steel is about 40 cm. 


: i i from an apparent source of contamination 
** A 1/r“ dose decay with distance trom an appare : cor 
is yi Sones ret not a sufficient condition for identifying the 


source of contamination unequivocally. 
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two regions of the beam. The surface normalized dose was measured 
for various conditions of fieid size (the largest field size was 
10 cm x 10 cm) and source to probe distance inside the field with 
the magnet in place and not in place at the accessory holder position. 
Air contributes most of the surface dose at 6 MV and 15 MV for these 
field sizes. At 15 MV, however, the increase in the dose, as a func- 
tion of field width due to electron contamination produced between 
the target and magnet, suggests that at a field size of about 15 cm 
x 15 cm for a source to probe distance of 100 cm, this contribution 
will be greater than the contribution due to air.* At 15 MV, when 
the SPD=100 cm for a field size of 30 cm x 30 cm, (which was the 
condition of the measurement of the distance distribution in Section 
2.4), the contribution due to air will be a minority component. 
Therefore, the results obtained at small field sizes does not 
contradict the results at large field sizes at 15 MV. At 6 MV the 
surface dose for large field sizes is mostly due to air contamination 
(with the possible exception of small source to probe distances). This 
is why the source of electron contamination could not be localized 
to a small region for the 6 MV large field beam in Sections 2.4 and 
ogg 

Qualititively, the SPD dependence of the contamination components 
are easy to explain. The greater the thickness of air intervening 
between the magnet, the greater the dose due to contamination produced 
in air. Conversely, increasing the distance between the region of 
production of contamination between the target and magnet will result in 


4 smaller normalized dose. 


* Based on a linear extrapolation of the increase in dose between a 
field width of 7 cm and 10 cm. 
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Plotting graphs of the square root of the reading due to contamina- 


tion produced between target and magnet will determine the position of 


the source of contamination if the source is confined to a small region. 
At small fields at both 6 MV and 15 MV with the cross-hair tray in place, 
the source appears to be the cross-hair tray. At 15 MV, with the cross- 
hair tray not in place, the source appears to be the collimators. This 
agrees with the observation that the surface dose appears to be produced 
by the field flattening filter at 15 MV for large field sizes. At small 


field sizes the contamination produced by the field flattening filter is 


largely blocked by the collimators. 

A general conciusion can be hypothesized by these findings. 
Electron contamination measured at a point tends to be produced by 
mass close to that point. Contamination produced further away iS 
stopped or scattered by mass lying close to the measuring point. 
Therefore, air is a significant contributor under most conditions. 
Accessories and the cross hair tray when present contribute. 
Collimators at small field sizes produce some contamination*. The 
field flattening filter can be the major source of electron contamina- 
tion at large field sizes (therefore not blocked by the collimators) 
if the cross hair tray or accessories are not in the field. 

All these possible sources of electron contamination result in 
a greater magnitude of electron contamination reaching the probe as 
a function of field size. A greater field size results in a greater 
interaction volume for air, accessory and cross hair tray sources. AS 
the collimator jaws are opened a greater amount of contamination pro- 
duced by the field flattening filter emerges through the collimator. 


* At small field sizes the perimeter of the collimator compared to the 
field area is much larger than the same ratio at large field sizes. 
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opening. 

The key to the separation of the contamination from the phantom 
generated dose is to know the contribution of contamination to the 
dose at one point. Except for the photon backscatter* which can be 
independently measured, the surface dose is due to contamination 
electrons generated outside the phantom. Knowing the contamination 
contribution at the surface and the penetration characteristics of 
the contamination component, the phantom generated dose can be ob- 
tained by the subtraction of the dose due to contamination. 

The same contamination penetration curve was used to determine 
the amount of contamination for each field width, SPD and accessory 
configuration. This assumes that different beam geometries or the 
presence of accessories in the field do not change the nature of the 
contamination. Justification for this assumption was provided at 
15 MV by determining the contamination penetration curve at different 
SPD's and accessory configurations. The contamination penetration 
characteristics were within 5% (see Table 3). The surface dose is 
usually less than 40% of the maximum dose so an uncertainty in the 
correct amount of contamination at the surface will be less than 2% of 
the maximum dose. At deeper depths the uncertainty will be even less. 

The validity of the contamination curve at 15 MV was verified in 
a number of other ways. The shape of the contamination penetration 


curve is characteristic of mixed electron and photon contamination. 


The surface contamination is established as being due to electrons. 


These electrons must penetrate to some depth beneath the phantom. 


* Some of the photon backscatter is due to scattered photon backscatter, 
but it need not be taken into account. The photon backscatter at most 
makes up about 20% of the electron contamination which usually makes 


up no more than 30% of the phantom generated dose, therefore, 


i i a 
ter will likely contribute less than 6% 
scattered photon backscat OT Re en 


photon backscatter. Photon backscatter 1 
scattered photon backscatter is negligible. 
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The peak and rapid fall off at shallow depths are sraygrawseusiure of 
electrons. The slow fall off at depths greater than 2 cm is character- 
istic of photons. Therefore, the contamination radiation appears to 
be a combination of an electron contribution primarily at shallow 
depths and photons contributing primarily at greater depths. 

The derived contamination penetration curve, using the increase in 
dose in a limited phantom versus the increase in field width (greater than 
the phantom dimensions) at 15 MV, is similar in shape to the contamination 
curves measured outside the field. The difference in shape (see Figures 
63 and 31) can be attributed to a difference between the relative amount 
of contamination electrons and photons present inside and outside the 
field boundary. Figure 27 demonstrated a Gaussian lateral distribution 
of the surface dose. Consequently, for a 30 cm x 30 cm field size there 
were about 45% fewer electrons present at 18 cm from the central axis 
as were present at the central axis. This largely accounts for the 
difference observed between Figures 31 and 63. 

The most convincing evidence that the phantom generated build-up 
curves are correct is their agreement with computer modelling obtained 
independently (see Figure 79). The computer model also verifies that 
the phantom generated TMR curves are independent of field size. 

The validity of subtracting the contamination component from the 
total build-up curve to yield the phantom generated build-up curve at 
6 MV was tested as well. The contamination penetration had its maximum 
at the surface and decayed eridty ath depth at shallow depths (see 
Figure 64). The curve levelled off at deeper depths. This, as was 
found at 15 MV, is indicative of electrons contributing at shallow 


depths and scattered photons contributing at deeper depths. The contam- 
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ination curves derived inside the field (Figure 64) and measured outside 
the field (Figure 32) are in approximate agreement if the lateral 
surface dose distribution is taken into account. Approximately 60% 
fewer electrons are present 3.cm outside of a 30 cm x 30 cm, «6 MV field 
as compared to the central axis. This accounts for the difference in 
the shape of Figures 64 and 32. The computer model of the phantom 
generated dose as a function of depth agreed with the experimentally 
derived phantom generated dose build-up curve at 6 MV (see Figure 78). 

The computer model established that the dose due to electrons set in 
motion in a thin phantom slab decays exponentially in the phantom as a 
function of depth from the slab. Therefore, Equations 1.5.8, 1.5.11 
and 1.5.16 should be valid descriptions of the dose in the build-up 


region*, 


Values of ut of 3.1 cm’! and 1.15 cm! for the mean energy absorption 


coefficient of electrons set in motion by 6 MV and 15 MV photon beams, 


respectively, were found (see Section 3.11) he value or u was 


i} 


found from Equation 3.1.29 to be 0.029 cm 4 for 6 MV and 0.021 cm” for 


15 MV photon beams. 
The due predicted by Equation 1.5.11 is 1.5 a@m for 6 MV and 3.6 cm 
for 15 MV, respectively. This compares with the experimentally determined 


d values (after the contamination component has been removed) of 1.6 cm 
max 


for 6 MV and 3.5 cm for 15 MV. 
Table 19 lists the TMR dose for 6 MV and 15 MV beams in the build-up 


region as predicted by Equation 1.5.16. This agrees well with the 


phantom generated build-up TMR dose. Therefore, the simple differential 


Equation 1.5.6 is physically acceptable and accurately predicts* the dose 


in the build-up region received by a homogeneous unit density phantom. 


* The model fails to predict the surface dose due to backscattered 


photons. 


207 


ahietue berwansi bas (hd saugilt), biel ath ot 
fervetak ori2 23 TAS te comet 
208 vieanmiandagA, steaeses ost ede aio 

bisit Wie 9m Of & wi» om 8 dom a | 
fr sonerstiib odd tol etmeagR aint - ies, dersim> ft oo he 
maumiq exit to aber AGT oft 5 bas ba sous 30 99 
iifetmeiteyss elt dtiw beetgs. dsent, Fee Mokau? 5 aaa 
(seat see). a te avi gerbil ius sant Setar, 
at tox enortuels 6? sub seah AM Tags hoiles fdmtes fence y3% agit 4s , 
e ze mutteda Gon i | viLebtekecge -eveoeh ¥en ae ‘ 
Hogek het enokcmupeh oadbantt dale ons most diqeb te 
qp-blind afd of seeb off to snoltgirwesb bHev od tyode a 
‘nial j Bs 
crproeda ‘trend masa ott 14 ap ar, i brs * mm, Le ae halt a ia 
2aped ootodg we Th Bes Ve ag nokton nt tee an a 
oo 0 YO soteveal RL mekgnee 228) be oTaw <) _ 
io 50.0 bake WLR mot fap WS. a as sb io a 

m) 0.0. bes VM aw "ria, fie 6 ak ae Medak eopemns vd iosatiiorq sol 2 

Loginrret$b 1) Pepsiabt saps Sele Mahe ad atdt Dunviciotiamh! ave me = 
48 2. to. (Bavanuet nestl 2ai Pestmagmes aaizeainstin edz ‘Oe 
| jo MA Rt a8 mw 2 ee 
co-titzig, a tse OW ‘siA pooh wer ads Fig a sais 
ods Ani“ Dew le od eat, er. 16 be ate wr sm 
tpi aneetRi ocenitz ans Ce ee) T cu-blind botereneg mod 
eeob at resnibeng, er nn si 
es ‘entails theo 2 “HS TAE 
ieradteseiond: ot ws 8 % ; 


a ae ene a; ne Te ae Wee Y i enh, ual v 


208 


Table 193 Tissue Maximum Ratios in the Build-Up 


Region as Predicted By Equation 1.5.16 


6 MY 1S MU 
Depth 
(com) TMR TMR 
8.88 @.808 @.080 
S216 @.408 @.181 
On ac @.654 6.331 
0.48 @.882 @.455 
8.64 8.8908 @.557 
4.95 @.3971 8.7288 
ie aie @.996 @.814 
2 es dl 9.995 Oe se 
2.54 = @.982 
3.18 = 8.998 


3.81 - 1.098 


ise of ime ene. i 


BLie.2 aoizauns ya 


ye 2 = 
RMT 
850 .G ae 
601. .& 
rée .8 
Zee .8 
{e728 
Sey 8 
hIG.8 a? is 
ser . * 
sae. 
ce = . 


ath oie 
| ik 
ir obi ane were a , bie 
: | as yn ae he 
An - : ai a) . 
a Tare it 


ry); ar a 
an r ih i rp) ny rhs 
ie a is ‘i 7 th 

y ah : rai ae a nth oat im, ya Pe 


hae ue aaa mn . 7. 


cl « 
ie r 
a. 
* it “a 

i 


. : 4 
mh see. “in? « P iit 


; Lo 
i! pas ne i / oy 


i _s : nt . Pie arp ; ii 


209 


region received by a homogeneous unit density phantom. 
The contamination electrons are set in motion by matter between 
the target and magnet by the same processes that set electrons in 
motion in the phantom. The major difference between contamination 
electrons produced between the target and magnet that reach the phantom 
surface and electrons set in motion by primary photons interacting with 
the phantom near the surface is that the contamination electrons will 
be more forward directed. Those contamination electrons that were set 
in motion with a large zenith angle will have left the beam. Figure 72 
illustrated that the forward directed electrons that are set in motion 
by a 15 MV beam have a higher average energy than those that are set in 
motion with large zenith angles. Therefore, contamination electrons have 
a higher average energy compared to electrons set in motion in the phantom. 
Having a higher average energy results in greater penetration of the con- 
tamination electrons compared to phantom generated electrons. Additionally, 
contamination electrons that are set in motion further away from the 
phantom are more likely to have been produced with a smaller zenith angle 
than contamination electrons produced closer to the phantom. Therefore, 
electrons produced further away from the phantom are likely to be more 
energetic* and, therefore, more penetrating than those produced close 
to the phantom. 
Figure 80 compares the relative dose imparted to a phantom from 
electrons set in motion in the top 0.1 cm of a phantom compared to the 


electron contamination penetration curve from a 15 MV photon beam and 


* The slowing down of electrons in air is insignificant. The stopping 
power of ve for 3 MeV electrons is about 1.8 MeV cm2/g (5). One , 
hundred centimeters of air has a mass thickness of about 0.13 g/cm+ so 
an electron will lose about 0.2 MeV or 7% of its energy in traversing 
this distance. The range of 3 MeV electrons in air is about 11 m. 
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Comparison between various penetration curves at 15 MV. 


Figure 80. 
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the electron contamination produced between the target ana magnet. 
The contamination penetration curve contains electrons as well as 
scattered photons so the comparison cannot be made fairly. In Section 
2.13 the scattered photon build-up curve was obtained. The electron 
contamination component can be obtained by subtraction of the 
scattered photon contamination component. An assumption is made that 
there are no contamination electrons penetrating to a depth of 4.4 cm 
of polystyrene so the dose at this depth is entirely due to scattered 
photons. Figure 80 shows the result. The electron contamination 
component is more penetrating than the electrons set in motion in the 
phantom at shallow depths. This is because the electron contamination 
component contains fewer low energy electrons compared to the electrons 
set in motion in the photon which are stopped at shallow depths in the 
phantom. The electron contamination component produced between the 
target and magnet is more penetrating than the total electron contamina - 
tion curve. The maximum range of the contamination electrons is 
approximately equal to the maximum range of the electrons set in 
motion by the phantom. This is because the maximum energy of contamina- 
tion electrons produced at any position in the beam is the same as the 
maximum energy of electrons set in motion in the phantom. All the elec- 
tron penetration curves in Figure 80 have their maximum at about 1.5 m 
depth in the phantom. 

Figure 81 is a comparison between the dose due to electrons set in 
motion in a thin slab at the top of the phantom and total dose due to 


contamination for the 6 MV beam. The contamination curve has not had 


the scattered photons removed. By comparing the total contamination and 


electron set in motion curves between 15 MV (Figure 80) and 6 MV (Figure 
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Comparison between various penetration curves at 6 MV. 


Figure 81. 
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81) it can be seen that the 6 MV curves would be very ence had the 
scattered photon component been removed. Therefore, the electron con- 
tamination component at 6 MV has similar penetration characteristics 
to electrons set in motion at the surface of a phantom implying that 
the contamination electrons are produced close to the phantom. This 
confirms the results of Section 2.10 which indicated that most of the 
contamination at 6 MV is produced in air Nene the magnet and phantom. 
Figure 81 also illustrates the penetration curve of electrons 
produced between the target and magnet at 6 MV. The penetration of this 
fraction of the electron contamination is much greater at shallow 
depths than the electrons set in motion at the surface of the phantom, 
however, this component contributes very little to the total contamina- 
tion at 6 MV. As at 15 MV, the maximum range of electrons set in 
motion at the surface of the phantom and contamination electrons 
produced between the target and magnet are almost identical. 
Decreasing the source-to-probe distance increases the fraction of 
detected contamination produced between the target and magnet and 
decreases the amount of contamination produced by air between the magnet 
and phantom that arrives at the probe. Therefore, the relative proportion 
of electron contamination that is produced by these sources changes as a 


function of SPD. Since these sources each have different penetration 


characteristics measured at the same SPD, the penetration characteristics 


of the total electron contamination component could change as the SPD 


changes. However, as the source to probe distance decreases, more 


electrons produced between the target and magnet are going to be 
arriving at the probe produced with larger zenith angles. As the probe 


gets closer to a source of electron contamination more lower energy 
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electrons will arrive at the probe so the contamination component will 


be less penetrating at shallow depths. More lower energy electrons 


arriving from a source between the target and magnet will tend to com- 
pensate for fewer electrons produced in air. Therefore, the total 
penetration characteristics of contamination electrons likely does not 
change significantly as a function of source-to-probe distance. 

At depths approaching d ax? most of the contamination dose is due 
to scattered photons. The Compton effect produces a scattered photon 
for every recoil electron produced. Therefore, the sources of electrons 
are also going to be the sources of scattered photons. 

Nilsson and Brahme, using computer Monte Carlo modelling, found 
that the field flattening filter and the collimator are the main 
sources of scattered photons (36). This can be qualitatively explained. 
The mean interaction distance for the scattered photons is of the order 
of the inverse of their mean energy absorption coefficient, lh peer 
Section 3.2 illustrated that one could define an effective energy 
absorption coefficient for the dose deposited by electrons set in 
motion. The mean interaction distance is of the order of the inverse 
of this coefficient which is /ut- Scattered photons will emerge from 
a depth from a source of the order Te while contamination electrons 
will emerge from a depth of the order 1/7,+. Both the collimators 


and the field flattening filter have thicknesses greater than W/Tan 


sO many more scattered photons will emerge compared to contamination 


electrons. Sources of contamination with a thickness of the order of 


1/i_+ , like accessories and air, will produce approximately equal 
e 


numbers of contamination electrons and scattered photons in the beam. 


Contamination electrons produce a higher maximum dose illustrated 
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by the surface dose of Figures 63 and 64 than scattered photons 
illustrated by the dose plateau in the same figures. The ratio of the 


electron energy fluence produced by scattered photons deposited at any 


point is from Equation 1.5.8 and is approximately equal to: 


acta ( y/e)x = Wgt/0)x 
Gisb)  cusiaiucl kyle 
ree) e Y 
-(Hy/p)dmax  - Het/0) dmax 
At ans the term (e =e ) is approximately equal to 
iyee Simece ms is much less than ut : 
anc hae) eer 
peer = (352.5) 
Wo uot 


It has been shown previously that contamination electrons penetrate 
somewhat further than electrons set in motion in the phantom by scattered 
photons, however, they are both of the same order of magnitude. There- 
fore, the relative energy fluence of contamination electrons compared 


to scattered photons is approximately: 


Vee contamination Dose (d=0) 3 


Bical s OSCR OE ait ree paren ee iar aa 

a 

Cu) scattered photon Dose (ca scattered photon) ue? 
(ec) 


Therefore, despite the fact that the maximum dose due to contamination 

electrons is greater than the maximum dose due to scattered photons, 

the energy fluence arriving at the phantom surface of scattered 

photons is greater than the energy fluence due to contamination electrons. 
The greater energy fluence of scattered photons compared to con- 

tamination electrons is due in part to thick sources like the 


collimators and the field flattening filter. Additionally, electrons 
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more readily multiple scatter* out of the beam compared to photons. 
The shift in q ax? especially when accessories are present in the 
beam, can be attributed to electron contamination. Since the build-up 
curve for scattered photon contamination is similar to the primary 
photon build-up curve, scattered photons are unlikely to contribute 
Significantly to the shift in qe The peak in the contamination 
penetration curve due to electrons occurs at very shallow depths so 
increasing the magnitude of electron contamination results in a cies 
shift towards shallower depths. This also explains why the shift in 
oe is accompanied by an increased surface dose. These findings agree 


with Biggs and Ling (41) who observed that a shift in eS. did not occur 


when electrons were swept from the beam. 


* For example, a 3 MeV electron in air has a mass angular scattering 
power of 0. 602 padians* = cm2/g (5) which, after passage through 100 cm 


; O 
of air, corresponds to a mean scattering angle of 16.0%. 
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3.3 Clinical Uses of the Results 

Low energy photon beams (isotope sources with an energy less than 
ex .equal to mc? or accelerator sources less than about 1 MV) produce 
a maximum dose at the surface of a patient. If these photon sources 
were used to treat tumors at a depth under the patient's skin surface, 
Sein reactions could occur. 

Skin reactions may be minor such as erythema (reddening of the 
skin) or more serious such as moist desquamation (peeling of the skin 
with blistering). The type of reaction, if it occurs, varies con- 
Siderably as a function of dose from patient to patient (7). Skin 
reactions are similar to sunburns in morphology and histology, however, 
the onset of the ionizing radiation induced skin reaction occurs several 
days after exposure rather than a few hours. 

Few skin reactions occur to patients undergoing megavoltage 
radiation therapy. The main reason is that the maximum dose is received 
at a depth well below the skin. The dose due to contamination electrons 
reduces this skin sparing effect. 

Another effect produced by contamination electrons and scattered 
photons is an elevated dose outside the field boundary. A situation 
where this may be clinically important arises when the eyes lie just 
outside the field. The lens of the eye is susceptible to the formation 
of radiation induced cataracts. Cataracts are likely to form if the 
total dose received over 30 days is 1500 rads (52). 

The amount of contamination increases with accessories in the 
field, greater field size and smaller SPD. Skin reactions or high dose 


Outside the field may result even at megavoltage energies if two or more 


of these conditions are encountered. 
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Passing the beam through the poles of a magnet could be effective 
in reducing the electron contamination arriving at a patient. To ensure 
that few electrons are produced in air between the magnet and patient, 
the magnet would have to be placed as close to the patient as possible. 
In a clinical situation the contamination electrons should not strike 
the patient outside of the treatment area, therefore, the electrons 
should be swept through at least 90°. A magnetic field of nearly 
1 Tesla would be required. Electron contamination increases as a 
function of field size so the pole gap should be as wide as possible. 
The field gradient should be steep enough to ensure that the dose 
distribution, due to electrons set in motion in the patient is not 
affected. These requirements could only be met with an electromagnet. 
In addition, such a magnet would likely have to be liquid cooled. 

When an accessory is in place, the surface dose inside the field 
can be reduced by about 20% if a lead (or other high atomic number 
material) sheet is placed underneath the accessory. The ''filter" 
scatters electrons out of the field so a reduction of surface dose in- 
side the field would be concomitant with an increase in the surface 
dose outside the field. The 'filter'' would then be clinically acceptable 
only if the increased dose outside the field would not cause complications 
and if the reduced dose inside the field significantly reduced the 
probability of the occurrence of skin reactions. 


A useful. immediately applicable feature of the removal of the 
b 


contamination component is that the TMR in the build-up region is 
independent of field size, SPD and presence of accessories and, therefore, 


is ameniable to data storage in computers. Instead Of Wisting allo 


the values of the TMR curves in the build-up region for field widths, 


SPD's and presence of accessories applicable to treatment planning, only 
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two TMR curves; one, the phantom generated TMR curve and the other, 

the contamination TMR curve and the surface dose as a function of these 
conditions need be listed. The surface dose dictates the magnitude of 
the dose received in the build-up region due to contamination. 

The separation of the contamination component of dose and the 
phantom generated dose may prove to be important for treatment planning 
in the build-up region in the presence of inhomogeneities. The treat- 
ment planning algorithms most often employed (developed by Dr. Jack 
Cunningham) separate the contribution to the dose to a point into 
primary and scattered photon contributions. In the build-up region, 
the phantom generated scattered contribution can be confused with the 
contamination contribution. 

The primary dose at a given depth is determined by extrapolating 
a curve of TMR versus field width to the y-axis. The TMR value at the 
y-axis is called the zero-area TMR value at that depth. The difference 
between the TMR value at any other field width is assumed to be due to 
scatter generated within the phantom. 

The usual measure of the amount of phantom generated scatter con- 
tributing to a point at energies of 6 MV or higher is the scatter 
maximum ratio or SMR. The scatter maximum ratio is the difference be- 


tween the TMR at a given field width and depth and the zero-area TMR at 


the same depth. 
Since the amount of contamination also increases with field width, 


some of the contamination is attributed to scatter in thisescheme. sinus 


is merely an academic distinction in treatment planning of homogeneous 


phantoms. In treatment planning of inhomogeneous regions, the scatter 


contribution from a region is increased if the inhomogeneity has a high 


electron density and it is decreased if the electron density is low. This 
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reflects the direct proportionality of the Compton cross-section with 
electron density. If the inhomogeneity is in the build-up region, then 
the contamination dose erroneously attributed to phantom generated 
scatter will be similarily modified. In fact, the electron component 
of the contamination dose should be reduced if high electron density 
inhomogeneities are present in the build-up region and increased if low 
electron densities are present, which is opposite to the treatment 


planning for phantom generated scatter. 
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4. CONCLUSIONS 


Conclusions are usually 
consolidated guesses. 
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The Mevatron -6 and Mevatron -20 accelerators produce photon beams 
which were found to be contaminated by electrons and scattered photons. 
The dose measured at the surface of a phantom was due mainly to electron 


contamination. The magnitude of the surface tissue maximum ratio varied 


greatly under the conditions studied between 0.05 (in an open 5 cm x 
wom oe MV t1eld) at. an SP). of 100 cm) to,.0.60 (in-a 50-em xes0rem 
field with a 3.2 cm Lucite accessory in place in a 6 MV beam). The 
other contribution to the surface dose due to photon backscatter was 
usually an order of magnitude lower. 

The electron component of the surface dose was clearly evident. 
The magnitude of the displacement of the TMR curves with increasing 
field width was greatest at the surface. Under the same conditions of 
field width, SPD and accessory configuration, the 6 MV accelerator had 
a greater dose than the 15 MV accelerator at the same relative depth. 

Adding additional thicknesses of accessories increased the surface 


TMR up to a thickness approximately equal to or greater than the electron 


rangepin the accessory (v1/u,*) when increases in the accessory thickness 


did not cause further increases in the surface TMR. 

Electron 'filters'’ reduced the surface dose when an accessory was in 
place by as much as 20% if the atomic number of the ''filter'' was equal 
to lead. The "filter" acted by scattering electrons out of the field. 


A reduction of dose inside the field accompanies an increase in dose 


outside the field. 
The surface dose increased rapidly when the source-to-probe distance 


decreased. At 15 MV the isocentric normalized dose dependence on SPD 


for a 30 cm x 30 cm field size was consistent with the field flattening 


filter being the source of contamination. At 6 MV the isocentric 
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normalized dose did not decrease according to the inverse square 
dependence so the source was not localized. 

The lateral distribution of the surface dose for a field size of 
30 cm X 30 cm at 15 MV was Gaussian. This further suggested that the 
surface contamination was due to electrons which emanate from a localized 
source. The surface dose in a 30 cm x 30 cm field at 6 MV was not 
normally distributed which was consistent with a source or sources of 
electron contamination that are not localized. 

The measurements were always taken with square fields except in 
Section 2.6 where the surface dose of a 15 MV beam was measured for 
rectangular fields. The surface dose was somewhat greater if the upper 
collimator defined the long axis compared to the situation when the 
lower collimator defined the long axis. The surface equivalent square 
field was very different from the equivalent square field at oR ae 

The surface dose has been shown to be a useful measure of the 
amount of electron contamination incident on the phantom. The amount 
of contamination at the surface of a phantom was more accurately assessed 
by determining the amount of dose due to photons backscattering from 
the phantom. 

Contamination electrons produced between the target and accessory 
tray were swept from the field by passing the beam through the poles of 
a permanent magnet positioned at the accessory tray. The surface dose 
with the magnet in place and the photon backscatter dose taken into 
account, was due to electrons produced in air between the magnet and 


phantom. This source of electron contamination was the most important at 


6 MV and for small field sizes at 15 MV. The air produced electron 


contamination increased rapidly as a function of field size at small 


field sizes and increased as a function of SPD. 
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Electron contamination produced between the target and magnet, as 
indicated by surface dose with the photon backscatter removed, increased 
approximately in direct proportion to the cross-sectional area of the 
beam. This contribution to electron contamination decreased as a 
function of SPD. 

The inverse square fall-off in the surface dose produced between 
the target and magnet along with the SPD dependence of the surface dose 
at large fields with no magnet in place suggested that electron contamina- 
tion tended to be produced by mass close to the phantom and this 
material attenuated or scattered electrons produced further from the 
phantom. 

The clinical aspects of an electromagnet to sweep contamination 
electrons from the field were discussed. The requirements for such a 
magnet were contradictory and would be difficult to achieve in practice. 

The penetration of contamination into polystyrene, both inside and 
outside the field, had similar characteristics. The contamination 
inside the field was obtained by equating the increase in the dose as a 
function of field size in a limited phantom (where the field size is 
greater in extent than the phantom) at various depths in the phantom 
with the amount of contamination penetrating to that depth. The con- 
tamination as a function of depth in a phantom outside the field was 
measured by tilting the phantom so that the phantom would not change its 
relative distance from the field boundary as the phantom thickness 
increased. 


Both penetration characteristics illustrated that electrons 


deposit their energy at shallow depths in the phantom and scattered 


photons contributed at deeper depths. The difference between the shape 


of the curves was due to a difference in the relative proportion of 


electrons and scattered photons inside and outside the field. 
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The penetration curve of electrons produced between the target and 
magnet and a penetration curve of electrons set in motion in a thin 
slab at the surface of the phantom were also obtained. Electrons which 
were produced further away from the phantom were more penetrating at 
shallow depths and, therefore, more energetic than those produced close 
to or at the phantom surface. However, the maximum range of the electrons 
was independent of where they were produced. 

The contamination penetration curves normalized to their own 
maxima were approximately independent of the configuration of the beam 
or the condition of measurement. Only the amount of contamination varied 
as a function of field width, SPD and presence of accessories. The 
phantom generated normalized dose at any depth was found by subtracting 
the amount of contamination penetrating to that depth. Since the 
shape of the contamination penetration curve was invariant, the amount 
subtracted at any depth was determined by the amount of electron con- 
tamination at the surface. 

The output factor for the phantom generated normalized dose was 
larger at small field sizes and plateaued to a smaller value at larger 
field sizes compared to the total output factor. The phantom generated 
normalized dose curves were divided by the maximum normalized dose of 
each curve to obtain the phantom generated tissue maximum FALIONCUIVeS 
The phantom generated tissue maximum ratios were independent of field 
width, SPD and presence of accessories. 


The phantom generated dose in the build-up region was computer 


modelled independently of experiment. The calculated results agreed 


with the experimental phantom generated tissue maximum ratios and were 


independent of field width. 
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The computer model showed that electrons set in motion in a thin 
layer produced a dose in the phantom which decreased exponentially. 
Therefore, an analytic model of the phantom generated dose in the 
build-up region could be adopted in which the primary photons and 
electrons set in motion fall exponentially, each with their own 


€haracterasticsattenuation constant. 


The separation of contamination from phantom generated contributions 
could reduce the data storage requirements of small computer treatment 
planning systems. Large computer planning systems could employ the 
results to produce a more accurate model of dose deposition when 


inhomogeneities are present in the build-up region. 


The energy fluence of scattered photons is greater than the energy 
of electron contamination even though the maximum dose of contamination 
electrons is greater than the maximm dose of scattered photons. The 


Shitt. in Shane a is due to contamination electrons. 
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APPENDICES 


All the mathematical sciences are founded 

on relations between physical laws and laws 
of numbers, so that the aim of exact sciences 
is to reduce the problem of nature to the 
determination of quantities by operations 
with numbers. 


James Clerk Maxwell 
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>.1 Appendix i Numerical Values For Constants 

=s5 

es 1.662% 18 Couloumb 
-34 

feet on 625.4, 10 Joule-second 
-34 

f Swi. OD On Ao be Joule-second 
= ia 

m= -3'..116r-< 10 kilogram 

: S 


civ 2.S3ss x 10 meter7secondad 


mee = @.5110 Mev 
7 -15 
i = 2.sls * 10 meter 
22 
hs 6.822 X 18 molecules/mole 


Meéwater) ~ 70 ev 
Wewater) ~ 34 eV~1i10n pair 


The quantity q(T-hvy,2Z) for water (see Equation 1.2.11) 
is listed below as a function of v (see Equation Lea 
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and Compton Attenuation 


Coefficients for Water (9) 


Total 
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>.3 Appendix 3 Total Stopping Power of Electrans 


in Water (5) 


Energy Stopping Energy Stopping 
(MeV) Power (MeV) Power 
(MeU/cm) (MeVU/cm) 

@.108 4.2202 1.02 i SoS 
O.15 3.304 2.88 ped ==) 
8.20 2.858 3.80 Lie Sort 

2 i rashes: 2.050 4.08 barat lide) 
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O35 2.289 6.88 (alge A BSN 
@.40 22190 7.88 2108s 
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@.5@8 Pasay a 9.08 Creeleoa 
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5.4 Appendix 4 Mean Square Mass Angular Scattering 
»~ 
Pawer of Electrons in Water and Lead 


Energy Water i Wo baal oy 
(MeV) (rad“g/cme) (rad@g7cm®) 
2 2 
8.10 dU ba Np ea He) ley a 
1 Se 
8.15 SSS3)'K ire 4,.68°*% "16 
1 2 
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Q@.68 too ST4S7xK 72¢ 
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@.8G 4.75 an S5 24912 
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1.08 Sagse he! Ca ee ae 
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li. SO ineSie 1.2.28 — xetO 
2.09 5 TE aT S.c2 
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3.88 6.88 X 18 4.63 
=4 
4.00 Siever a eae aah, 
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5.28 2.55 A Le Ade, 
cil 
6.08 TOT Uk 28 1.46 
25) cb 
8.08 LiL Sea ee 8.87 xX 12 
mae <5! 
18.02 7.61 X 18 6.88 X 128 
ee) =e 
15.98 3.67 X 18 2.91 X 18 
=e =f 
28.98 rie i ou ee 1273" “10 
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S12 rAppeNnd:x —5> Listing and Documentation For 


“Se primary ahoTr « 
program seprimary 


ROROROROROROROROROROROROROROIOROICORORONOIOINOROIOROROICOIROIC HONOR IC IOROROROIORORORORORGORORORO RG RO ROR ROR IG CI RCI ROR 2 KORCOK 


* sd 
«x This program subtracts the contamination normalized dose from x 
*x the total normalized dose to yield the phantom generated norm- * 
* alized dose. The amount of contamination arriving at the phan- * 
*x tom is the proportional to the difference between the total * 
* and photon backscatter normalized doses. The phantom generated * 
*x TMR curve is calculated after a search has found the maximum * 
* normalized dose. The error in the normalized dose and TMR * 
* values are calculated. * 
* * 
* OB 


OOO IOK IGOR IGIOIOIOI I IOIOIOIOIOIOIOIOR OR III IOIIOK III IOIOIOI OI IOI IO FOR IOI ROK IOI ICR KK ICR 


integer i-n 

real deptnh(2@),ctmr(20),w,nd(28),bsnd 
real ernd-ebsnd,cmax,pnd(2@),epnd(2ea) 
real recmax,max,remax,ptmr(2@),eptmr (20) 


YOIOIOIOIOROIOIOI IOI IOI OIOFOIOIOIOIOIIOIOIOICI IO III IOIOIDIOIOIOIOIOIOIOIOI IOI IOIOIOI IORI OIOIOIOR OCI OK ROOK 


« 

* deptn(28) is the depth in the phantom. 

*x ctmr(28) is the contamination tissue maximum ratio,TMR, as a 
* function of depth. 

* wis the field width. 

* nd(2@) is the total normalized dose as a function of depth. 

* bsnd is the backscatter normalized dose. 

x ernd is the error in the normalized dose. 

* ebsnd is error in the backscatter normalized dose. 

*x Cmax is the value of maximum contamination normalized dose 

*x pnd(28) is the phantom generated normalized dose as a function 
* of depth. 
* 
x 
« 
« 
* 
« 
‘« 
« 
« 
x 
* 


% NK HH 


epnd(2@) is the error in the phantom generated normalized 
dose as a function of distance. 

recmax is the relative error in Cmax. 

max is the maximum phantom generated normalized dose. 

remax is the relative error in Max. 

ptmr(28) is phantom generated normalized dose as a function 
of depth. 

eptmr(20) is the error in the phantom generated normalized 

dose as a function of depth. 
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NOS CICIEIE ISIE OGRE IOIOICIOI IOI ICICI III IC IR JOR IOI IRI NOR SIA ROE SOE AOS SS 


Bpentuniteirstatus= Ola". file=‘contemio-dar 
Bhemtunit=2, status: "Old™,file="ndisnole- ca... 
opentunitzgrscatose ned ss Mest prumsvedal © 


“ce RTE ee RE AE ERE RHEE EES SR EE EAEL ELS SETS 


contam.dat stores contamination TMR values as 4 function of 
Gepth. 

ndiSnoi6.dat stores the total normalized dose curve as a func 
tion of depth and the value of the photon back- 
scatter normalized dose for a particular config- 
uration of field size and SPD at 15 MU. 

is the output file that lists the phantom generated 

normalized dose and TNR values along with their 


associated errors. 


primiS.dat 
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48 
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188 
S88 
6a0 
780 
1008 
1188 
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read (1,188) n 
Go SS vetsn 
read (1,508) depthn(id.,ctmr ci) 
read (1,600) ectmr 
read (2,708) w 
do 18 We Ur 9} 
read (2,188@8) ndci)d 
read (2,1180) ernd,bsnd.,ebsnda 


ROROIOROIOIIOI IORI OIOROIOROIO HORI IOIOIOIOIORORO ORR OROIORORRRORORCORORCOIC IIRC ROR ICRC CIC ICR ICR RO ICICI KOK OK 


* 
x 
x 
* 
x 
« 
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The following statements calculates the maximum amount of con- 


tamination and the amount of contamination occurring as a 
function of depth. The contamination is removed from the total 
normalized dose. The relevant errors are calculated as well. 


SK OKOROFOROIOROKOIORORO ROR ROROROROROFOR OR OK OR OIOIOK IOI OIOIOIOIOROIOKOIOIOI KR OIOIOK ROK IOK ROR OKI IKOK I IOROROK OOK 


Cmax=(nd(1)-bsnd)7ctmr (1) 
recmax=(ernd+ebsnd)/7(nd(1)-bsnd)+ectmr7ctmr (1) 
pnd(id=bsnd 
epnd(1)zsebsnd 
do 28 i=e-n 

pnd(id=endlid-cmaxxctmr (i) 


epnd(i)dsernd+(recmax+ectmr7ctmr (i) )*xcmaxxctmr (i) 


* 
* 
* 
* 
»* 
* 
* 


OROROIOIORORORORO ROR UI OIOIOIOIOIOIOIOIROROR ROI ICIOICIOROROIOROI UI IOIO IOI IO IOROI I ROIOIORO IORI OIOICIORO IC IKK KICK > 


*“ 
x 
x 
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The maximum phantom generated normalized dose is found and the 


phantom generated TMR is calculated along with the associated 


errors. These quatities are then tabulated. 


SOOO OIOOIOIOIOIOK IOI OOIOIOI IOI IR OR I III IORI IIR IOI IOIOIOI IOI IOI ROK IOI ROK IOI RICCI ROK OOK KOK 


max=pnda(i) 
do 3@ isi-n 
if(pndCi)d .ge. max) then 
max=pnd(i) 
remaxzepnd(i)/pnd (i) 
end if 
continue 
do 48 izi-,n 
ptmr (i d=spnd(id7max 
eptmr(id=(remaxtepnd(id7pnd(id)*ptmr (i? 
write (3,120@) ’Field Size =’»w,’cm X’,w,’Cm’ 
write (3,-1308) ’Photon Backscatter N.D. =’,0Snd 


write (3,*)’ ” 


write (3.1580) "Depth’,’ND’,’Error in NDR poh Or 


write (3.1608) ’(cm)’ 
do 5@ iz=1-n 


in TMR’ 


write(3,1708) depth(id-pndtide,epndtid-ptmr(id,eptmr (id 


format(’ ’°,i2) 

format(’ °f4.2:2x/f4.2) 
format(’ ’,f4.2) 

format(’ ’,f4.1) 

format(’ °,f5.3) 

format(’ °,3(f5.3/2x)) 

format(’ »,a12,#4.1,a4/f4.1-a2) 
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format(’ 
format (’ 
format(’ 
format(’ 


close (i) 
close(2) 
close (3) 


stop 
end 
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5.6 appendix 6 Listing and Documentation For 


SBuauUp Ss. Ors 
program buildup3 


i OROCIOROIOIOROROIOIK IORI ICRORORO ROI OIRO OKO ICRC ROR ROROROROROIOIROIOIOIO CORO ROR IORI RO ICROIOIOROIIOIOHORROKOROR OK ROK ICOHOK 
c * Electron Monte-Carlo program to determine the electron dose 
c * in a semi-infinite water phantom for polyenergetic electrons using * 
c * the Rossi Gaussian lateral scattering approximation. 
ic ORO ROK OK IOI OROROKOR ROR OROR ROK ORO KOR IK IIR RIOR OROROROIOK IOIROR ROK RRO IORI IORI ORO KORO ROR ROI ROR OR OROROK KOROK KC 2K 2K 
[a ROROROROIOROROIOROIOROIUIOI ORCI ROROIOIOIOIOIOIOIOROROIOROIOROIOIOIUIOROIOIO ROI OROROI ROR IOI RO IOIOROIOIOIOIC OKI OKOK 
(S *x*x Complete variable documentation appears in file Buildup3.doc * 
c eOk OOK ROK OR FOROROROIOIO ROR ROR ROK ROK ROK RRO RRO UR KK ROK ROR ROK IOROIOIOROIOROIO RIOR ROR OROIOROIOROOIORCOK ROOK 2K 
integer l,il-nerg 
real r,zs-,Cal,ssal,cbeta,sbeta 
real stop(179),tna(i79),ep,cay 
real eo-e,dose(260),pl,gaus(1i8@),Qg 
real z-max,ang,en(36),spectr (9G, 36) 
integer i»j»m-p»kay 
real xnrg(36),mu(36),Sig(36),1l0cal,que(10) 
common “blkivxnrg-en-mu,sig»local,spectr,que 
common /b1lk2/zs,z-ep,dose 
c FOIOIOIOOIOIOIOIIOIOIOIOIOIOIOIOOIOIOIOIOIOIOIOI IOI IOI OR OIGIOI IOI OIOIOIOI IOIOIOR IOI OIOICIOIOR IOICICI IORI OROIOK ACK OK KF: 
c * Gaus.dat stores a table of x’s as a function of erf(x). * 
c x Quadint3.dat stores values of interpolated scattering and stopping * 
c * power data in step sizes of @.1 Mev. * 
c «x Dose.dat stores normalized absorbed dose as a function of depth * 
c * in millimeters. * 
c * PhotoniS.dat stores the normalized number of photons~bin-total * 
cS * absorption and Compton absorption coefficients as a * 
c * function of bin energy. At the end of this file 15 * 
c * a coefficient array which is used in Stearn’s * 
é * determination of a pair production electron trajectory.™ 
¢ NORCROIOICIICIOIOIOIOICIOIOI IOI II IIEIOIOI GIGI SIOIOIOIOIDIOI OI I OIOI IOI OI IOIGIOIOIOIOIOIOI GIGI OIOR IOI IOI IORI IRR OK 
open(unit=i,status=’o0ld’-file=’gaus.dat’) 
open(unit=2-status=’old’-file=’quadint3.dat”) 
open(unit=4,status=’new’,file=’dose.dat’) 
open(unitz=i1,status=’old’,file=’photon6.dat’) 
c FOIOCIGIOIOI IOI IOI IOI DIGI DIO IOI IO IOIOIOIOI IDI IOI DIO IO IOIIOII GIO ICI IOI II AIOE 
c x The following statements read the look-up tables and initializes. * 
c SAA AEA A CR RHE NO I OR AOR IO HORAK I AR IE ES AMA BI EAS SS 
do z t=i,i73 
read(2,42@) tha (l),stop(l) 
2 continue 
do 10 i=1,36 
read(11i,340) xnrgcidsencids»mucid,sigli? 
18 enci)dzenci)*1080.8 
read(11,350) (que(i),i=i,1@) 
do 30 1=1,188 
38 read(1,422) gaus(1) 
1122650819 
cay=8.98 
Cc ee Sa A RR A HORT HRA AK AARC HN BIER OAD RD ee Nae i 
c «x Subroutine Impulse is called to determine the electron-set-in-motion e 


c * number spectrum, 


’spectr(t,thneta) ’-whnere "t’ is tne electron energy 
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x and ’theta’ is the original direction with respect to the photon * 
x Girection. The program iS repeated at all electron energies and * 
* angles for all histories as specified by ’spectr(t»,theta)d’. The * 
* following statements count histories and initializes the electron * 
* energy and angle. ne 
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do 1280 kay=18,188 
call impulse 
dose(kay)=dose(kay)+local 
do L2@j=1,.75 
eo=floatjCjo75.8+.1 
GOwtceO p-isi6 
ang=floatj(p)*5.0-2.5 
do 120 l=1,jnint(spectr(j,p)) 
e=e0o 
zs-floatj(kay)/10.8-.85 
cal=cos(ang/57.296) 
sal=sin(ang/’57.296) 
r=6.2832*ran(il) 
cbheta=cos(r) 
sbetaz=sin(r) 


OIOROIOIOIOOIOIOIOIOROROIOROIOIOICIOIOIOIIIOIOROIIOROICIOIOIOIO IORI I RCIOIC IC ICICI IOI I ICICI ICI I IO iC OK 
x ’Tha(nerg)’ and ’stop(nerg)’ are arrays that contain the 

“x mass angular scattering and mass stopping powers respec~- 

x tively. ’E’ is the kinetic energy of the electron and 

* ’ep’ is the amount deposited in the patune Vengt hs spl a. 

x ’Gaus(rgti) contains equally probable scattering angles 

* for a GausSian distribution with mean=8 and standard dev- 
* jationz1. When this distribution is multiplied by the 

* mean scattering angle the result is ’g’,»the zenith scat- 

* tering angle. The final value of ’r’ in each scattering 

«x event represents the scattering azimuth angle. The Sub- 

* routine Geom determines the angles in the coordinate 

* system fixed with respect to the phantom. 

SOI GICIOIOIOIOIOIOIGICICICICIDI IOI IOIDIOI OI IDICICIOIO OOO CII IO IO IOI TT IA AOA AOE 
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do 118 whilete .gt. 8.2) 
ep=(1i-cay)*e 
nerg=ex*10.8-ep*5.@ 
plzsepvstop(nerg) 
e-e-ep 
gzsqrt(tha(nerg)*pl )*gaus( jifix(100.@erancil) +1) 
r=6.2832"ran(il) 
call geom(r-g-cal,sal,cbeta,sbeta) 
z-epl*cal 
zs=zSt+zZ 
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* The following statement checks for the edge of the phantom. * 
a RAIA AAAI AAAI AEE T TT TEESE EEL SEES 


jp Sdzes *s1t. 8.8 ) go to aS 


4 PARRA AAA SAA AINA IAA EAA II OAT AIS I IIIA I AAR ID 
«x The energy ioss is assigned to the cumulative dose * 
x received by the phantom using the Subroutine Assign. * 


iol ASA ARAL ASIA TEATS A IAL TI SIDA IA IAS ESS 


call assign 
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continue 
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* At this point the primary electron has an energy less than 


*x* 86.2 MeV. The remainder of the energy 
* at last depth interval. 


ROROIOROIOROROROIOROIOIOIOROROIOIOIIOROIOIOIOROROIOR RIOR IORI IORI IOIOROROIORCOROROIOIOK ROI ICICI IO IOROROIOIOIC IC OK KC 


call assign 
continue 
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is assigned to the dose * 


* The following statements normalize the dose to the maximum 
*x dose encountered and write the normalized dose. 
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max=8.8 
do i22 i-1,180 

if (dose(i)d .gt. max) maxzsdose(i) 
continue 
do L2Saui=17 168 

write(4,250) i,dosecid/max 
continue 


format(’ °,i13+2x,f6.2) 

format(’ »,€#5.2,2x,f5.322x,26F5.2-2x)?) 
format(’ °,18(f4.2,2x)) 

format(’ °,e18.4) 

format(ei@.4) 


close(1i) 
close(2) 
close(3) 
close(4) 
close(il) 


stop 
end 


subroutine geom(r,g-cal»sal,cbeta,sbeta) 


real r,g-cal,sal,cbeta,sbeta 


real ss,cc,oldcal,oldsal,oldcbeta,oldsbeta 


YOR OIC OIDIOIOOIOIOIOIDIOIGIOOIOIOI GIO III IOI IORI OR IOI IR IOI III I IEA 
* ’Sal’ and ’cal’ are the sine and cosine respectively of * 
* the zenith angle in the phantom system. "Sbeta’ and * 


* cbeta are the sine and cosine respec 


* muth angle in the phantom system. 


tively of the azi- * 
«K 
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oldsal=sal 

oldcal=cal 

oldcbeta=cbeta 

oldsbeta=sbeta 

calzcal*cos(g)+salxsin(g)*cos(r) 

if Cabs(cal) .ge. 4.8) then 
calzii.@*sign(1.@,cal) 


sal-9.9 
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ORR OR ROK OIRO ROI OI RIOR 
* The following statements determine the number of photons that 


mx 
* 


WOOK 


cheta=cos(r) 

sbetazssin(r?} 

go to 488 
end if 
Salzsqrt(1.6-cal«xx«2) 
if (oldsal .eq. 8.8) then 

cbheta=scos(r) 

sbetaz=sin(r) 

go to 488 
end if 
ssz=sin(g)*sin(r)7sal 
ec=(cos(g)-oldcal*cal)voldsal/’sal 
if Cabs(ss) .ge. 1.0) ss=1.8*sign(1i.@8-ss) 
if Cabstec) .ge. 1.8) ccz1.0xsign(1.B,cc) 
cheta=cc*oldcbeta-ssxoldsbeta 
sbeta=ssxoldcbeta+ccxoldsbeta 
if Cabs(cbeta) .ge. 1.8) cbheta=1.8*sign(1i.@-cbeta) 
if (abs(sbeta) .ge. 1.0) sbeta=1.0*sign(1.8,sbeta) 
continue 


return 
end 


subroutine impulse 


real en(36),spectr(9@, 36) 

real xnrg(36),mu(36),sig(36), local 

integer i»j 

real nteract(36),ncomp(36)-npp(36) 

real nrgmax,sum,que(1@),q 

real kinrg(20)-nrgdif,a,b-c,dsigdt(e8),arg 

real theta(2Q@),v,vee-,gi-g2-g(20),nbin,»angle-prob 
integer binang-binum,vint 


common /blki~xnrg-en,-mu,sig, local,spectr,que 


ROR IOI ICICI IOI IOI III IOI OR IOI ICIOR IORI IOFIOK ROI IOIOIO IORI NOK IC OK ACK 


* 


interact in the phantom slice and the number due to the Compton x 


and pair production interactions in each photon energy bin. 
YOO IOOIOIOIOIOIOR IOI IOIIOIIDICIGIDIOI IOI UII IOI IOI IOI ICICI OI IOI OICI ICRI ROK IORI ICICI I AOKI 


do 5 #17906 

do 5S Hees 
spectr(i,j)=8.8 

local=8.8 

do 4@ i=1,36 
qzen(i)d*exp(-muC(id*. 681) 
nteract(idsen tii-q 
encid=q 
ncomp(id=sig(id7muCi)d*nteract(i) 
npp id=sex(nteract(id-ncomp i> 
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« The following statements generate the electrons-set-in-motion 


*x due to the Compton interaction. 
YO IOI OOICI IOI IOI IOI IOI IOI III I IO ROOK 
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48 


sum=8.@8 
do 18 j=1,28 
kinrg(jafFloatj(jiI-.S)*nrgmax72a.8 
nrgdif=xnrglid-kinrgij) 
az(.S5iixkinergljd7Cxnargli) )«x2)xx2 
b=ex(nrgdif“xnrg(i))**2 
c=((kinrg(ji-.511)**2-. 261) *nrgdif7(xnrglid)**3 
dsigdt(j)=Catbt+c)dvnrgdifxxe 
Ssum=sum+dsigdt(j) 
arg=(nrgmax-kinrg(j) *C2exnrglid7.Siil+1d~kinrglj) 
theta(j)257.296*atan(saqrt(argd7(itxnrglid7.Siid?) 
continue 
do t2 j=i.28 
dsigdt(jdsdsigdt(j)d/“sum 
Nnbin=dsigdt(j)*ncomp li) 
binmang=jifix(thetaCj)*.2)+1 
NP CIR Rina CHD oiiea 67) AoCtd 
local=local+nbin*kinrg(j) 
else 
binumz=jifix(kKinrg(j)*5.8) 
spectr (binum,binang)=spectr (binum,binang)+nbin 
end if 
continue 
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x The following statements generates the electrons-set-in-motion * 


*x due to the pair production interaction. x 
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nrgmax=xnrglid-1.822 
sum=8.9 
if (npp(i)d .gt. @) then 
ele Bley esloae 
kinrg(j)s(float(j)-.S)*nrgmax728.8 
ve(kKinrg(j)+.511)7xnrglid 
veez=1i-v 
gizvxxZ+veexxZ2+2. Bxumvee 73.8 
g2=alog(2.@«xnrg(i)xvavees.511)—-.5 
g(jd=gixge 
sum=sumt+g Cj) 
vintsjifixCv*l18.8>+1 
thetal j)=57.296*que(vint walog(xnrg(id7.Sii)*.S1i7xnrgis> 
continue 
do 17 jei-20 
gt j=gCji7sum 
nbinz=g(j)*nppli)d 
binumzjif ix (kinrg(j)*5.8) 
do 17 binang=1-,36 
anglezfloatj(binang)*.@8727-. 84363 
drob=.c69es7thetakj exp (= Canglecene vac) iE tegeee 
Spectr (binum, binang)=spectr<Pinum:Gamaqg) “nbi apr e® 
continue 
end if 
continue 
return 
end 


subroutine assign 


integer irj»k»kp,»dk 
real z2s,z,dosefrac,ep, dose (200) 
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common /blk2/zs,z-ep,dose 


ROR ROR OROK ROR OFOROIGIOFOR ROR RO FOR KORO IO K IOIOIC OK 2 9K 


YOROROROIOROR ROR ORO IOROROROROIOI OIRO IOIOR ORI OIOR OR IOKOK 


x The fractions of dose to be assigned to all elements traversed * 


* are calculated. 


* 


_YOIGIOICICIOIOIOIIOIOIOIOIOI I IDIDIOIOIOIOIOI IIOIIOIOIDIOIOIOIGI I IOIGIOIIDIOIDIOIOIOI I III IOI OI IK FON CI I IK KF 


k=186.0*zs 
kp=10.0*(zs-2) 
TG Keen) tO unt Mem 
do S ist,kp 
dose (kp+i-i )=dose(kp+i-i) 
dose (kp+i)edose(kpt+id+epx( ( 
go to 55 
end if 
dk=k-kp 
j=2 
if (jiabs(dk) .gt. 1) then 
jzjiabs(dk)-1 
do 18 izi>j 
dose(kptitsign(i-dk))=da0s 
end if 
if (dk) 48,45-58 


-ep/’z/18.8 


z-zs)+floatj(kp)718.80)72 


et(kpti+sign(i-,dk})+abs(ep718.072) 


dosefracz(zs-floatj(k+1)718.0)/7Z 


dose (k+1)=dose(k+i)+epxdose 
dose (kp+i1)=sdose(kp+i)+ep* (i 
go to 55 


dose(k+1i)=dose(k+id+ep 
go. toxs> 


frac 
-floatj(j)718.87z-dosefrac) 


dosefracz=(zs-floatj(k)710.0)72 


dose (k+1)=dose(k+i)+ep*xdose 
dose (kpti)zdose(kptidrep*(i 


continue 
return 
end 


frac 
-floatj(j)718.87z-dosefrac) 
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Variable Documentation For Buildup3.For 


Main Program 


Integer 
i»je>l»p,» Dummy indices. 
il, Random number seed. 
kay», Slab number (1 is the top slab) in which the photon beam 
interacts . 


nerg, Energy parameter specifying locations in look-up 
tables. 


Real 
tna(i79), A look-up array storing the scattering power 
of water in radian**e/(cm*x2/Qg). 
stop(179),» A look-up array storing the stopping power of 
water in MeV7(cm«x*27Qg). 
xnrg(36), An array storing the photon bin energy. 
en(36), The normalized number (in % ) of photons in the bin. 
mu(36), The total attenuation coefficient for water. 
$ig(36), The total Compton absorption coefficient for water. 
que(1@), Stearn’s coefficients used in the determination 
of pair production electron trajectory. 
gaus(1@@), Values of the inverse of erf(x). 
cay» Represents the fractional amount of energy retained by 
the electron after each step. 
dose(200), The accumulated dose received by the slab. 
local,» Dose deposited in the slab by electrons set in motion 
with an energy < @.2 Mev. 
eo,» The initial energy of the electrons. 
e, The energy of the electron. 
ang» The zenith angle (degrees) at which the electrons are 
set in motion. 
spectr(98,36), Generated by Subroutine Impulse,» this array 
contains the number of electrons set in motion 
as a function of energy (9@ bins) and zenith 
angles (36 bins). 
zs, The depth of penetration into the phantom. 
cal, The cosine of the zenith angle in the phantom coordinate 


system. 
sal» The sine of the zenith angle in the phantom coordinate 


system. 
r, The azimuth scattering angle (radian) in the electron frame. 
cbeta» The cosine of the azimuth angle in the phantom coor- 
dinate system. 
sbeta,» The sine of the azimuth angle in the phantom coordinate 
system. | 
ep,» The energy absorbed by the phantom in a step. 
pl,» The path length (g7cmax2) travelled by the electron during 
the step. ; 
g, The zenith scattering angle (radian) in the electron frame. 
z» The change in the penetration during the step. 
max» The maximum dose received by the phantom. 


Subroutine Geom (Variables listed only where different from the 
Main Program. ) 


Real 


oldsal» The value of sal at the end of the previous step. 
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olscal» The value of cal at the end of the previous step. 

oldcbeta, The value of cbheta at the end of the previous 
step. 

oldsbeta, The value of sbeta at the end of the previous 
step. 

SS»-cc, Parameters used to simplify geometric calculations. 


Subroutine Impulse (Variables listed only where different from 
the Main Program.) 


Integer 
binang, Parameter specifying angle bin. 
binum, Parameter specifying energy bin. ; 
vint, Parameter used to specify Stearn’s coefficient. 


Real 
q- The number of photons transmitted through the slab. 
nteract(36), The total number of photons in the bin interact-— 
ing in the slab. 
ncoomp(36), The number of electrons set in motion by photons 
in the bin that interact in the slab due to the 
Compton effect. 
npp(36), The number of electrons or positrons set in motion by 
photons in the bin interacting in the slab due to pair 
production. 
nrgmax, The energy that is available to transfer to charged 
particles. 
sum,» Sum of the cross-section parameters. Used for normalization. 
kinrg(2@), Electron kinetic energy. 
nrgdif, Difference between photon and electron kinetic energy. 
a»b,rc,» Terms used to calculate the parameter proportional to the 
Klein-Nishina differential cross-section with respect 
to electron energy. 
dsigdt(28), A parameter proportional to the Klein-Nishina Ont 
ferential cross-section. 
arg,» The argument used to simplify the Calculation of recoil 
electron zenith angle for the Compton effect. 
theta(2@), The charged particle recoil zenith angle. 
nbin- Product of normalized cross-section and number of charged 
particles set in motion. 
vu, Ratio of total electron energy to photon energy. 
vee, Ratio of difference between total electron energy and 
photon energy compared to photon energy. . 
gi+g2,» Terms used to calculate the parameter proportional to 
the pair production differential cross-section with 
respect to charged particle energy. 
g(28), A parameter proportional to the pair production differ- 
ential cross-section. 


(Variables listed only where different from the 
Main Program) 


Subroutine Assign 


Integer ; . 
k» The slab number (counting from the top down; the first slab 


is labelled @) the electron is in at the end of the step. 
kp» The slab number the electron was in at the end of the 


previous step. 
dk» The number of slab boundaries crossed. 
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Real 
“ gGosefrac., 


The fraction of dose deposited to Slabs across bound- 
aries. 
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>. APPendix 7 Listing and Documentation For 
*‘Normdist.Ffor ’ 


program normdist 


FOCIOIOIOIOIOIOIOIOIOIOIOICICICIOIOIOIOIOIOIOIOIOICIOIOOIOIO I IOI ICI IOIIOR IOIOIOI IOI IOIOIOI IO IOI IOI K III ORK IOIOK I ROR Keak 
* This program calculates erf(x) using a modified approximation of the x 
* formalae listed in Dwight’s Tables of Integrals And Other * 
x Mathematical Data page 136. * 
ZCIOIOCIOIOIIOIOIOIOIOIOIOIOIOIOIOIOI IOI IOI IOI IOIOI IOIOIOIOIOIOIO IOI I IOI OIOIOIOIOIOR IR OK IOI IOI IOI IR IK II I I I tO ok 


integerx*x4 factor : 
do 18 tea=8.8001,2.6,8.0881 
nike 
sum = 8.8 
test = 1.0 
do 9 while (test.ge.@.8001) 
term=(-1)**x(n-1)x(teaxxn—-1)/nxfactor(n—-1) 
sum=sum+term 
testz=term/sum 
n=n+i1 
Sum=sum/1.25331 
print 11,tea,sum 
format (f8.5,10x,f8.5) 
stop 
end 


function factor (1) 


ORR OIOKOROROROR ORO ORK RIOR ROKR IO OIRO OK OIOIOIOIO IORI IORI ORI ROIOROROR ROI IO ROR OK ROR ROKR OK ICRC KK OK OK ROKK OK 


« The subroutine Factor calculates the factorial of an integer. * 
ROROIOROROROIOROROROIOROR OIRO IO ROROIOOROROIUI ICIOIOROIOIO IOI UIOROROIOIOIOIOIOIOIOROIOROI ICR ICHOIC IIRC ROROIOOIOIOIC IO IORI KICK IC Ic 


integer factor 

factor=1 

Uf) Ci. Newt) Ce CUR 

dots) jiscne 
factor=factor*j 

return 

end 
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5.8 Appendix 8 Listing and Documentation For 


*‘Addose.For’ 


program addose 


c OIOIOROROROIOROIOIOROIOIIOIOROROIOROROIOIOIUIOIOIOR OF OIOIROR IOI NIORGROI ORO IOI ICIOI IO RORCROI IC IO ORCIOI IC ORONO HC IK KCK 
c * Addose.for calculates the dose as a function of depth in a * 
c x phantom from electrons set in motion in @.i1cm slabs. The photon * 
c * energy fluence is attenuated exponentially. The result is the * 
c * broad beam TMR build-up values at @.icm intervals of depth. * 
[my YOOOIOROIOIOIOIOIOIOIOIOFOIOI GIO OI OIOIOR ICICI IOI OIOIOIOIOI IOI IOIOI I I IOIOIOI IOI IOI OIC IOI IOK NOI IK AICI ACF 
integer irjrk 
real dose(i1@Q),max,adose(18@) 
open(unit=1,status=’ola’,file=’dose.dat”) 
open(unit=2,-status=’new’-file=’addose.dat’) 
do 18 i=1,1088 
168 read(i1,108) k-dose(k) 
¢ NOIQIOICICIOIOIDIOIOIOIOIOIOIOIOIDIOIOIOIOIOIOIDIGIOIOIOIOIIOIUIDIDIOIIOIOIOIIOIOIOIDIOI IOI IOI I IORI I IOI II IOI KOK 
c x The following statements attenuates the photon energy fluence * 
c *x at the deptn at which the electrons are being set in motion. * 
c « The attenuation coefficient is 8.829/”cm which is the average * 
c x energy absoption coefficent at 6 MV. % 
c SOOO IOIOIOIOIGOO IO IOIOIOI IO IOIOIDIOIIOR I OIOIIIOIOIOIGI I IDIOIIIOIOI IOI IOI IOI III I ROR I FOF 


do 45 i=1,188 
do 48 jz=i-i 
adosei )sadose( i)+dose tj) xexp (=, 0290 Cf loatyCity=1)ent=. 057) 

42 continue 
45 prints, ’iz’,is’adoseti)=’,adose(i) 

max=8.0 

do 50 k=1,188 
58 if Cadose(k) .gt. max) maxzadose(k) 

do 60 k=1,188 
62 write(2,188)k,adose(k)7max 


188 format(’ °,i3,2x,-f6.2) 


close(1) 
close(2) 
stop 

end 
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The 15-MV photon beam of a linear ac- 
celerator (Siemens Mevatron 20) was 
studied for electron and scattered photon 
contamination. The surface dose, attrib- 
utable almost entirely to contamination 
electrons, has a Gaussian lateral distribu- 
tion, a linear dependence on field width 
for square fields, and an inverse square 
dependence on distance from the bottom 
of the fixed head assembly. This geomet- 
rical dependence is consistent with the 
proposal that the field flattening filter is 
the main source of electron contamina- 
tion when accessories are absent. A tis- 
sue-maximum-ratio curve in the build-up 
region for the electron and photon con- 
tamination was produced utilizing the 
linearity of dose with respect to field 
width. The derived contamination curve 
inside was similar to the measured build- 
up curve outside the field. The primary 
photon component, obtained by subtract- 
ing the contaminant contribution, 
showed no dependence on field size, 
source-to-probe distance, or presence of 
accessories. 


Index terms: Electrons e Linear accelerator ¢ 
Radiations, measurement « Therapeutic Radiology, 
instrumentation 


Radiology 144: 403-409, July 1982 


1 From the Physics Department, University of Al- 
berta and Cross Cancer Institute, Edmonton, Alberta, 
Canada. Received July 24, 1981; accepted and revision 
requested Nov. 13; revision received Feb. 16,1982. jr 


Contamination of a 15-MV Photon 
Beam by Electrons and Scattered 
Photons! 


ya understanding of the role of electron and secondary photon 
contamination in megavoltage therapy machines is required 
to characterize dosimetry in the build-up region. Clinically, high- 
energy x-ray build-up produces a skin-sparing effect and contami- 
nation tends to reduce this effect. 

Electron contamination has.been identified as the major source of 
dose a‘ the surface of a phantom for ®°Co (1, 2) and a wide range of 
accelerator energies (1-7). Most of these studies have employed 
thin-window parallel-plate chambers to measure the surface dose. 
Padikal and Deye, and Biggs and Ling, have measured the electron 
contamination directly by sweeping the electrons magnetically from 
the beam (3, 4). 

There has been widespread disagreement concerning the origin 
of the contamination electrons. The collimator jaws (1-3, 5), field 
flattening filter (1, 2), and intervening air between the source and 
detector (2, 8) have been cited as possible sources. It is generally 
agreed that the presence of assessories in the field of the beam will 
increase the magnitude of electron contamination (1,5, 6, 8-10). Bagne 
as well as Nilsson and Brahme, have derived expressions for dose in 
the build-up region by assuming that the surface dose is due to the 
presence of contamination electrons (8, 11). 

Although electrons provide most of the contamination at the sur- 
face, contamination should be defined to include secondary photons 
generated outside the phantom. Throughout this paper the contam- 
ination component of the beam is defined as electrons and scattered 
photons produced by interactions of the primary photon beam with 
material outside the phantom. The primary photon beams are photons 
emerging through the collimator which have been produced by 
bremsstrahlung interactions in the target and modified by the fixed 
head assembly (Fig. 1). 


I. OPEN-FIELD BUILD-UP CURVES 


Materials and Methods 


Build-up curves were measured for a 15-MV photon beam froma 
Siemens Mevatron 20 linear accelerator. A schematic diagram of the 
beam defining head is shown in Figure 1. All measurements were 
made at constant source-to-probe distances (SPD) using 0.16-cm 
(4e-in.) thick and 0.64-cm (1/4-in.) thick square polystyrene slabs as 
sllustrated in Figure 2. The phantom slabs used in the determination 
of the open-field build-up curves each had a cross-sectional area of 


625 cm?. 
The detector used was a Capintec 192A electrometer with a PS-033 
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Figure 1 
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Diagram of the head assembly. 
Figure 2 
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CENTRAL AXIS armor 


BUILD-UP PHANTOM SLABS 


PARALLEL PLATE IONIZATION CHAMBER 2 
i EAE. PHANTOM 


BACKING 


ee 


re) “yO em SCALE 


Diagram of the experimental arrangement. 
When determining the contamination depth 
dose, the phantom slab cross-sectional area 
was 100 cm2 and the slabs were always within 
the field boundary. The phantom slab cross- 
sectional area forall other measurements was 
625 cm?. 


thin-window parallel-plate ionization 
chamber with an effective volume of 
0.5 ml. The entrance window was alu- 
minized polyester film 0.5 mg/cm 
thick. The tissue maximum ratio (TMR) 
is defined to be the dose at some depth, 
d, in the phantom normalized to the 
maximum dose occurring at a depth 
dma, for the same field area A and en- 
ergy E. 
Dose (d,A,E) 


TMR = [1] 
Dose (dwawAcE) 
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Throughout this paper the field size is 
referred to as the lateral field dimen- 
sions defined at 106 cm. Only square 
fields were used. 


Results 


The central axis build-up curves for 
various field sizes are illustrated in 
Figure 3. These build-up curves dem- 
onstrate a dependence on field size. 
The TMRs for all depths below dyax 
increase with increasing field size. The 
greatest relative difference between 
these curves occurs at or near the sur- 
face. 


Discussion 


The spread in the curves is greater 
than can be expected from lateral scat- 
tering of photons and electrons gen- 
erated in the phantom. Indeed, in the 
absence of backscatter or contamina- 
tion, the TMR at the surface for all field 
sizes should approach zero. Thé field 
size dependence and the finite surface 
dose can be attributed to contamination 
since backscatter is negligible at the 
energy investigated. 


Il. ACCESSORY BUILD-UP 
CURVES 


Materials and Methods 


Various thicknesses of Lucite (den- 
sity 1.18 g/cm) alone or Lucite and 
lead (density 11.4 g/cm*) slabs were 
placed in the beam to determine their 
effect on the build-up curves. The dis- 
tal surface of the slabs was placed at the 
accessory tray holder position 56 cm 
from the source. In all cases a Lucite 
surface was distal to the source. The 
build-up curves were measured as in 


Parl 


Results 


Figures 4 and 5 show the build-up 
atan SPD of 100 cm and 75 cm, respec- 
tively. The field size in both cases is 30 
cm *.30..cm. When accessories are 
placed in the field at the tray holder 
position, a larger TMR for all depths 
less than dyax 18 observed. There is also 
a shift in the point of dmar to shallower 
depths. This 1s more pronounced at 
SPD = 75 cm (Fig. 5). There is virtually 
no difference in the build-up curves 
between a 3.2-cm (13/,-ins)s Luveiter ac 
cessory and a 0.30-cm lead slab on top 
of a 0.64-cm lucite accessory. 

At smaller SPDs, there is a greater 
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source-probe distance for various field sizes. 


Figure 4 
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Effect of accessories on the build-up curves for 
a field size of 30 cm X 30 cm. Source-probe 
distance = 75cm 


difference between the TMR with and 
without accessories. Therefore, as the 
distance to the source of contamination 
decreases, the dose increases. The TMR 
for SPD = 75 cm (Fig. 5) is greater than 
that for SPD = 100 cm (Fig. 4) at all 
depths and for all accessory conditions. 


Mackie and Scrimger 


VR: vai : Fert tnabee kr peyacied ‘e 


eb p= ha yma 7 pi nee eeardionte waa : 
2 ‘ 7 ; i 
. : 
P L le + 5 tet ncn aa at 
: i ( ’ abe a0) 
. td: wach al AE mseret Vnont val 
} | week ety, é : 
i so a ‘- ia Annee | 6 
an aie +) tale writen fe or oe 
. 1 CG, asta ielaeorty, pen 
’ te 1 400 AO CoN a Coed =(! ‘ me 
Aart] 
Suiveowlly: | 
7 > oof @ th erys Tf 
iy gg! CD Ad eam MHA. 
F ng Sr ak Bi mite y, ny aas 
yi, Aate iaey Gi le Sina é 7 ioe 
ro) 4 a) Vor Nand eg 2 penueds J ietias 
P % Pid ily Wh wiriele iu} in iy ory rytert i ‘y ‘ae 
0%. 547 aw nie egy Pid id A garam =P 
' A Le ee Oe ee) ae eee 
ete 's 9 rele a 
a7 ToLutil Pie ai nee? ww 712le 
ts i Seto py Av) Vata, 
: ' 4 a; ¥ pir 
Pe Mt eed ses eee 
: aa eA ; | | hd as 
. es hres pacaneney : ay 
rcarted ‘phic } i) ow 
chee gitucl 4p ap cline a af; sant a Pe 
‘ya vabeie (2th 
We nl aaa oe as “a My 
- some ee 
‘ 7 
r ve, ae tele a) 
i, he nen teat ef 
: a Weal bine A a 
‘ "tice 
- ° b| i 
el dae) 
4 ; 
Gy >. 
grt he 
he Pree / 
ee Pe) * i 


Ts Geet: A yey & Te a 


Pee re) ry at ial 
VeVi) Ce La en 
ae en ee ee 
. nays dian 


apt nae . 


(al. five 
Anis haemo ty) ria 


The difference between the TMR asa 
function of SPD is most pronounced at 
the surface. 


Discussion 


As the thickness of accessory in- 
creases, the amount of contamination 
it produces increases. However, con- 
tamination produced at shallow depths 
within the accessory is shielded by the 
deeper layers. At a thickness equiva- 
lent to 3.2 cm (1! in.) of Lucite, there 
is a balance between production and 
absorption of contamination. 

These observations agree qualita- 
tively with those of Velkley ef al. (5) for 
25 MV x-rays collimated to a 10-cm X 
10-cm field and Gray (1) for 4 MV x- 
rays. 


Ill. DISTANCE 
DISTRIBUTION OF SURFACE 
DOSE 


Materials and Methods 


The surface maximum ratio is de- 
fined as the TMR measured with the 
chamber at the surface (i.e. d = 0). 
Whenat the surface, the chamber used 
has a sufficiently thin window to en- 
sure that the build-up depth is ap- 
proximately zero (5 um polystyrene 
equivalent). The window thickness has 
an equivalent mass stopping power of 
less than one centimeter of air. Only 
very low energy photons have a rea- 
sonable probability of interacting with 
such a window. In order to be detected, 
such photons have to be produced 
within a few centimeters of the cham- 
ber, otherwise they will be rapidly at- 
tenuated in air. Primary photons of low 
enough energy to interact with the 
window are completely attenuated 
when the beam emerges from the beam 
defining head. Therefore, most of the 
dose at the surface is attributable to 
contaminant electrons. 


Results 


Figure 6 illustrates that the surface 
maximum ratio in the open beam in- 
creases rapidly as the SPD decreases for 
both 20-cm X 20-cm and 30-cm X 30-cm 
fields. 

The curves fora collimated field size 
of 20 cm X 20 cm show the effect ofa 
thin accessory. The Siemens Mevatron 
20 is equipped with a removable, l-mm 
thick, Lucite cross-hair tray which can 
be inserted at 40 cm from the source. 
When the SPD is less than 75 cm, the 
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Distribution of the surface dose as a function of source-probe 


distance. TMR = tissue maximum ratio 


surface dose with the cross-hairs in 
place is greater than when it is re- 
moved. However, the cross-hairs, 
when in place, reduce the surfacé dose 
in the clinical region of SPD greater 
than 75cm. 


Discussion 


The reduction in dose when the 
cross-hairs are in place is a notable ex- 
ception to previous observations that 
accessories in the field increased the 
surface dose. If it is assumed that all the 
electrons set in motion in the thin Lu- 
cite accessory by the incident photon 
beam emerge, they will be fewer than 
the number of incident electrons 
stopped. This is because the probability 
of production of electrons increases 
with increasing slab thickness whereas 
the survival of electrons already pro- 
duced decreases with increasing 
thickness. 

As the SPD increases, the volume of 
air between the target and probe in- 
creases. For ©Co, Nilsson and Brahme 
(8) predicted an increase in dose as a 
function of SPD due to ®°Co photon 
interactions with air. In these experi- 
ments, the reverse appears to be the 
case and a decrease in dose as a func- 
tion of increasing SPD has been ob- 
served. This implies that at 15 MV the 
interaction of the primary photon 
beam with air is not the major source of 
electron contamination. 

The surface maximum ratio is nor- 
malized to the dose measured at diay 
within a polystyrene phantom. The 
maximum dose in a phantom is due 
almost entirely to primary photons for 
which the variation with distance falls 


off as 1/(SPD)*. However, the electron 
contamination does not necessarily 
arise at the source of primary photons. 
Instead of normalizing the surface dose 
ata given SPD to the dose at diay, each 
surface dose at the same SPD is nor- 
malized to the dose measured at a 
convenient point inside a phantom 
using a standard field size and photon 
SPD. The normalized dose is defined 
as: 


Dose to point at arbitrary d, 
SPD, field size, 
Dose to standard point 


ND 


[2] 


The standard point chosen was at the 
dma Pont ator Di 100 cm and field 
size = 10cm X 10 cm. 

The apparent source of contamina- 
tion electrons can be located assuming 
that it isa point source and there is no 
attenuation of electrons by air. There- 
fore, the normalized dose at the surface 
is directly proportional to the inverse 
square of the distance to the source of 
the electron contamination, or: 


ND(d = 0) [3] 


ot 

(SPD — do)? 
where do is the distance of the source of 
contamination electrons below the 
primary photon source. 

Figure 7 is a graph of the inverse 
square root of the normalized surface 
dose versus SPD. The graph yields a 
straight line except for very large SPDs. 
Therefore, the source must be small 
and the inverse square approximation 
valid for SPDs less than 120 cm. The 
y-intercept of the graph yields a dis- 
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The Gaussian distribution of the surface dose. The inset is a 
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Graph illustrating the inverse square depen- 
dence of the surface dose (ND) on the distance 


to the source of contamination (SPD). 


tance do of 12 cm below the primary 
source of photons. This corresponds 
approximately to the position of the 
bottom of the fixed head assembly 
which consists of the field flattening 
filter and the beam monitor. The bot- 
tom surface of the field flattening filter 
is 3.1 cm wide; at distances greater than 
50 cm this subtends an angle no larger 
than 3.5°, which approximates a point 
source. 

The field flattening filter has a 
maximum thickness of 5.4 cm of stain- 
less steel, which represents the largest 
interaction cross-section of any com- 
ponent in the beam defining head. The 
beam monitor is the last major compo- 
nent completely in the beam before the 
phantom, so it is likely to contribute to 
the dose due to contaminant elec- 
trons. 


IV. GAUSSIAN LATERAL 
DISTRIBUTION OF SURFACE 
DOSE 


Materials and Methods 


The magnitude of surface dose as a 
function of lateral displacement from 
the central axis was measured for sev- 
eral SPDs. The central axis angle (CAA) 
was defined as the angle between the 
central axis and a line joining the probe 
to the apparent source of electrons. 


July 1982 


This established parameters for the 
lateral direction in order to compare 
the distribution at various SPDs (Fig. 8, 
inset). In each case, the surface dose at 
any CAA and SPD is normalized to the 
dose at dyax on the central axis (CAA = 
0). 


Results 


The curves in Figure 8 approximate 
a Gaussian distribution. At all points 
across the field, a smaller SPD results 
in a larger electron contamination 


dose. 


Discussion 


The surface dependence measured 
along the central axis is in agreement 
with the finding at all points across the 
field that a larger electron contamina- 
tion dose is the result of a smaller SPD. 
Both the lateral distribution and the 
dependence on the distance from the 
source of the electron contamination 
agree with the work of Almond (12). 


V. BUILD-UP CURVES 
OUTSIDE THE FIELD 
Materials and Methods 


The primary photon dose outside 
the field boundary 1s very small (ap- 


Volume 144, Number 2 


schematic diagram of the experimental arrangement. 


proximately 1% of the central axis dose 
at dmay). Contamination electrons and 
scattered photons provide the largest 
contribution to the dose outside the 
field. This was studied experimentally. 
In order to get as close as possible to the 
field boundary, the probe was placed 
near the edge of the phantom (Fig. 9). 
The phantom was tilted at an angle? so 
that its surface would not move closer 
to or farther from the field boundary 
with the addition of build-up layers. 
This also reduced the number of par- 
ticles incident at oblique angles upon 
the phantom. 


Figure 9 
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Diagram of the experimental arrangement to 
determine the build-up curve outside the 
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Build-up curve obtained 3 cm outside of a 
30-cm X 30-cm field (see Fig. 9). 


Results 


Figure 10 illustrates the depth dose 
curve 3 cm outside a 30-cm X 30-cm 
field. The TMR has a peak at a shallow 
depth (1-2 mm), then falls off rapidly 
at depths beyond about 2 cm. 


Discussion 


The general shape of the curve 
agrees with the work of one of the au- 
thors (10) for 8-MV x-rays. This curve 
represents the contamination depth- 
dose due to electrons and scattered 
photons beyond the field edge. The 
peak at shallow depth followed by a 
rapid fall-off in dose is typical of the 
depth-dose curve for electrons. The 
relatively slow fall-off at depths be- 
yond 2 cm is due in large part to scat- 
tered photons. 


VI. LINEARITY OF DOSE vs. 
FIELD WIDTH 


Materials and Methods 


The dose as a function of field 
width was measured for a number of 
depths with the limited polystyrene 
phantom of the type illustrated in 
Figure 2. Since there was unit density 
material on top of the ion chamber, it 
could not be assumed that all of the 
contamination contributing to the dose 
was due to electrons. 

The phantom was placed completely 
within the field so an increase in de- 
tected signal due to increasing field 
size could only be attributed to beam 
contamination. To obtain the dose as a 
function of depth for as many field 
widths as possible, the phantom lateral 
dimensions were kept as small as pos- 
sible (i.¢., 10 cm X 10cm), which in turn 
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Linearity of the surface dose with respect to field width under 


a variety of conditions. 


dictated a minimum field size used in 
the determination of the contamina- 
tion depth-dose curve of 15 cm X 15 
cm. The maximum square field size 
attainable for the machine studied was 
30 cm X 30 cm. These field sizes de- 
fined the range of the study. 


Results 


The normalized dose (as defined by 
Equation 2) for d = 0 is plotted as a 
function of field width in Figure 11. 
Under a wide range of conditions there 
is a linear increase in the surface dose 
with respect to field width. Note that 
either decreasing the SPD or placing 
accessories in the field increases the 
slope. Since the surface dose for any 
given field size is proportional to the 
slope of the line, the slope is also a 
measure of the electron contamination, 
but with the advantage that it is not 
dependent on the field width. 

Figure 12 illustrates that for various 
depths less than djyjqy there is a linear 
dependence of dose on field width. 
Qualitatively, the slope is greater at 
shallow depths than at deeper depths. 
As has already been seen for the dose 
at the surface, the slope of normalized 
dose versus field width is a measure of 
the magnitude of contamination pen- 
etrating to the depth specified. There- 
fore, there is a greater magnitude of 
contamination at shallow depths. 


Discussion 


A number of authors (3,5, 6,8) have 
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The increase in the normalized dose with re- 
spect to field width for various depths within 
the phantom. The increase can only be attrib- 
uted to contamination because the phantom 
(lateral dimensions = 10 cm X 10 cm) was 
completely contained within the field 
boundary. 


commented on the dependence of the 
surface dose on the square root of field 
area both with and without accessories. 
Velkley et al. (5) normalized the surface 
dose at each field size to the dose at dna, 
for the same field size. However, since 
the dose at dja, increases with field 
size, such a plot would not unambigu- 
ously express the dependence of sur- 
face dose on field size. The concept of 
normalized dose, however, eliminates 
the ambiguity because the dose is nor- 
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malized to a quantity that does not 
change as a function of field size. 

The magnitude of contamination has 
been shown to be proportional to the 
slope of normalized dose plotted 
against field width. To obtain a 
depth-dose profile for the contaminant 
contribution, the surface slope is set 
equal to the surface dose. This estab- 
lishes a conversion between the slope 
at some depth and the dose due to 
contamination at the same depth. An 
assumption is made, as was previously 
discussed, that the dose measured by a 
thin-window ion chamber without any 
build-up is due to contamination elec- 
trons. 

The primary photon contribution is 
obtained by subtracting the contami- 
nation component from the total 
build-up curve. Figure 13 shows the 
separate contribution from the con- 
tamination and primary photon com- 
ponents of the total build-up curve for 
a 20-cm X 20-cm field. The contami- 
nant contribution curve dominates the 
exposure dose at depths less than about 
1 mm. At depths near day, the con- 
tamination contribution is less than 5% 
of the primary photon component. 

If the contaminant contribution is 
renormalized to its own maximum, the 
depth-dose profile becomes more evi- 
dent. Figure 14 shows the renormal- 
ized contaminant contribution curve. 
The peak in this depth dose curve oc- 
curs at about 1-2 mm. There is a rapid 
fall-off at depths slightly greater than 
the day. At depths greater than a cou- 
ple of centimeters, the curve exhibits a 
slow fall-off with depth. This curve 1s 
similar qualitatively to the build-up 
curve measured outside the field (Fig. 
10). 

Since the contamination contribu- 
tion was obtained by attributing the 
increase in dose to an increase in field 
size, the primary photon component 
should be independent of field size. 
TABLE I confirms this. Thus, the mag- 
nitude of contamination depends on 
the field size but the spectral charac- 
teristics of the contamination do not. 

The derived primary photon TMRs 
in TABLE I are independent of field 
size. TMR was not, however, designed 
to be independent of field size for 
conditions other than an open field 
and SPD = 100 cm. The primary pho- 
ton contribution is independent of SPD 
and the presence of accessories. The 
primary photon build-up curves for 
various conditions of SPD and acces- 
sory states are shown 1n TABLE Il. The 
conditions of SPD were tested for 75 
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The contributions of the contaminant and 
primary photon components to the total 
build-up curve fora field size of 20 cm X 20 
cm. 
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The contaminant and primary photon com- 
ponents in Figure 13 renormalized to their 
separate maxima. 


TABLEI: Primary Photon Tissue Maximum Ratios for Various Field Sizes (No 
Accessories, Source-Probe Distance = 100 cm) 


Tissue Maximum Ratio for Field Sizes 


No. of Depth 5 cm .,  kOvem 15cm 20 cm 25 cm 
Slabs (cm) x 5cm x 10cm xX 15cm x 20 cm 25 em 
0 0.00 0.00 0.00 0.00 0.00 0.00 
1 0.16 0.19 0.19 0.18 0.19 0.18 
2 Ons2 Om38 0.33 0.33 0.33 0.33 
3 0.48 0.44 0.44 0.44 0.44 0.44 
4 0.64 0.54 0.54 0.54 0.54 0.54 
6 0.95 0.68 0.68 0.69 0.69 0.69 
8 127, 0.79 0.78 0.79 0.79 0.79 
12 OH 0.91 0.91 0.91 0.91 0.91 
16 2.54 0.97 0.97 0.97 0.97 0.97 
20 3.18 1.00 0.99 0.99 0.99 0.99 
24 3.81 1.00 1.00 1.00 1.00 1.00 
TABLE II: Primary Photon Tissue Maximum Ratios under Various Conditions 
(Field Size = 30 cm X 30 cm) 
Tissue Maximum Ratios 
SPD = SPD = SPD = SPD = SPD = 
No.of Depth 75cm 100 cm 140 cm 100 cm 100 cm 
Slabs (cm) No Access. No Access. No Access. 0.64-cm LTA*  3.2-cm LTA* 
0 0.00 0.00 0.00 0.00 0.00 0.00 
1 0.16 0.19 0.19 0.18 0.19 0.18 
2 Ons? 0.34 0.33 O82 O38 0.33 
8 0.48 0.45 0.44 0.43 0.45 0.45 
4 0.64 0.54 0.54 0.53 0.56 0.55 
6 0.95 0.68 0.68 0.67 0.70 0.69 
8 127 0.78 0.78 0.77 0.80 0.79 
12 1.91 0.90 0.90 0.90 0.91 0.91 
16 2.54 0.97 0.96 0.96 0.97 0.98 
20 3.18 0.99 0.99 0.99 1.00 1.00 
24 3.81 1.00 1.00 1.00 1.00 1.00 


* LTA = Lucite tray accessory 


cm, 100 cm, and 140 cm, whereas the 
electron contamination curve Was de- 
rived for SPD = 100 cm. The indepen- 
dence of the primary photon build-up 
curve for these wide-ranging condi- 
tions illustrates that contamination 
accounts for the difference in depth- 
dose curves encountered. By defini- 
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tion, the primary photon contribution 
to the total build-up dose is due to that 
portion of the beam that does not in- 
teract before reaching the phantom. 
Changes in the primary photon 
build-up could only be due to changes 
in the phantom size, density, or con- 
figuration. 
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CONCLUSIONS 


The contamination depth-dose curve 
(Fig. 14) and the build-up curve outside 
the field boundary (Fig. 10) have sim- 
ilar shapes. The surface dose has been 
established as being due to contami- 
nant electrons. The contaminant elec- 
trons evident at the surface must pen- 
etrate to some depth beneath the 
phantom. The peaks and rapid fall-off 
at shallow depths illustrated in Figures 
10 and 14 are characteristic of electrons. 
The slow fall-off at depths greater than 
2 cm is probably due to scattered pho- 
ton contamination. Thus, the contam- 
ination radiation appears to be a com- 
bination of an electron contribution 
primarily at shallow depths and pho- 
tons contributing primarily at greater 
depths. 

Figures 10 and 14 differ in shape due 
to a difference between the relative 
amount of contamination electrons and 
photons present inside and outside the 
field boundary. Figure 8 demonstrated 
a Gaussian lateral distribution of the 
surface dose. Consequently for a 30-cm 
x 30-cm field there were about 45% 
fewer electrons present at 18 cm from 
the central axis than were present at 
the central axis. 

Although the contaminant electrons 
are produced by the fixed components 
in the beam defining head, the size of 
the collimator opening dictates the 
relative amount of electron contami- 
nation emerging from the beam de- 
fining head. Low energy electrons are 
scattered by air through wide angles. 
The product of mass angular scattering 
power, density, a nd the electron path 
length for low energy electrons in air 
is of the same order of magnitude as the 
central axis angle.* Therefore, the angle 
of an electron when emerging through 
the collimators is not critical in deter- 
mining whether it will arrive at the 
detector. The greater the collimator 
opening, the greater the flux of con- 
tamination electrons arriving at the 
detector. It has been shown empirically 
that the magnitude of electron con- 
tamination reaching the detector at the 
central axis is directly proportional to 
field width for square fields. 

The contamination attributable to 
photons was derived from the increase 
in the dose as a function of field width. 


3 For example Rossi (13) states that a 3-MeV 
plectron in air has a mass angular scattering 
power of 0.602 radians2-cm2/g which, after pas- 
sage through 100 cm of air, corresponds to a mean 
scattering angle of 16.0". 
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Therefore, the amount of material with 
which the primary photons could have 
interacted to produce scattered photons 
must be dependent on field size. The 
surface area of the collimator faces and 
the amount of air in the field both in- 
crease with field size. The contamina- 
tion dose due to photons and the sur- 
face area of the collimator faces are 
both approximately proportional to 
field width, so it is likely the contami- 
nation photons are photons scattered 
by interaction of the primary photon 
beam with the collimators. 

The shift in day, especially when 
accessories are present in the beam, can 
be attributed to a greater magnitude of 
electron contamination. Since the 
depth-dose profile for the scattered 
photons is likely to be similar to the 
primary photon profile, they are un- 
likely to contribute significantly to the 
shift in day. The peak in the contami- 
nation curve due to electrons occurs at 
very shallow depths, so increasing the 
magnitude of electron contamination 
results in a dma, Shift towards shallower 
depths. This also explains why the shift 
is accompanied by an increased surface 
dose. These findings agree with Biggs 
and Ling (4) who observed that a shift 
in dmax did not occur when electrons 
were swept from the beam. 
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